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Sommaire. — Le but de ce travail a été de calculer la perméabilité apparente y, et le facteur de 
qualité Qx d’une poudre dont les grains magnétiques el conducteurs assimilés à des « particules » ou 
- «grains » sphériques se trouvent noyés dans un milieu isolant. L'ensemble est soumis à un champ magné- 
tique périadique et uniforme dont la fréquence peut atteindre plusieurs mégacycles. 

Le problème a consisté à calculer d’abord les champs résultants extérieurs et intérieurs à une seule 
particule. Puis, en adoptant l'hypothèse de Lorentz sur le champ efficace et uniforme, à appliquer 
les résultats obtenus à un milieu dans lequel la distribution des particules est supposée uniforme. Les 

E résultats donnent évidemment une perméabilité apparente complexe d’où l’on déduit aisément le 


coefficient de qualité Qx = — É et le facteur deréduction de flux n = Fr. Les caleuls sont longs et assez 
ù 


délicats, mais ils permettent d'obtenir des solutions, fonctions de la fréquence, de la perméabilité, de 
la conductibilité et de la grosseur du « gain », en très bon accord avec les résultats expérimentaux. 


Considérons une sphère de conductibilité y, de 
perméabilité y: et de pouvoir inducteur spécifique & 


placée dans un milieu isolant (y; = 0, p4, 1). Sur 


cette sphère agit un champ magnétique uniforme 
et sinusoiïdal de pulsation « telle que la longueur 
d'onde des oscillations correspondantes 


(e = 3.101 cm:s) 


soit grande par rapport au diamètre de la sphère. 
Nous nous proposons tout d’abord de calculer 
les valeurs du champ électrique et du champ magné- 
tique à l’intérieur et à l’extérieur de la sphère. 
Pour nous permettre d’expliciter les conditions 
aux limites, nous adopterons les coordonnées sphé- 
riques ayant pour centre celui de la sphère, le champ 


inducteur à — H 0! étant dirigé suivant l’axe 
des z (fig. 1). 
Nous appellerons 


> ; 7 à 
3€ — H ei”! eb b —)}eiwt 


le champ magnétique et le champ électrique (4 et R 
sont indépendants du temps). Les inductions corres- 
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pondantes seront 


+. > > > 
GB = LÜe et bre be 

Nous supposerons que les quantités y et e sont 
constantes et le champ inducteur assez faible pour 
éviter la saturation si la sphère est ferromagnétique. 


A . 
Si o est la densité de charge et à la densité de 
courant, les équations de l’électromagnétisme qui 
s’écrivent. en U.E:M:: 


He 
10h 
rotb =— — i 
dt? (1) 
La 
F2 > 1 db > > 
rotdt = (ET or avecet.—= Yh; (2) 
c? dt 
ES 
div = 0, (3) 
Mie 4 
divé = 4rcp (4). 


se réduisent évidemment aux seules équations (1) 
et (2), si l’on suppose les solutions harmoniques 
et si l’on tient compte de la condition de conser- 


L4, 


- A 
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vation du courant 
| % 


divi=— —+ 


ot 


D'autre part, les conditions de symétrie per- 
mettent d'écrire 


= 


: < 5e 
br = 0, bp = 0, En = O. 


Ho 


Fig. 1. 


Les équations (1) et (2) se réduisent alors aux 
suivantes 


I d 


r sinÿ mp Ce sm0) +/ou/,= 0, . (3) 
= à L(rhy) + J'ou H6— 0, F 1 (4) 

ee 10H RTE JE ie 5 
r RUE r OÙ (: CAE c? }hr= 0: (5) 


Les équations (3) et (4) donnent H, et H5 en 
fonction de h;. En tenant compte de cela (5) se 
transforme en une équation (8) à laquelle doit satis- 
faire le seul champ A: 


I ‘0 I 
A Ans) et ) 
11, Fan0 (hs sin me a (6) 
1 0 I 
| ME Tor ho) co 


1 9 Fo(rhe) 1 0 I dy n0 
mor or | Tr 00 | 7600 06 Ve M0) 


$ DEN... 
_ (54 FUYUW + ee Je o. (8) 


Nous poserons 


PSE 
a /'O2UE Ë À 
QE (VA — — jJATUUY, (0) 


ec? 


et nous nous placerons dans le cas du régime quasi 
stationnaire, c’est-à-dire que nous négligerons le 


courant de déplacement devant celui de conductior. 
Comme nous l'avons déjà vu, cela est parfaitement 
admissible jusqu'aux fréquences les plus élevées 
(f = 1011 par exemple, c’est-à-dire : — 3 mm), N’ou- 
blions pas cependant que la longueur d’onde corres- 


pondant à ces fréquences doit être nettement plus … 


grande que le diamètre 24 de la sphère, de façon 
à pouvoir définir une région dans laquelle on puisse 
admettre l’uniformité du champ [1]. Ainsi nos 
exemples porteront sur les diamètres de l’ordre de 
quelques microns qui sont utilisés dans la technique 
haute fréquence. | 

Remarquons qu’à l'extérieur de la sphère on a 


& =, E — €). 
» 


Donc 


c; est la vitesse de la lumière dans le milieu extérieur. 


Posons également 


rhy sin0=S. 
L'équation (8) devient 


a + AS — 0. 


(10) 


Suivant la méthode classique, remplaçons S par 
le produit des fonctions R,,.et 0, qui ne dépendent 
chacune que de r ou que de 0 et divisons tous les 
termes obtenus par R,, 05. On obtient 


1 d'A 
R-dr2 


1 1 d?0 TT 
0 r?.d0? | 0 72 


ae 


Comme on sait, cette équation se décompose en 
deux équations qui sont 


d20 dû 

LATE CARE L=- 19 
TE cotg + L0—=o, LA 
dR Nes 3 
a +(e- ms) = 0. (13) 


Il est bien connu que la première équation n’a de 
solution régulière en 0 que si 


LE MUNIE GED TE 


Si L a l’une de ces valeurs, une discussion assez 
simple de‘cette équation montre que sa solution 
qui est liée aux polynomes associés de Legendre est 


0/=— sin0P}(cos0). (14) 

La seconde équation ; 
arf 9 ]r = 0 (15) 

dr? r? 


peut se transformer en posant 
1 


R=Yÿ (DT? el 


ESS Ge 


Job #6 


en une “équation de Bessel 


1 df 


u du 


def. 


du? 


La seule solution de cette équation d'ordre L + = 


qui reste régulière pour r — o, est J,.:!(Kr). La solu- 


tion de l’équation (15) est donc 


1 
R— MrèT, \(Kr) (AT étant une constante). 
[+3 5 
La solution générale à l'intérieur de la sphère 
est donc 


S; = — Mr2J 
L+ 


\(Kr)sin0P}(cos0), (16) 


avec 
S=7rh; sin0. 


A l'extérieur de la sphère, si l’on néglige le courant 


de déplacement, on à K—o. Donc (15) devient - 


d? À 


Her), 2 


dr? r? 


qui admet pour solution 


"2 


Le : 


2 


Gs 
R = Birti+ Re 


A l’extérieur de la sphère, la solution générale 
est donc 
Se=— (frise À À) sin0 P} (eust), (17) 
avec 
S = rh sin. 


Remarquons que là continuité de hk; exige d’autre 
part la même valeur pour ! dans (16) et (17). 


Calcul de / et de B,. — A l'infini on doit trouver 


. pour À, et A5 des valeurs indépendantes de r qui 
sont respectivement 


H,.= H, cos0 et Hy=— H, sin0. (18) 


Or (6) et (7) donnent, avec (17), 
2S 
r2 sin0/ou; 00 


I PRE e 
= ——— |Gr he 
Joux sin 0 : 


He Her 


| Line PI (cost) 
I 2S 
r Sinüj7ou or 
— P}(cost  fo2 
LÉ) entres |: 
" 2 rl? 


- Jop 
Pour que H, et H5 soient indépendants de r à 
l'infini, il faut done 1 — x. 


Dans ces conditions, comme P} (cos 0) — 


Ho = 


— sin 0, 
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| on déduit finalement 


on en tire 


Se (ar =) sin? 0. 


Enfin, avec: l —5r et par l'intermédiaire de la 
condition à l'infini (18), on déduit 


PET 0 4 
Dit SNS 


On a donc finalement les solutions suivantes : 
A l'intérieur de la sphère : 

L 

Si= Mr°J,(Kr)sin?0. 


A l'extérieur de la sphère : 


cos À | s 


2 sin Æ 7° L 
Ve M 20S z LAS 110" 
Si= M4/ | Kr cos Kr | sin 0, (19) 
Go 
SEE (er: 2 È) sin? 0, (20) 


Calcul des champs. — Dans la sphère l’équa- 
tion (19) nous donne 


S2 nr. DE Ed sin X r 
r sin 0 TA ETE Kr 


Puis par (6) et (7) 
ER kr| ) 


VE 
) 


— JM ee Sin 0 

LL) \ KZ Le 
L. cos À r I 

A l'extérieur de la sphère des calculs identiques 

portant sur S, donnent 


He 


— cos r| 


sin À 7° 


E 27 
He y Tr 


0 


(21) 


re 
Kr + sin À7 ess 7 


e Se Hi: ; Bo in 0 
n = = = | — D 6) 7 + “— }- sin 
si: r sin Ü ee. Ju r? î 
2 7 f9 0 99 
H,=  Hcos0 + ty cos D, (22) 
ee] 
= TE n0 


5 sin 0 + 
es 
Posons dès maintenant 


et LE —= BH. 


M \ / 2 =N HS 


Déterminations des constantes N et B. 
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L'expression de la continuité des composantes 
tangentielles À; et Ho à l’intérieur et à l’extérieur 
de la sphère de rayon a fournit les deux RU 
suivantes : 


Continuité de À, : 


HiN(sinKa—KacosKa)= BHK— Ji so, 


DS 


NS " Lu © a K 


. Li 
Bi jÉvalieshr 


Utilisation des fonctions de Bessel de 


et =. — En partant des relations suivantes : 
» à 4 à n° & 
a is à 
. DE 
sin À7° 4 rANNE 
K7 — cos À7 ={ : ) Ja AC 


3 A5 À 
== 1) sin À 7: — = SK r = pra 


3 
kK? 7? 


on obtient pour les composantes des champs inté- 


rieurs à la sphère Her 


1 
2. 


ho=— ne sin 0 7r Ji(Kr) + sr 


(25) 


1 


He Treo] Ja( Kr)+J; sk); 
É 1 ï FES 
Hy = CT sinon iv, Go LU : 
avec k 
VENT INK £ à TEA 
T= Fe CFA 


Chaque fois que cela sera possible, nous suppri- 
merons l'indice (Xr) ou (Ka) des fonctions J afin 
d’alléger l'écriture. Nous aurons donc pour 4 N, 


B les valeurs suivantes : SES 

j ; Le 

A 31 a? ES 
T= 1 1(28 
CERN PEN TETE ES nn” 

+ 4: 2 à ee 

97 Lui © a? ï ; 
IN re (29) 

KG3KY SE DAS fee TE 
B = Ji 0 & 2 Cu + ET (30) 

2 (2m + p)Ji + (up) à 
2 2 


Remarquons dès maintenant que pour des valeurs 
de X très faibles on peut utiliser les dectpents 
en série de AS : 


Atkr)r (2 =) =) à 

1 AU ! 

J;(Kr) w (=) = ù 5 

Ê 2 Re | 
RE 


Ke @ ta smKa+ (D — m){: snka-—kKa cos À « ] é 
csinKa—KacosKa à 
Ÿ era sin À @ + (u — Bo PTE cosKa]\ 


RS UE HN @—1)snKa+ Kacoska). 
= — BA Ur 
Ces deux équations sont sufisantes s pour déte 


miner “é 2 B. On trouve Me 


# 


. Ce qui donne 


Te A Su = 
2K \? 
| CE) 
êt < 
| # 
= DRE ES JUi7e 
De NN Ale rTenre 
PA M tt Ur 
D te 
2h41 + 4 
— 3} è er : 
H5 = pau 
de SES AU 
ee F2 PCR 


valeur 


H:=— 7 sin “ Hrcosû Æ Du. 
‘ crie Le 
Pour = 1: dindicbon Re donc 
HS 5 30 _X 
si 3 De f 
| Fa . je HS +2 . 


une sphère isolée de Dana à Fe: 


GA d6 RE nte —— re D. 
résultante suivant Oz des champs Het Ho: - 


H:=— Ho sin 0 + H,cos0 LES 
&- 4 ASS = = ? : « % 
= HTree A) (ere sn ete 


LT? aimantation Correspondante suivant (2 e 
par unité de volume £ 


rue Hs. dv. 


L’ aimantation pour toute + sphère äe rayon 
est donc 


TE et 1) 


47 e ÿ “ 
a ñ & ‘4 s 
< $ [ rinanr- . sine | +2, (kr). 
0 02180 20° CMêe 
ce x r? sin0 dr dde 


eo 
soit ; 
fee & à 


rJ,(Kr) dr 


3 7 4 1 FE 
» ee | kr | 


0 


qui, en tenant compte de la relation 


« 


ê 7 px AE Ty 
MERE Sr (Kr) + J:CKr)] 


3 


donne sans difficultés 


5 e : 
Er Hu ss 1) GE JACK &), 
K | Eà 
En posant 
M —1) a J,(K a NS. 


on aura donc 
G 4H, 


le coefficient « caractérisant à la de de consi- 


dérée le degré d’aimantation que la sphère est 


susceptible d’acquérir. 


Application. Considérons un milieu isolant 
de perméabilité unité dans lequel sont noyées 
un grand nombre de granules ou particules sphé- 
riques semblables à la sphère déjà étudiée. Leur 
densité de répartition est supposée constante et leur 
diamètre identique pourra être compris entre une 
fraction de micron et quelques dizaines de microns. 

Il est bien évident .que si les distances mutuelles 
des particules sont extrêmement grandes, l’aiman- 


tation totale est la somme des aimantations indi- 


viduelles. . 

Le problème est différent, si, les distances sont 
petites. Il faut alors tenir compte des influences 
qu’exercent sur la particule étudiée toutes les autres 
particules de la substance qui sont, elles aussi, 
polarisées. x 

Or, le champ réel dit efficace F qui agit sur une 
particule est celui qui existerait en ce lieu si l’on 
supposait que cette particule en fut enlevée sans 
troubler l’état des autres particules. Si nous considé- 
rons l’ensemble de la substance, hors de la particule, 
comme aimanlée d’une façon homogène 5, 
n’est qu’une première approximation, puisque les 
_ particules sont en nombre fini quoique très grand; 
la quantité de matière extérieure à la particule 
sphérique agit comme si sa superficie portait une 
charge de densité —#, (n normale extérieure). 

Un calcul classique donne pour le champ engendré 
à l’intérieur de la PAGE par une telle distribu- 


tion Az à, le coefficient ‘. 7 est le coefficient de Lorentz. 


0 
Nous admettrons de comme valable pour notre 
cas le raisonnement utilisé pour Becker sur les 
diélectriques : aussi longtemps qu’on peut assimiler 


# 
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(4). 


ce qui- 
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les effets mutuels des particules à ceux de dipôles 
ponctuels, le champ efficace qui agit sur une particule 
individuelle est donné par 


Fer 5-0). t 
0} 

Tant que la répartition des champs ne sera pas 
perturbée par des contacts entre particules résultant 
d’une densité de répartition trop élevée, cette hypo- 
thèse sera d’autant plus justifiée que le diamètre 
des particules sera plus faible et leur nombre par 


- unité de volume plus grand. 


Calcul de la perméabilité apparente et des 
facteurs de qualité et de réduction de îlux. — 
Appelons N le nombré de particules par centimètre 
cube, a leur rayon, 7 le volume total qu’elles occupent 
dans r cm. On a 

37 


LT 
ce 


NE SL (z varie de 0 à 1). 


D'après le raisonnement précédent, nous consi- 
dérerons F au lieu de H,. Donc 


L’induction apparente est donc 


Ban— Ho (TS — Hi 


Et la perméabilité apparente 


ME E Ba D 
« 
. Ho 


a+ 274 


ad? — Ta 


Remplaçons « par sa valeur tirée de (28) et (31) 
en remarquant que 


Ja(Ka)= SEA CK&) + J3 EAT 


_ En supprimant l'indice (Ka) pour no écriture 
on déduit 
(u—1)(1+27)J,; 


LE 
£ £: (32) 
È 


Eur) GT): 


[u+2+92r(h—1)] 


Ha — 


J 
[u+2—7(u—1)]J 
Nous pouvons dès maintenant remarquer que 
pour de très faibles valeurs de Ka (basse fréquence) 
les développements en série limités au second ordre 
des fonctions J1 (Ka).et J: (Ka) donnent pour la 


() à est l’aimantation présente dans la matière. Gelte 
formule s'applique rigoureusement pour le réseau régulier 
et en moyenne pour des corps amorphes dans le cas des 
diélectriques (8). 
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perméabilité apparente la valeur 


3+(u—71)(27t +1) 7e 
Tr ; X + (33) 


4 3+(m—1)(1—7) 


Elle est indépendante de Ka. 


En particulier, si 7 — 1 on retrouve bien p,, = 
CbISL = 001 A 0 —1ù1. Si est très grand avec 77, 
Lo, tend vers 


2% 


C'el 


Ho, = 


Si Ka n’est pas petit, z, présente une partie réelle 
et une partie imaginaire. 
Posons 


————— - 3 , te . # 
$ = V2ru0Yy a (ce qui revient à écrire Ka = s V2). 


J;(Ka) = ee | ( — 1) [sins coshs— 7 cosssinhs] 
— . (1+ 7)[cosscoshs +7 sins sinhs]. 
28 4 


En transportant ces valeurs dans l'expression 
générale de y, et en séparant les parties réelle et 
imaginaire on aboutit après un long calcuk aux 
expressions suivantes 


{ (u—1) (ir) (i+2r)[1+02s8R—05sP] |) 
l _ SU NS JO RE EN 
fe LH LER = ES PI 
lo +{+is(una-1)0 
= : I—Ss2 ne 
RER oruer ((u—1)(1—T}{1+2s2R—9sP| 0) 
l Æ 4St+ 255(u —1) (17) Q 


ÈS 


5N2(34) 


expressions dans lesquelles on a posé 


sinh2s + sin2s 


P — » 
cosh2$— cos2s 
sinh2s— sin2s 

O — s 

Ë cosh2s — cos2s 
cosh2$ + cos2s 

= 


cosh 25 — cos2s 


Les trois fonctions P, Q, R tendent chacune vers 
l'unité dès que la variable s est supérieure à quelques 
unités (pratiquement à partir de 2s — 4). 

Les formules (34) et (35) permettent déjà d’ob- 
server que pour des valeurs de s supérieures à 
quelques unités, si p est grand et + — 1 (remplis- 
sage complet), y, varie comme (\/u, ce qui est bien 
conforme aux résultats conuSs concernant le cylindre 
plein [1]. 

Sir;<1 pour les grandes valeurs de y, la per- 


27T 


, Er z : I 
méabilité 1; tend vers zéro, mais u, tend vers : 


ni T 


Ds 


N° 5. 


æ 


Ce résultat est tout à fait de l’ordre de grandeur 
des perméabilités apparentes observées dans la 


. . û LL 
pratique. Ainsi pour T—=-:, On a Ua—=4 et 
3 4 
pour 5 — >; on a p, — 10. Ce sont là des valeurs 
4 
courantes, 


Sur les formules donnant y. et u:; on peut également 
constater 


19 Que pour r = o(u 1) quel que soit s: p, = 
et ü— 0. 

Il en-estde mêmes =r( Tr) 

20 Que pour o<Tr<1, lorsque s croît, y; tend 


vers 1 et 4; tend vers zéro. 


Nous mettons ainsi 


en évidence un résultat 
important. 


Lorsque la fréquence augmente(s — \2rpwya), la. 


perméabilité apparente tend vers l'unité. 

Dans les mêmes conditions, nous avons vu [1] 
que la perméabilité apparente tend vers zéro pour 
le tube ou le cylindre plein. 

Les expressions (33) et (34) permettent d'évaluer 
les caractéristiques fondamentales d’une bobine 
dont le « noyau » est constitué par une poudre 
magnétique noyée dans'un isolant. 


Ce sont : 


10 Le facteur de réduction de flux : n — ra qui est 
le rapport du flux inductif appliqué à a Tu de 
même amplitude, mais de basse fréquence. C’est aussi 


le facteur de réduction de la perméabilité apparente. 
On a 


IUT 3+(u—1)(1— 7) 
œ IS PE) For) 
(u—x)(1—x)(1+2t)[r-H2sR—2sP]} 
4si+ 28(u —1)(2 +7)Q \ 
Gui QaS<} lire R— 2sP]] 
| VAE Au = T) (Er) O 


(36) 


Les développements en série de P, Q, R (jusqu'à 


la puissance 4) montrent que pour les faibles valeurs 
de s, n est très voisin de l’unité. 

Le facteur de réduction 1 ne commence à décroître 
que pour des valeurs de s beaucoup plus élevées 


que dans le cas du cylindre. D'ailleurs la limite de r . 3 


pour les grandes valeurs de s re en à 
u + 


dep et T. Ainsi pour r — _ A+ ? et vaut | = 


2U,+ 


pour les grandes valeurs de p (n=1 si p = à: 


Roues r — 0," + 1 et enfin pour + très voisin de &, 


n—> = = (zéro si est grand). 


Ces résultats indiquent une différence profonde 
avec le cas du cylindre plein ou creux. Essentiel- 
lement elle est due à ce que la perméabilité y inter- 


vient de façon différente (contrairement à « et 0 


A NS n 
Has SE 


WT : 4 


2 * 


dans les équations des champs qui concernent le 
_ Cylindre indéfini et la sphère. 


20 Le facteur de qualité 


Lo et R, sont l’inductance et la résistance de 
pertes par courants de Foucault [1] 


ee. 
\ —iŸ Gr) G+or is R—05P] 


{ + 4si+2si(u—1)(2+7T)Q 
% 37T(u—1)s2(sP— 1) RE) 


Br 


L'étude de cette expression montre que @;, 
très grand pour de faibles valeurs de s, décroît 
lorsque s augmente, passe par un minimum aux 
environs de 


ne le r)(1+ ONE (2) (38) 
ou en posant d — 2 4, 
uw d? = Ph re) (39) 


(2) La dérivée de Qr a quatre racines positives ce qui, 
conduit à deux maxima et deux minima. Mais seul le premier 
minima est à considérer. Au delà du second Qr croît comme s. 
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Le facteur de qualité reste donc important jusqu’à 
des fréquences d’autant plus élevées que le diamètre 
des particules est plus faible. Les grandeurs H, ©, *, 
r étant fixées, il y aura donc intérêt à utiliser des 
particules de diamètre nettement inférieur, à celui 
qui est donné par l'égalité (39). 

En fait la fréquence correspondant au minimum 
se trouve en général au-dessus du domaine des 
fréquences d'utilisation. à 

Suivant les usages auxquels on destine les noyaux, 
le diamètre des « particules » ou « grains » varie 
aux environs de 2 à 30. Ainsi la poudre de fer 
de 204 est utilisée pour des noyaux toriques jus- 
qu'à 100 kc. Celle de 34 est utilisée en très haute 
fréquence jusqu’à ro Mec (#). Mais la pratique montre 
bien que quelque soit la grosseur du grain, le facteur 
de qualité diminue rapidement au-dessus d’une 
certaine fréquence. Ceci est bien conforme à nos 
calculs et justifie nos approximations. 


Tous les calculs de cet exposé sont en très bon 
accord avec l’ensemble des résultats expérimentaux 
obtenus dans la technique des poudres. 


* (#) La valeur de Q, ne concerne que les pertes par courants 
de Foucault. Le facteur de qualité total doit tenir compte 
en plus des pertes hystérétiques et-des pertes résiduelles. 
Pour les ferrites, ce sont ces dernières les plus importantes. 
Mais il apparaît qu’on ne peut les réduire qu'aux dépens de 
la perméabilité [10]. 


Manuscrit reçu le 26 janvier 1950. 
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PROPRIÉTÉS MAGNÉTIQUES ET STRUCTURE DE LA PHASE QUADRATIQUE 
DU FERRITE DE CUIVRE | k 


. Par-L.-WEIEÆRSBERTAUT LS TE BOCHIROL. 


Laboratoire d’Électrostatique et de Physique du Métal, Grenoble. 


Sommaire, — On montre que le ferrite de cuivre prend par refroidissement lent, une structure quadra- 
tique à grand champ coercitif, à saturation magnétique plus faible et dont le rapport des 


(a » # # 4 GE Be 4 : 
axes atteint 1,06. À chaque température correspond un état d'équilibre caractérisé par une valeur 


Pac c : 
différente de rs comprise entre 1 et 1,06. 


L'étude aux rayons X, à l’aide du rayonnement Co X &, montre que les variétés cubiques et quadra- 
tiques sont proches de l’état de ferrite inverse, mais que la fraction À des ions Fe+*t croît sur les sites 
tétraédriques au fur et à mesure de la décroissance de la température. On propose une surstructure 
comportant une répartition ordonnée d'ions de Fe et Cu sur le sous-réseau des places octaédriques. 
Une telle surstructure, si elle existe, serait très difficile à mettre en évidence pour les autres ferrites, 


car les ions de Co, Ni, Mg, ete., déforment moins les réseaux que les ions de Cu**. La variation de FE 


aurait alors deux causes, d'une part, l'établissement d’un ordre entre les ions Gu*+ et Fe*t* dans le 
sous-réseau des sites octaédriques et, d'autre part, la variation de À, c’est-à-dire les changements de 
concentration des ions du sous-réseau : on FRIPR en effet que,pour toutes les surstructures, les paramètres 


varient en fonction du titre. 


1. La structure des ferrites. — Dans les ferrites, 
de formule générale Fe,0, MO, où M est un métal 
bivalent, les ions d'oxygène, beaucoup plus gros 
que les autres, forment un empilement à peu près 
cubique compact; les ions métalliques viennent 
s’insérer dans les trous de cet empilement. Des 
considérations de valence conduiraient à penser 
que les ions bivalents devraient se placer sur les sites 
tétraédriques avec quatre atomes d'oxygène voisins 
alors que les ions ferriques iraient sur les sites 
octaédriques avec six atomes d'oxygène voisins. On 
obtiendrait alors un ferrite normal [1, 13]. Mais 
d’autres facteurs doivent intervenir puisque, à part 
le ferrite de Zn, les ferrites ont toujours une fraction À 
des ions ferriques sur les sites tétraédriques À où 
ne reste donc qu’une fraction 1 —2À des places 
disponibles pour les ions Cu. Une fraction 2 des 
ions Cu se trouve sur les places octaédriques B en 
même temps que la fraction y — 1 — 2 des ions Fe 
restants. On dit que le ferrite est ‘inverse 
lorsque À = 0,5 

Nous avons pensé que c’est à des variations de À 
qu'il fallait attribuer les changements de propriétés 
magnétiques qu’on observe sur de nombreux ferrites 
en changeant les conditions de traitement thermique : 
le ferrite de cobalt par exemple présente une forte 
augmentation du champ coercitif lors d’un recuit 
prolongé un peu en dessous de 3500 C. Le ferrite 
de cuivre a un Champ coercitif beaucoup plus grand 
lorsqu'on le refroidit lentement que lorsqu'on le 


en, Ô Ze. 


trempe depuis 1 0000 C. En même temps sa structure 
se modifie profondément : il devient quadratique [2] 
au lieu de rester cubique comme à l'état trempé. 
De plus, examinés avec la raie Co K «, les ions de 
fer et de cuivre ont un pouvoir diffusant relati- 
vement différent, ce qui permet de les localiser 
dans la structure. C’est pour ces raisons que nous 
avons consacré cette étude au ferrite de cuivre. 

2. Préparation et caractères généraux des 
échantillons. — La préparation des ferrites a été 
faite soit suivant la méthode classique de mélange 
des oxydes suivi d’un traitement mécanique et 
thermique d’ homogénéisation, soit par la méthode 
mise au point à notre laboratoire [3, 4], utilisant 
la décomposition, suivie d’oxydation, de sels orga- 
niques mixtes. On obtient la phase cubique par 
trempe à partir d’une température supérieure 
à 8000 C; un refroidissement lent fait apparaître la 
phase quadratique. Nous avons vérifié aux rayons X 
la pureté des produits obtenus; il faut seulement 
éviter que, dans un séjour d’homogénéisation 
à 1 0000 trop prolongé, ne se produise une ségré- 
sation de CuO. 

Indépendamment de sa structure aux rayons X, 
structure relativement longue à déterminer, la 
phase quadratique est caractérisée par ses propriétés 
magnétiques : le champ coercitif est plus grand [5] 
que dans la phase cubique. Pour deux des échantillons 
il a par exemple passé de 12 et 6,4 à 4o5 et 190 Oe 


Je 


respectivement au cours d’un refroidissement lent. 
La perméabilité initiale, d'autre part, diminue 


de 45 à r,3. Enfin la saturation magnétique, mesurée ! 


par M. Pauthenet, est très différente : égale à 42,9 
et 54,1 CG. G.S. à 2940 K et 770 K respectivement 
pour la phase cubique, elle n’est que de 23,5 et 34,5 
pour la phase quadratique; la loi d'approche est 


très différente; linéaire en ï pour léchantillon 


cubique dès r 000 Oe, elle donne encore une forte 
courbure à 5 000 Oe pour la phase quadratique et 
traduit l'existence d’une très grande dureté magné- 
tique. 

C’est la mesure du champ coercitif que nous avons 
utilisée pour étudier les conditions de formation 
de la phase quadratique. Nous avons opéré avec 
des vitesses de refroidissement variées et, des prélè- 
vements trempés à partir de diverses températures 
nous ayant montré une variation continue du champ 


coercitif (et du rapport , des axes de la phase 


quadratique ) nous avons été amenés à rechercher 


- s’il existait un état d'équilibre fonction de la tempé- 
rature. 


Fig. 1. 
1. {oo C. Départ phase quadratique. 
2. 4oo0 C. » » cubique. 
3. 7000: Gt» »  quadratique. 


Les courbes 1 et 3 de la figure 1 montrent par 
exemple, qu’à des températures supérieures à 3000 C, 
on peut en des temps raisonnables, en partant d’un 
échantillon refroidi très lentement jusqu’à tempé- 
rature ordinaire, obtenir un état d'équilibre. Il n’en 
serait pas de même en partant d'échantillons 
cubiques (courbe 2) : le refroidissement lent permet 
sans doute l'apparition de germes qui favorisent 
l'établissement de l’équilibre, alors que dans les. 
échantillons trempés ces germes n’existent pas. 


3. La variation du champ coercitif. — Nous 
avons reporté sur la figure », en fonction de la 
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température à laquelle on a laissé l’équilibre s'établir, 
les valeurs du champ coercitif. Les points à tempé- 
rature ordinaire donnent les valeurs atteintes lors du 
refroidissement lent. Les différences entre les valeurs 
observées pour une même température, dans les 
deux séries, sont dues à une différence de granulo- 
métrie : le frittage des bâtonnets à 1 0000 C a réduit 
leur porosité et donc leur champ coercitif de Néel. 
Les irrégularités de la courbe 2 sont associées à 
une ségrégation de CuO, visible aux rayons X. 


300- 


100+ 
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Fig. 2. 
F (T) poudre. 
F (T) bâtonnets. 


ES 

2. H. 

4, La variation du rapport _. 
(courbe 1) 
mètre (£ — 1) en fonction de la température. Partant 
de x 


à 7600 C, . atteint 1,06 à l’ambiante. On constate, 


— La figure 3 


représente la variations du para- 


pour une température un peu supérieure 


. Lt à DENTS € 
d’une part, le fait très intéressant que e est une 


fonction de la température, et, d’autre part, que 
tout se passe comme si à température ordinaire la 
valeur d'équilibre n’était pas encore atteinte. Les 
énergies d'activation correspondant aux migrations 
d’atomes sont trop grandes pour qu’on puisse 


£ . C Fe $ 
observer une évolution vers un x Supérieur, à 
(4 


température plus basse. 


être considérée 


La quantité s —=7r peut 


commé ayant un rapport certain avec l’état d’ordre 
associé à l'établissement de la structure quadra- 
tique. L’allure générale de la courbe représentant $ 
en fonction de £ est bien celle qu’on a l'habitude 
de rencontrer lors de l'établissement de surstructures 
dans les alliages, mais il ne semble pas qu'on ait 
signalé jusqu'ici pour ces derniers une variation 


, Ce , € 
concomitante aussi importante [6] de ci 


Mes. 


\ 
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Nous donnerons dans la suite une interprétation 
possible de ce phénomène. 

La courbe 2 de la figure 3 donne les valeurs de a, 
en fonction de la température: la dilatation de la 
maille suivant c est donc accompagnée par une 
contraction suivant a. 


T: 


D 
A 
8,40 
F a FT 
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Fig=3: 


9. La valeur de À. —- Les intensités de réflexion 
dépendent de À. L'influence de À se fait évidemment 
d’autant plus sentir que les pouvoirs. diffusants 
des ions Fe*** et Cu*+ sont.plus différents. Avec le 
rayonnement de Co X « cette différence correspond 
à environ à électrons. Dans le calcul des intensités, 
le paramètre u qui fixe la position des atomes 


. Naturel- 


lement un écart & par rapport à cette Sa joue 
un rôle; nous avons choisi des plans pour lesquels 
l'influence de l'écart « est du second ordre sur les 
intensités, Pour le dosage de } nous avons comparé 
des raies dont les intensités varient l’une dans le 
même sens, l’autre en sens contraire de À. Enfin, 
pour nous placer dans des conditions optimum de 
photométrie, nous avons comparé des. raies d’inten- 
sité voisine. 


d'oxygène [11] a été supposé être égal à 


On trouvera en annexe un tableau de classifi- 
cation des facteurs de structure. 


En définitive, les couples de plans suivants ont été. 


comparés : pour la phase cubique (400) et (220), 
(40o) et (422); pour la phase quadratique (400) 4(004) 
et (220) + (202) d’une part, et (oo4) et (422), 
d'autre part. Le tableau ci-dessous donne les valeurs 
moyennes obtenues. 


Température 


d'équilibre (°C). Structure. A: 
CODE TÉL ES Cubique 0,34 ie 
SOON Cubique 0,33 
DODRTAr t Quadratique 0,35 

FAO EME NNE ST » 0,39 
DO Rene ». 0,43 


Ces valeurs ont été obtenues en corrigeant les 
facteurs atomiques d’après Hônl [7]. La correction 
de Bradley et Hope [8] aurait donné des valeurs 
différentes. Nous avons donc vérifié notre méthode 


sur le ferrite de Zn qu’on est d’accord pour considérer 


comme normal. Nous avons obtenu À — 0,06. En 
conclusion, si l’on ne peut accorder une confiance 
très grande aux valeurs mêmes données dans le 
tableau en raison de la mauvaise détermination 
actuelle des pouvoirs diffusants, du moins peut-on 
considérer que ce tableau donne un classement des 
valeurs de À et affirmer qu’au fur et à mesure que 
la température s’abaisse le nombre d’atomes de Cu 
sur les sites À diminue, autrement dit que le ferrite 
de cuivre tend vers l’état de ferrite inverse. 

L'évolution de À que nous constatons concorde 
bien avec la variation du moment magnétique : à 
l’état inverse, les ions ferriques des sites À et B, 
à coefficient d'interaction négatif, donneraient un 
moment résultant nul et seul resterait le moment 
du cuivre qui est faible. A l’état normal, au contraire, 
les ions ferriques additionneraient leurs moments 
et le faible moment du cuivre s’en retrancheraïit. 
On conçoit que le moment soit plus petit quand À 
est plus grand (cf. $ 2). 


6. Conclusion. 
quelle est la surstructure dont la formation donne 
lieu à la phase quadratique. Nous pensons que 
l'ordre s’établit sur le sous-reéeau B entre les ions 
ferriques et cuivriques. La figure 4 montre en projec- 
tion suivant l’axe c, dans l’hypothèse de À — 0,5, 
la structure que nous proposons sur le sous- 
réseau B, les plans de cuivre alternent avec des plans 
de fer; dans chacun des plans les atomes forment 
des files parallèles et, dans deux plans consécutifs, 
ces files sont perpendiculaires. On trouvera en 


annexe les données cristallographiques. 

On peut se demander pourquoi pareille structuré 
quadratique, qui est supposée exister à basse tempé- 
rature pour Fe,0, [13] n'apparaît pas pour d’autres 
ayant des 


ferrites caractéristiques magnétiques 


— Il nous reste à comprendre 


n 


du Es ee dr 


comparables. Un raisonnement par analogie peut 
en fournir une explication. Dans des chlorures tels 
que Co CL, Mg CL, etc., les cations sont entourés 
d’un octaèdre d’ions CI, tous équidistants du cation 
central, alors que dans Cu CI, [9], le rapport des 


_ distances des CI, voisins proches est de 1,28. Il y 


a tout lieu de penser qu’une structure ordonnée 
peut également exister pour 
mais que, les ions bivalents étant plus symétriques 
que ceux de Cu‘*, elle n’est pas quadratique. 


S Z 1: 3 
LE 83 DS 28 
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pates + OO © 
ue é © 
Fig. 4. — Projection de la maille suivant l’axe c. 


r . . fé 
Dans ce schéma, la variation de a? deux causes : 


on observe une variation du rapport des axes en 
fonction de létablissement de l’ordre et, simul- 
tanément en fonction d’une véritable variation de 
la concentration puisque les ions Cu viennent 
remplacer les ions Fe dans le sous-réseau B au fur 
et à mesure de l’abaissement de la température. 
Cette variation des paramètres en fonction de l’ordre 
et de la concentration est également classique dans 
les surstructures métalliques [10]. 


ANNEXE CRISTALLOGRAPHIQUE. 


Classification des facteurs de structure. — Les 
facteurs atomiques sont notés par leur symbole 
chimique. Les facteurs de structure en module sont 
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de la forme 
| F= Ars X ABtp, 
où 
NC AFeEEt Tr. 

Dans le Tableau I les différentes catégories de | F ! 
sont abrégées par les lettres a à f de la première 
colonne. La deuxième colonne contient les valeurs 
de la somme À? + X? +71? caractérisant les raies 
dont les | F| ont même expression. Les valeurs 
de Aisxn et Binxn sont consignées dans la troisième 
et quatrième colonne respectivement. 


TABLEAU I. 
Catégorie, R+ k+ EE, Au): Exp 
PR CR Drt) 10 +2Fe + Cu 0 
L Cu — 
b...: 16n Ær1 Fe + de (V2 1) 
2 
c.... I16(27 +1) 40 +1Fe— Cu f 
dE NAITIEEIS —410+2Fe : 2 
e,,.. 167 + 8 Cu —.2 
; £ Cu 5 
f 2e ROME Fe (V2 1) 


V2 


(RERO 120, SE) 


Cette classification permet d'interpréter qualita- 
tivement les différences d’intensité entre raies Debye- 
Scherrer de Fe,;0;, CuO cubique et de Fe,0;, ZnO 
(spinelle normal À = 0). - 

Les raies f, éteintes dans Fe,O, Zn O [A,sxn = 0] 
apparaissent nettement dans Fe,O;CuO Fr >o 
et Pire ol; raies (ur), (337), (531), (551), (gr). 

Les raies e sont plus fortes que leurs voisines de la 
catégorie c dans Fe,O;, ZnO; raies (220) > (400), 
(422) > (400), (620) > (444), (660) + (822) © (840). 
Dans Fe,0;, CuO par contre, les raies «e sont affaiblies 
[Burn <ol et les raies c renforcées [Biz > ol, 
de telle manière que les inégalités ci-dessus sont 
renversées. 


Dosage de À}. — En posant u — : + es(u est le 


paramètre de position des atomes d’oxygène [11]), 
l'incertitude sur A,,41 due à l’écart : est, aux faibles 
indices, du premier ordre en € pour les raies b et f 
et du deuxième ordre, donc négligeable pour les 
raies €, d, e, utilisées pour doser . Les intensités sont 
définies par 

Tn,k,5=p|F}, 


où p est la multiplicité du plan (Akl). L'intérêt de 
la méthode réside en ce que le rapport d’intensités 
de raies telles que c et e ou d et e varie largement 
avec À. Dans l’exemple numérique suivant, il varie 
de 0,7 à 2,45 pour 0 € À < 0,5 [raies (400) et (422) 
de Fe,0; CuO cubique] : 

L(s001 __ 6(27,1+19,4X}? 


Tusay — 24(16,2—8,8À} 


Les intensités observées sont corrigées de la manière 
usuelle des facteurs de Lorentz, de polarisation et 
d'absorption (Claassen). 
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Le ferrite de cuivre quadratique. -— Le Tableau II pre Eee AR SG ARR CS 
résume les valeurs de sin?@ observées et calculées Indice Sir © si 0 He NOR : 
(K « de cobalt) et la catégorie des raies de Fe,0;, CuO É AU: | CODE NORGE TES RON ne PAC AUERS 
quadratique, décrit ici pour des raisons de commodité (351, 711).:.. 0;6020 0,6016 f f 
comme un spinelle pseudo-cubique. Les atomes d’oxy- CAD 0,6148 0,612 Î se 
gène dépendent alors de deux paramètres u et », (RATER 10 ,6336 7050898 me 
in dé 3. (ORNE 0,6570 0:6668/7RE e 
8 2 (535)....... -0,6660 _ - 0,6649 m b 
= Test LE c£ (DOSYALES  0,6740 0,6734 Î a 
F'errite de cuivre quadratique… (55 713).. 0 ,6863 0,6857 m+ b. 
a=8,216 À; = 1,060. | (TRS 076964 0 SOU ES SERA 
: a ne Ÿ (SOIF 0,7970 0,7965 m+* a 
acteur — 1Q4 
UE 250 rt ae (157)...:.., 6,8230 0,8229 m b 
(AK). (obs.). (cale.). Intensité. structure, (660)... E = ..0,8510 0,8511 ff à e 
CASSER 0,0321 0,0341 mn f : (555, 115)... o,8547 0,841 F b 
COOPER 0,0879 0 ,0894 m ne (AUS) 0,8627 0,8625 m , c 
(90) SRE: 0,0939 0 ,0046 f e (ÉUR TE © 0,8849  o,8852 m to 
(413) OT RTR AVES O,1183 . F ae (SOA) 0,9245 0,9249 m e é 
E (ALES 0,1270 0,1287 FE: b (840.2. 21050495 : 0,9456 * mn £ k 
COPIER" 0,1393 0,13067 2 FE d She A 
(ODA ALERT 0,1670 0,1683 mm ( Notations des intensités :. di faible très 
(AOD) FETE Es 0, 1880 0,1891 m+ ti faible, m = moyen, F = fort, FF = très fort. Se 
COS) ER 0,2120 0,2129 ff fers < L 
(331)....... 0,2225 0,2233 ff f La due ordonnée, ne. ici (fig. 4), dérive 2 
COIN RUE 0,2626 0,2629 f -e d , | LINE 
; Here ue u groupe d’es Le NE 
(GONE EX 0,2774 0,2785 mn e= Ê - pate Dre Le [2] (ou : ddm 
(44B)...1.:. 0,2897 0,2867 Rs Ve par Île dédoublement du complexe de réseau en 8 d. 
(338). 2. 03087 0:3075 à k Le vrai groupe d'espace serait C2° — 7 mm. La dite 
É ee 03183 o STD F b petite maille est celle formée Re c et les moitiés des 
(404)... 0,380 0,359 FF a diagonales de la figure 4 « = © Vo. 
Sas re Sol 0,3783 = F er Les raies de surstructure Se le. Héaoes ë 
(931). ...... 0,4130 0,4124 ff Î  blement des sites B donnerait lieu n’ont pas encore … 
COOP re 0,4265 0,4261 ff Ge été observées, soit à cause de leur très faible intensité, 
(620) 727278 0,4719 0,4728 ff è soit parce que l’ordre à grande distance n'est pas 
(BIS NES, 0,4742 0,47b8 FE b ÉRRRER 
(833) 0: 2976 0 :4966 ee > Lex de plusieurs états ordonnés au- HA É 
(892). He. D i50 SO Es d dÉRE 05 ne peut être exclue a priori. ee LS ce 
(MIN ee A ir ce RE 
(HAS VER 0,9711 0,5704 Î f Manuscrit reçu le 3 janvier 1950. 
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CONTRIBUTION À LA THÉORIE DES CHAMPS NON LOCALISABLES 


Par BERNARD KWAÏ.. 
Institut Henri-Poincaré, Paris. 


Sommaire. — Les champs non localisables © dépendent à la fois des opérateurs quantiques æ, et p". 
Les équations relativistes des champs scalaires, bâties par -Yukawa conformément au principe de 
réciprocité de Born, introduisent une nouvelle constante fondamentale à, jouant le rôle d’un rayon fini 
| des corpuscules élémentaires. En posant m = 0, mais en laissant À 0, on obtient une Mécanique 

ondulatoire réciproque de celle ordinaire des corpuscules ponctuels; le corpuscule réciproque étant 
caractérisé par une masse nulle et par une extension spatio-temporelle finie. Les équations de Yukawa 
admettent comme solutions des produits de certaines solutions de la Mécanique ondulatoire ordinaire 
et de celles de la Mécanique ondulatoire réciproque. 

A l’approximation de l Optique géométrique, les équations de Yukawa conduisent à une sorte de 
Mécanique classique, caractérisée par une fonction de Jacobi, qui dépend de deux groupes de variables X, 
etr, par Pintermédiaire de l'expression S, (X.) f (r,), où f (r,) est, dans le référentiel propre une tonction 
arbitraire des points de la surface d’une sphère de rayon À, tandis que S, (X) est une fonction de Jacobi 


ordinaire. ; 
Dans l’espace de la nouvelle variable Ty S (X4, Fr.) est donc partout nulle sauf à la surface d’une 


sphère de rayon À. 

En ce qui concerne la théorie des champs spinoriels, il semble possible a priori de linéariser les équa- 
tions scalaires et d'introduire le spin indépendamment dans les deux espaces, celui de configuration et 
celui des moments. Les équations du type Dirae, Valables à la fois dans les deux espaces, sont possibles 
et compatibles entre elles, lorsque æ* est un quadrivecteur du genre « temps », ee solutions corres- 


pondantes n'étant nulles qu’à l'extrémité d’unintervalle de temps propre égal à 2. Pour obtenir les 


équations relatives à un quadrivecteur æ* du genre « espace », on fait appel à une règle, énoncée par 
l'auteur en 1936 et qui permet de déduire les équations d’onde, relatives à une grandeur tensorielle 
attachée à un corpuscule élémentaire, de la transformation de Lorentz, écrite dans la représentation 
cayléenne, dans le passage au référentiel propre du corpuscule en question. Un quadrivecteur de genre 
« espace », se réduisant dans ce dernier référentiel au vecteur spatial de composantes r4, T», Ty 
avec r? + r; + r; = }?, on trouve que les équations d'onde correspondantes s’écrivent 


e (CŒREErEE Fr) S 


= 0; 


TV, (ue = 1, 2,3, 4)etT% (x, 2, 3) étant deux systèmes commutant des matrices de rang 8, satisfaisont 
chacun aux mêmes relations qui régissent le système des matrices de Dirac, et S = (S;, Ss, Sa, Sa; 
Si,83, Si Si). Grâce à cette représentation, il est possible de développer une théorie des champs spi- 
_noriels non localisables S(X*, p*) qui dépendent d'un quadrivecteur du genre « espace » æ* el d’un 


quadrivecteur du genre « temps » p*. 


L'objet de la théorie des champs non localisables 
est d'étudier les champs d’onde plus généraux que 
ceux dont s'occupe la Mécanique ondulatoire ordi- 
naire, en les considérant précisément comme des 
grandeurs qui dépendent à la fois des +’ et des p”. 
Ces fonctions d’onde de l’extension en phase doivent 
obéir à certaines équations dans l’espace-temps 
qui se réduisent aux équations d’onde ordinaires 
de la Mécanique ondulatoire, lorsqu'on les suppose 
ne pas dépendre des p'. Elles doivent, d’autre part, 
obéir à certaines autres équations dans l’extension 


1. Introduction. — En Mécanique quantique 
les coordonnées de l’espace-temps zx, et celles de 
l'extension en moments ne sont pas mesurables 
simultanément avec une précision infinie, car elles 
sont représentées par des opérateurs non commu- 
tables satisfaisant aux relations de non-commu- 
tation 

a pu — pad =[A", pu] = ho. (44) 
Toute grandeur physique, fonction des 2’ seuls, est 
localisable dans l’espace-temps, car l’on peut faire 


une mesure simultanée portant sur cette grandeur, 
et sur l’endroit où s'effectue l'observation, mesure 
qui n’est affectée d'aucune limitation théorique 
quant à sa précision. 

Il n’en est plus de même lorsqu'une grandeur 
physique est fonction 
car elle se présente alors comme un opérateur qui 
ne commute plus avec æ’ (ni avec p,) et sa mesure 
ne peut plus se faire avec une précision infinie 
en même temps que celle de la coordonnée +”, qui 
la « localise » dans l’espace-temps. 


la fois des +’ et des pu. 


en moments, et la manière la plus simple d'établir 
ces équations consiste à transposer les équations 
de l’espace-temps, en prenant comme guide Ja 
symétrie entre les x’ et les p, que révèlent les 
équations de la Mécanique analytique classique- et 
qui se retrouve dans la théorie des transformations 
de la Mécanique quantique. C’est à cette propriété 
de symétrie que Born [1] donne le nom du prin- 
cipe de réciprocité, et c’est à ce principe que fait 
appel Yukawa [2] pour bâtir les équations de la 
théorie des champs non localisables. 


ee 
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2. Le principe de réciprocité de Born. — 
M. Born part de la constatation que les équations 
de la Mécanique classique sous la forme canonique 


dx, 0H dp, ol 


= = é (2.1) 


dé 0py’ dé DES 


sont invariantes par rapport aux transformations 
du groupe de Lorentz mais également par rapport 
aux transformations plus générales de contact, et, 
elles restent inchangées lorsqu'on effectue la subs- 
titution 


Ly> Py; PT — y: 


Cette sorte de symétrie apparaît également, en 
Mécanique ondulatoire, où le mouvement d’une 
particule libre est représenté par une expression 


exp (PDA symétrique en x, et p,. La même 


symétrie apparaît dans les lois fondamentales de la 


Mécanique quantique, dans les lois de commutation 
(1.1), dans l’expression de l’opérateur-tenseur du 
moment de rotation 


(2.3) 


ainsi que dans la théorie générale des transformations 
dont le cas particulier est représenté par la trans- 
formation de Fourier, permettant de passer de 
l’espace des x, à l’espace des ps. 


My Lypy— Pal; 


n i y 

Dp)= |. WCr)er "ar. (2.4) 

Ce sont ces considérations de symétrie qui ont 
suggéré à Born le principe de réciprocité, selon 
lequel les lois de la nature seraient symétriques 
par rapport aux espaces de configuration et des 
moments, 
invariantes par rapport à la substitution 


LD; Dire by: 


Remarquons tout de suite que le principe, ainsi 
formulé, est très général et peut se prêter à des 
applications aussi variées que différentes. Il suffit, 
pour s’en rendre compte, d’examiner l’évolution 
des différentes tentatives théoriques, qui ont été 
faites. par Born lui-même, depuis sa théorie de 


l’univers clos dans l’espace des p qui serait « réci-, 


proque » à l'échelle microscopique de l'univers 
sphérique d’Einstein dans l’espace des x à l’échelle 
cosmique, jusqu’à sa théorie récente des masses 
des corpuscules élémentaires. Une autre idée de 
Born est que la constante de Planck À doit être 
considérée comme le produit d’une longueur mini- 
mum & par la grandeur réciproque, le moment 
maximum b (selon Born, la grandeur «a serait voi- 
sine du rayon classique de l’électron). 


3. Équations des champs non localisables 
de Yukawa (Champs scalaires). 
localisables ordinaires W(x) doivent être considérés 


(2.9). 


ou, plus précisément qu’elles seraient 


comme cas particuliers des champs non localisables 
Q(x, p), et sont représentés par des matrices dia- 
gonales dans la représentation où les x le sont. 
Dans ce cas, 
d'ondes de Gordon-Klein. 

Inc 


(Ou 2) W(x)=0, . x = Fe (31) 


qui peut être écrite encore sous la forme suivante : 


[Pa PAR w]] nc —0, 
[A4, B] = AB — BA. 
Yukawa admet que le champ scalaire non locali- 
sable satisfdit une relation de commutation de la 
forme précédente et comme on a besoin en outre 
d’une relation de commutation entre Q(x, p) et 
les x, il se guide par le principe de réciprocité et 
aboutit ainsi au système d'équations suivant : 


| Pus LP b, 0]] + CAO 
[es [æb0 |] —2= 0. 


Remarquons tout de suite que si l’on devait appli- 
quer strictement le principe de réciprocité, on eût 


(3:3) 


h 
posé À = ne? mais ce n’est pas le cas ici, car À doit 


figurer, comme nous allons le voir tout à l’heure, 
le rayon du corpuscule et ne peut pas être identifié 


avec la longueur d’onde de Broglie 


On voit immédiatement que, rs Q ne dépend 


pas des p et commute donc avec les x, on doit 
avoir À —0o, de sorte que l'intervention de cette 
longueur attachée au corpuscule élémentaire est 
intimement liée au caractère non localisable du 
champ ©. 

Dans la représentation où les opérateurs x, sont 
des matrices diagonales, chaque élément de la 
matrice (x, |{2]\x,) peut être représenté par une 


fonction de deux variables X, = (a + ay) et 


Ty = (ty — ty), et, dans ce cas, le système prend la 
forme que voici 


d? mn? c? > 
EE EP EE EE Q) ) — 
ce OA TE ) (AY, r#)= 0, 


(3.4) 
(rarb— \)Q0(4u, 7B)= 0, 
système qui admet les solutions suivantes : 
OA, TE) = o(4, Ky)d(rurb— 22), 
CA hs À [UE Ku.) ( u ) | (3.5) 
K,,KW+ m°c*= 0, 


où à est la fonction singulière de Dirac. On voit 
donc que les fonctions d’onde € dépendent de deux 
variables, l’une X, décrivant pour ainsi dire le mou- 
vement d’ensemble (du centre de gravité du cor- 
puscule), tandis que l'autre r, correspond à un 
degré de liberté supplémentaire, interne, de la 
particule. 

Mais le problème: qui nous intéresse est de pouvoir 


N°5. 


ils doivent satisfaire à l'équation 


alba tan dc Es ANSE 


is COR EAN D‘ El < Me RES 


identifier la longueur À avec le rayon propre de 
la particule, il faut donc que dans le référentiel 
propre, caractérisé par X, = K, — K3 = 0, K, =imc, 
les solutions soient nulles partout sauf sur la sphère 
de rayon À. Pour qu'il en soit ainsi, force nous est 
d'ajouter aux équations (3.3) l'équation supplé- 
mentaire 


[zy, [pt, Q]| = [ pv, [ty 9 Le 10 (3.6) 
qui, dans la représentation HE en employée, 
s'écrit 

- dO( AS, 74) ° 


1 (3.7) 
MIRE 


} et conduit aux solutions du système total qui ont 
Ja forme que voici : 


OA, rb) = (Au, Ku)0(rurt— 2) 0(ru KV), (3.8) 


dont le dernier facteur d(r, K;) contribue, dans le 
référentiel propre de la particule (K, — K, = K3= 0, 
K;= imc) à annuler r;, de sorte qu’il ne reste que 
les composantes Ti lo €t T3, AVEC Ti tr +ri = 
Ainsi, est-on arrivé à bâtir un système d'équations 
dont les solutions jouissent des propriétés que nous 
leur avons voulu faire attribuer. 


4. Mécanique ondulatoire réciproque [3]. — 
La Mécanique ondulatoire ordinaire d’une par- 
ticule ponctuelle est contenue comme cas limite 
dans le schéma précédent. Les solutions £(x) de 
la Mécanique ondulatoire ordinaire satisfont en 
effet aux équations du schéma précédent où l’on 
_ a posé À — 0. Il existe un autre cas limite en consi- 
- dérant des solutions I[(p), qui correspondent 
à moc — 0 et À; <o. Ces champs, qui ne dépendent 
pas des coordonnées de l’espace de configuration, 
mais uniquement de celui des moments, décrivent 
une particule réciproque (de celle qu’envisage la 
Mécanique ondulatoire ponctuelle), caractérisée par 
l'absence d’une masse propre et par une extension 
spatio-temporelle finie À. L’interaction des par- 
ticules réciproques ne peut se faire que par l'in- 
termédiaire des fonctions d’interaction U qui ne 
dépendent que des moments p, et, dans l’espace 
des p l'équation d'ondes s'écrit 


{ Cl LS 2 32 l of = 
lapepr al [U(p) + »] {2(P) 0. 
Le mouvement de la particule réciproque 
libre sera donc représenté par des ondes du 


type exp | (3) p| avecr,r’— }?= 0,se propageant 
dans l'extension en moments. Et de même qu’en 
Mécanique Ondulatoire de M. Louis de Broglie, 
au moment P, est associé l'inverse d’une longueur 
d'onde, caractérisant un champ périodique dans 
l’espace de cénfiguration, on a ici,’par réciprocité, 
associé au vecteur r, représentant le « rayon fini » 
de la particule, linverse d’un moment p,, cet 


7 
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| ne A MES, 
IAA 


u 
J 
3 L I , A 
inverse ss = E Tv jouant lé rôle du vecteur 
Ë à 
onde dans l'extension en moment de la même 


5e I I à 
manière que K,— rP le joue dans l’espace 


SET 
de configuration. 5 

En particulier, si nous enfermons la particule 
réciproque dans une «boîte» dans l’espace de p, 
c'est-à-dire si la valeur absolue de son moment ne 
peut varier qu'entre deux limites dont la différence 
est égale à P, on obtient une quantification du 
vecteur interne r, de la particule, grâce à la condi- 


tion cyclique 


n, étant des entiers, satisfaisant la relation 


NNNV = :X2 ë . 
4 h 


Soient alors £(x) et r(p) des solutions de ‘la 
Mécanique ondulatoire ordinaire et de la Mécanique 


ondulatoire réciproque, alors les fonctions 
Qi(x, p}=E(x)r(p) et (x, p)=r(p)E(x) 


satisfont aux équations (1), à condition que l’on ait 


E(æ}r(p)=:0 ou r(p)E,(æ)= 0, 
où :, est le résultat de l’opération p, effectuée 
sur £(x), et 7'’(p), le résultat de l'opération X” 


effectuée sur 7(p) 


E(z)=EnE(e), 2 (p)= fe xp) 

5. Approximation de l'optique géométrique 
des équations de Yukawa. — On sait que les 
équations de la Mécanique ondulatoire ordinaire, 
en approximation d'ordre zéro, lorsque le quantum 
d'action À est petit devant les grandeurs d’action 
mises en jeu, fournissent des équations de la Méca- 
nique analytique classique sous la forme de Jacobi- 
Hamilton. Nous allons étudier les équations qui 
résultent d’une approximation analogue dans le 
cas de la théorie de Yukawa. | 

A cette fin, nous allons nous servir de la repré- 
sentation (X,, r,) dans laquelle on a 

m° c° 


Je CE ri) = 0 | 


(ryrk=—%2)0= 0, 
En posant alors 


O( Au; l'y) © exp [ SA, Fr) Si, (5.2) 


on est conduit, en approximation d'ordre zéro. 
par rapport à h, aux équations que voici : 


[Sen fe ses } 2 0, | 


l OX, 
Rte D) 


PyI*— = 0 r'y - 
V : 04, 


= 
© 
mm. 
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donc la fonction $S dépend de deux tu >. 
et r, et peut être mise sous la forme S$(X;) f(r,). 
La fonction f(r,) n’est définie, dans le référentiel 
propre du corpuscule (dans l’espace X,) que sur 
une sphère de rayon À et, peut donc être représentée 
par un développement -en série des fonctions sphé- 
riques. En considérant que /f(r,) est normalisée 
sur la surface de la sphère et en intégrant sur cette 
surface la première des équations (5.3), on trouve 
que S,(X,) peut être identifiée avec la fonction de 
Jacobi-Hamilton de la Mécanique classique. 


6. Les champs non localisables du corpu- 
SUR ASE ? D'or 
scule à spin [4].— Nous généralisons avec Yukawa 


des équations de Dirac de la manière que voici : 


[pus YHQ] +imcO—0; VV YYYR = 20 (GE 1) 
d'où, par itération, on obtient bien 
[pv [Puy o]| — m20c2Q = 0: (6.2) 


Dans la représentation où les x sont diagonaux, 
on peut considérer l'équation entre opérateurs t#, p 
et Q(xs, p) comme agissant sur une certaine 

: : ee 
fonction f (Ru). Lorsque l'opérateur (x, pu) 
dépend pas de p,, et se réduit à W(xt), on a 


fépa(Wp) —EWpu(f)= YECpe WT) 
et l'équation (6.1) devient 
J(æ)(py 


d’où résulte l’équation de Dirac ordinaire, relative 
à la fonction d'onde W(xr). 

Il est plus délicat d'établir d’une manière correcte 
les équations réciproques. Tout d’abord, est-ce 
que la réciprocité doit nécessairement s'appliquer 
aux spins ? Nous avons vu que dans la représen- 
tation où les 4 sont diagonaux, tout se passe 
comme si le champ dépendait de deux variables, 
l’une x décrivant le mouvement d'ensemble de la 
particule, et l’autre r" représentant des degrés de 
liberté supplémentaires, internes de la particule. 
On peut donc très bien concevoir à priori que la 
particule possèdé un spin dans l’espace de X# mais 
non dans celui des r#, Nous avons d’ailleurs constaté 
que la théorie du champ scalaire, admet des solu- 
tions de la forme (x, p) —(x) t(p),. où (x) 
et z(p) sont certaines solutions de la Mécanique 
ondulatoire ordinaire et de la Mécanique ondulatoire 
réciproque, de sorte qu'on doit pouvoir effectuer 
la. linéarisation de l'équation d'onde, introduction 
du spin, indépendamment dans l’espace de æ et 
dans celui des p. Et, en effet, on vérifie immédia- 
tement que, par exemple, les équations suivantes 
sont parfaitement compatibles entre elles 


YEW + mc) = 0, 


[Ps YHQ] +imcQ 


[zvLpu, 01] = 


(l 
ds 
3 
E 
ra 
Æ 
|æ 
as” 
| 
© 


Considérons maintenant les équations 


Le, 21 + EX0 =, | (6.4) 


Montrons, qu’en effet, on peut les rendre compa- 
tibles avec les équations (6.1). À cette fin appli- 
quons à l'équation (6.1) l'opération [x”’, y,,...] et 
à l'équation (6.4), l'opération [p,, y+, 

Te [2 Lu 01]= 
Mais puisque 


Fev, LPu: Ed >]}= 


on doit donc avoir 


EN — véy,[ pu FæŸ: 0]|: 


pu Le, ef] Go 


Gt tas) [ru Le 07] = 0, 


d’où, on tire, en vertu des propriétés des matrices y, 
la condition de compatibilité suivante : 


bre EX 0]] — Lr*, [æu: o | 
=|4,, [pb 0]| — | æv, [P; o]| 10: 


qui, dans la REA RUES (ZX r), prend la forme 


(6.5) 


° que voici 


|: AUX, Fr) :  dU(X r) 
(EE ÉNSEE 


Considérons alors le référentiel propre du corpu- 
scule, dans lequel son mouvement d'ensemble est 
représenté par une onde de vecteur-onde 


K; = Ko = K3 = O0, K, = im C. 


Les équations montrent alors que l'en a 
=== 0etr, 0, de sorte que la solution 
des équations s’écrive sous la forme 


u(0, 0, 0, Emoc; 0, 6, 0, 7»)0(— r$ + A)errcX, (6.7) 
et ne diffère de zéro qu’à l'extrémité d’un intervalle 
: ; X 
de temps propre égal à + ne 
On peut se demander si de telles solutions ne 
puissent convenir à des corpuscules essentiellement 
éphémères, pour lesquels lx donnée de la durée 
de vie jouerait un rôle aussi fondamental que celle 
d’un rayon fini dans le cas des corpuscules essen= 
tiellement stables, 
Il reste donc à trouver de quelle manière doit-on 
écrire les équations - d'onde relatives à l'opérateur æ, 
lorsque celui-ci est un quadrivecteur du genre- 


(1) Il suffit d'envisager l'identité ! 
[4, 2, ci]+[2,fc, 41]+[c,14,281]=0 


el la relation 
Fa p,)=in8,. : 


linéa- 
risées en æt. Lorsque celui-ci est un quadrivecteur 
du genre « temps » comme l'est p,, il n'y a aucune 
difficulté pour écrire les équations d’onde corres- 
pondantes et nous pouvons poser 


1], 1sviénts 


=0: (6.6). 


et, qui permet de trouver, 


espace. Pour résoudré ce problème, nous allons 
faire appel à une règle, que nous avons établie 
en 1936 [5] à l’occasion de l'équation de Dirae, 
d’une manière seule 
correcte croyons-nous, les équations d’onde rela- 
tives à une grandeur tensorielle quelconque, atta- 


x 


chéè à un corpuscule élémentaire. 


7. Règle fondamentale de formation d’équa- 
tions d'ondes relatives à une grandeur tenso- 
rielle, attachée à un corpuscule élémentaire. — 
Nous allons formuler cette règle de la manière que 
voici 


Les équations d'ondes relatives à une grandeur 
tensorielle, attachée à un -corpuscule élémentaire, se 
déduisent de la représentation cayléenne de la trans- 
formation de Lorentz que subit le tenseur en question 


lors du passage du référentiel de l'observateur au 


référentiel propre du corpuscule, en remplaçant dans 
la formule ainsi obtenue les composantes du tenseur 


dans le premier référentiel par les opérateurs quan- 


tiques qui leur correspondent. 


Ainsi, en considérant le quadrivecteur quantité 


RES W) 
de mouvement énergie P(P, PS PRE _ }» la 


formule de Cayley s’écrit-elle sous forme du produit 
quaternionien suivant 
PS= SP, (TA) 


S étant le quaternion de transformation. Dans le 


référentiel propre P'— (0, 0, 0, im,c), etla formule (7.1) 


se simplifie 


PS = iMmocS*, (T9) 
cette équation jointe à sa conjuguée —— inopcS 
est équivalente à l'équation de Dirac et à se conju- 
guée 

YYP,Ÿ = imocW, PEU mo ce (18)2 
comme on le constate en posant 
L'on = Si + S4 ne AE US S > } 
Wa=— Si Si +i(— S>+S5), 
LR — S3— SJ m5 t(So— 19% }. 


En appliquant notre règle à un vecteur du genre 
espace 2%, qui dans le référentiel propre du corpu- 
cule est représenté par un vecteur spatial de compo- 
santes 7, Ta» ls AVEC Ti +13 + T3 
de Cayley nous conduit à la relation que voici 


AS=S"R,  R=(r, re, rs 0). (7.5) 

Mais cette fois-ci, il n’est plus possible comme 
dans le cas de l'équation de Dirac, de trouver une 
représentation complexe à quatre lignes et à quatre 
colonnes, qui permet la séparation des équations 
complexes conjuguées. Force nous est donc de 


— }?, la formule” 
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recourir à une représentation à huit lignes et à huit 
colonnes, englobant les fonctions complexes conju- 


guées, ou, Ce qui revient au même, envisager une 


représentation réelle de rang 8. 

Nous allons montrer comment on peut trouver 
rapidement une représentation de rang 8 dé la for- 
mule de Cayley, en se basant sur nos anciennes 
recherches sur la représentation des quaternions 
au moyen des matrices de rang 4. Nous avons 
montré [6[ qu’il en existe deux, qui commutent 
entre elles, chacune introduisant deux unités quater- 
nioniennes indépendantes, de sorte que leurs produits 
de deux à’ deux définissent un système hyper- 
complexe à quatre unités de base (comme l’est le 
système des matrices de Dirac). Avec nos notations 


antérieures la formule de Cayley s'écrit 
(Ps. .1+emPE)S =(P;,.1—exP#)S"*, (7.6) 


e'k et ek étant les deux systèmes d'unités quater- 


_nioniennes indépendantes et telles que 


Cp —= €? = 1 
x — — n 


QUPEENTE 


Fig 


e162 = €3, RRQ eRe}— Re;er RG) 


(J £ = 1,2, 3). 


Le système de 16 matrices £, =e,e.(v,u=1,2,3,4) 
permet la construction. des formes bilinéaires, qui 
définissent différentes grandeurs tensorielles, selon 
le caractère relativiste des fonctions qui servent à 
les former. Ainsi, en partant des composantes 
d'un vecteur ou d’un tenseur antisymétrique gauche 
forme-t-on des tenseurs du second rang, le dernier 
étant celui, bien connu de Maxwell. 

Dans le cas des demi-vecteurs $S, la théorie 
conduit à envisager deux êtres géométriques, selon 
que l’on forme des formes bilinéaires selon le schéma 


SLR (7.8) 


ou selon le schéma 


SRE. (7.9) 


Dans le premier cas, on forme un quadrupode 
vectoriel dont les composantes f,, et £,, forment 
le quadrivecteur courant et le pseudo-quadrivecteur 
de spin, bien connus en théorie de Dirac. 

Dans le second cas, nous formons un tripode 
bivectoriel dont un des membres est le tenseur 
antisymétrique gauche, représentant le « moment 
magnétique et électrique >» de l’électron de Dirac. 

Nous aurons besoin de tenir compte de l'existence 
de ces quadrupode et tripode, lorsqu'il s’agira 
d'écrire les équations d’onde relative à une grandeur 
quadrivectorielle du genre espace, attachée à un 
corpuscule de spin 1. 

En revenant à la formule de Cayley et à sa conju- 
guée nous voyons qu'elles peuvent s’écrire | 


(P,.1+ ex PE) S = (PI :1— ex PA)S* (r.= 1 ie) 


| 
(Pr ex PH) S"=(Phi+er POS (Pi=—P:), | 


4 


équations qui peuvent s’écrire à l’aide des matrices 
à huit lignes et à huit colonnes, de la manière que 


voici 
f < : e}, 0 l S 
| P,;,.1- = ax er Pr si e 
70% 0. —ex |, S ; 
à. Bacs | + ei (&): A1) 


(L'application de l'opération 


! 
ek O0 


conduit bien à la relation 


P? + Pr; P# = Pi ES P' a 


Dans le cas d’un vecteur 


genre « temps », nous 
avons donc : 


1e e’. 0 à l D 
| Pi Le Fr e Px ï S* 
= 1m)C Ses | : . (7.12) 


Et dans le cas d'un vecteur genre « espace » :. 


la Fe, e'k 0 LES 
due Ft pee ee 
l One ; Si 
OI Er S : ; 
VE He 
x à 2 | s*P (7.43) 
avec 
TI HIS +ri = D». 


Nous allons condenser l'écriture de ces équations 
en écrivant S à la place de la colonne S,, S,, Szx Sa: 


Si, S:, S;, S; et en introduisant les nat de 
2 PA ELES — 
rang 8, suivantes [es [pb 1] Rs | FR 
HET ES \ prenant, dans la représentation (X, r), la forme 
rs = | Fe l : SE 7) TR 
Ty = e— = 0. (8.5) 
e} 0 O0 —e ë JA da 
ENS Re ms ; (7.14) Les EE ; 
0 —e% €k. 0 Et c’est ainsi que nous parvenons bien aux solu- 
M OR nd , tions, qui, dans l’espace de r sont nulles partout, 
< € 0 sauf sur la sphère de rayon À. 
Manuscrit rech le 18 janvier 1950. 
AMOGRARHIE, 


[1] Born M. — Proc. Roy. Soc., 1938, 165, 291; 
Born M. et GREEN H. S. — Proc. Roy. Soc. Edinburgh, 
1949, 62, 470. 
[21 YukawaA H. — Phys. Rev., 1949, 76, 300 et 1781. 
[3] KwAz B. — C. R. Acad: Sc., 1950, 230, 184. ? 
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qui vérifient les RS nantes CE 
LyTy+ Tulv= 20g 
él + 17 le = 20/4, 
NA DES De 
2,3,4; K, 41,2; 3). 
Cela étant, nous pouvons écrire les équations 
précédentes de la manière que voici : 


F0 


Go 


(pince) S =0, l 


Fe 7.46 
(Déry— rl) = 0, KF:19) 


TRIEZ. 

8. Champs spinoriels non localisables dans le 
cas d'un vecteur réciproque 1, dugenre «espace». 
— Il nous est facile maintenant d'écrire les équations 
auxquelles doit satisfaire un champ È(t,, p') non 
localisable, généralisant le champ demi-vectoriel S(x,). 
Nous poserons 


[pY TS, E]=EmocE = 0, 


(8.1 
[ap TEE] — 0, (En 


rRrt= }2: ( 


Pour obtenir les conditions de compatibilité de 


ces équations, nous allons appliquer à la première 
d’entre d'elles l'opération [x,, Lv, 
seconde [p", P,, ..: 


obtient ainsi 
CET, [au Ep*,.5 1| = imcril'E =— ryru] p*, [zu2]}: 
D'où on tire, compte tenu du fait que 
[Su Cr 21] =?" [ew 51} 
quels que soient : et », et des relations 
PET, + D,Tu= 2 ôf, 


læicondition de compatibilité suivante : 


[4] Kwaz B. — C. R. Acad. Sc. dose 230, 276. 

[5] Kwaz B. — C. R. Acad. Sc., 1936, 203, 983; J. Phy- 
sique, 1937, 8, 81. 

[6] C. R. Acad. Sc, 1934, 198, 1582; J. Physique, 1934, 5, 
445; Bull. Sc. malh., 1935, 59, 328. 


1 


7] cet: à; 108 
J, les matrices r pouvant d'aik 
leurs être mises en dehors de la parenthèse. On 


| 
| 


SUR LES ÉQUATIONS INTÉGRALES DE L'HYDRODYNAMIQUE 


Par ANTONIO GIAO. 


Sommaire. — Les équations intégrales de l'hydrodynamique, donnant en chaque point du domaine 
occupé par un fluide les valeurs des inconnues (vitesse, pression, densité, etc.) en fonction des conditions 
initiales dans tout le fluide et des valeurs aux limites sur les surfaces frontières, ont été déduites jusqu’à 
présent sous certaines conditions restrictives (mouvements lents où à tourbillons négligeables, fluides 
incompressibles, viscosité constante) faisant disparaître les termes non linéaires des équations aux 
dérivées partielles. Nous exposons dans ce Mémoire une nouvelle méthode qui généralise la méthode 
d'Oseen et permet d'écrire les équations intégrales de l'hydrodynamique pour des mouvements 
quelconques et quelle que soit la compressibilité ou la viscosité du fluide. 

Nous pensons que les nouvelles équations sont susceptibles de nombreuses applications, mais nous 
nous bornerons ici à une application à la théorie des perturbations, afin de montrer comment la 
connaissance des actions extérietires qui agissent sur un fluide permet de résoudre le problème de la 
prévision de l’évolution du fluide par la déteremation de la perturbation en fonction des seules 
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conditions initiales. 


I. — Le problème général de l’hydrodynamique. 


1. Les équations fondamentales. — Quand on 
étudie le mouvement interne d’un milieu fluide 
au point de vue eulérien (variations des propriétés 
en un point fixe du référentiel choisi) les équations 


_ fondamentales aux dérivées partielles de l’hydro- 


dynamique sont les suivantes en notation vectorielle 


FT ANRT PR RES ES 


ca Pie : 
d. mo VERNRTV = a NP ae 
dz 


S (1) 
dt 


œ 


_ Are CE 


Voici la signification des différents ee É 


Ÿ, vitesse des particules fluides par rapport à un 
référentiel R (x, y, 2); 
pression; 


Te nur L 
volume spécifique (densité po = :); 
forces extérieures (en général la pesanteur); 


vitesse angulaire de rotation du référentiel R par 
rapport à un référentiel galiléen R,; 

constante arbitraire (que nous prendrons par 
exemple égale à la valeur moyenne de 2 dans 
tout le domaine occupé par le fluide, & étant 
le POPICQUE ie viscosité) ; 


4 + 
F=F—: AY, où F est la force de frottement par 


unité de masse; 
d 
£, temps, ? ER ; °pérateur de dérivation partielle par 


LA: au temps; 


V, opérateur  hamiltonien ou opérateur  gra- 
dient:(V.= div, V-K= rot; VV A,; apla- 
cien); 

X et ., multiplications vectorielle et scalaire. 


A l’exception du cas particulier d’un fluide incom- 
pressible à viscosité constante et pour lequel on 
admet qu'il est possible d'exprimer le tenseur des 
efforts de frottement en fonction du tenseur des 
déformations, le système (1 &, b) n’est pas complet. 
Il faut lui ajouter l’équation qui exprime le premier 
principe de la thermodynamique et les équations 
d'état du fluide. Ces équations sont les suivantes : 


CR | 
dt dé (2) 
TA CL PE) Eos RP), 


Q étant la quantité de chaleur gagnée ou perdue 
par les particules, E; l'énergie interne rapportée à 
l'unité de masse et 7° la température absolue. L’opé- 


rateur 


d ) es 
== +VN (3) 
dé ol 


est l'opérateur des variations individuelles (varia- 
tions sur les éléments matériels). Si le fluide est 


un gaz parfait (au sens thermodynamique), on- 
a E;— C;T et les équations (2) prennent la forme 
HET 
CRE. | 
dt dé (3) 
aP =RT | 


C, étant la chaleur spécifique sous volume constant 


et À la constante du gaz. 
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Dans le cas du fluide compressible où l’on peut 
Le 


exprimer 4 en fonction de V et de coefficients de 
viscosité supposés connus, le système (1) + (2) 
n’est un système complet que si, par exemple, la 
fonction Q est une donnée du problème. Il y a 


alors les sept inconnues VV, V5 V.), PScRErel 
les sept équations (1) + (2). 


2. Position du problème général. Nous 
donnons le nom d'équations intégrales de l’hydro- 
dynamique aux équations qui expriment la solution 


> : 
PÉLOeE M) EP (CR 
RAS) 


du système (1) + (2) en un point M(x, y, z) du 
domaine D(t) occupé par le fluide pour toute valeur 
de {=-{, après un instant initial arbitraire # en 
fonction des valeurs initiales 


RUE 
TRAME EEE 


(4) 


AL, Vs Jr 23 L): 


> 
Vo(z, y; z, Co); PL Z, lo); | (5) 
402, Z7; lo }; To(x, y, 3; lo), (B)o(æ, y; Z; lo) 
dans tout le domaine D({,) et des Valeurs aux limites 
> oŸ )P 
Fr, RE 
(6) 
d2 01 


sur la ou les surfaces frontières S(t) de Do pour 


toutes les valeurs de {=> 44. L'opérateur + 2. désigne 


une dérivation suivant les normales n aux S(t) et 


o? 
2Pe des incon- 
nues jouent le rôle de valeurs aux limites. Il va 
sans dire que les valeurs aux limites ne peuvent pas 
être choisies arbitrairement, car on a la condition 
cinématique 


éventuellement les dérivées secondes 


io 
Ver 


(7) 
sur S(), Vy désignant la composante normale de 
la vitesse des éléments des S(. D'autre part, sur 
une surface libre on a la condition P = o. 

S'il est possible d'établir de telles équations 
intégrales correspondant, sans aucune restriction, au 
système (1) + (2), alors on peut dire que le problème 
général de l’hydrodynamique est résolu. En effet, la 
théorie des équations aux dérivées partielles nous 
apprend qu’il est de l'essence même du sys- 
tème (1) + (2) que sa solution générale se présente 
sous la forme d'équations intégrales. On voit donc 
que la thermo-hydrodynamique classique, dont le 
contenu est exprimé par le système (1) + (2), ne 
peut conduire, dans le cas général, à la solution du 
problème fondamental de la prévision de l’évolution 
d’un milieu fluide, puisque les valeurs aux limites (6) 
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Ne 


sont en général des inconnues pour { > {,. La prévi- 
sion des valeurs aux limites est même la partie la 


plus importante de certains problèmes (le problème : 


atmosphérique par exemple). Pour pouvoir résoudre 
le problème de. la prévision il faut donc ajouter 
aux principes de la thermo-hydrodynamique un 
autre principe permettant en quelque sorte d'éli- 
miner l'influence des valeurs aux limites dans les 
équations intégrales, en les remplaçant par des 
expressions ne faisant intervenir que les conditions 
initiales. C’est ce que permet de faire la théorie 
des perturbations dont il sera question plus loin. 


Cette impossibilité foncière où se trouve la thermo- 
hydrodynamique classique de résoudre le problème 
de la prévision de l’évolution d’un fluide (quand 
les valeurs aux limites ne sont pas des données) est 
d’ailleurs toute naturelle. Il est clair en effet que 
lPévolution d’un fluide ne peut être complètement 
déterminée que si l’on connaît les actions extérieures 
qui représentent les liaisons (mécaniques, thermo- 
dynamiques, électromagnétiques) entre le fluide et 
le reste du monde. Les équations de la thermo- 
hydrodynamique doivent donc être insuffisantes 
pour traiter le problème de la prévision puisqu'elles 
établissent une séparation en quelque sorte arbi- 
traire entre le fluide et le milieu extérieur, et c’est 
précisément en se donnant les valeurs aux limites 
qu’on peut suppléer à l'ignorance des actions exté- 
rieures. Au point de vue des équations de ia thermo- 
hydrodynamique, le phénomène hydrodynamique 
décrit par le système (1) + (2) se comporte comme 
un phénomène aulonome dans “une certaine région 
du monde et les actions extérieures sont repré- 
sentées par les valeurs aux limites. 


‘Pour arriver à résoudre un problème de prévision 
pour un fluide il est alors indiqué de décrire l’in- 
fluence des actions extérieures non pas à l’aide des 
valeurs aux limites (qu’il est en général impossible 
de considérer comme des données) maïs en consi- 
dérant dans le domaine de l’espace occupé par 
le fluide un phénomène hydrodynamique auxiliaire, 
satisfaisant donc à toutes les équations du sys- 
tème (1) + (2) et représentant, en chaque point 
de ce domaine, l'effet direct des actions extérieures 
sur le fluide. Si ce phénomène hydrodynamique auxi- 
liaire (qu’on peut appeler le mouvement entretenu par 
les actions extérieures) peut être considéré comme 
une donnée, alors l’évolution du mouvement total 
(réellement observé) doit devenir complètement 
déterminée. En d’autres térmes, la différence entre 
le mouvement total et le mouvement entretenu (la 
perturbation) ne doit plus dépendre que des condi- 
tions initiales dans tout le domaine. L'intérêt de 
l'introduction du phénomène hydrodynamique auxi- 
liaire réside dans le fait que ce phénomène peut 
souvent être considéré comme quasi permanent 
alors même que les valeurs aux limites relatives au 
‘mouvement total sont très variables (avec le temps). 


PR AR .4 = | 2... LA PC LT JA TL ” EN CS a 
2 2 v- e : É > c' 


Nat 


Nous verrons plus loin comment l'introduction du 
mouvement entretenu permet d'utiliser les équations 
intégrales de l’hydrodynamique pour la prévision 
de la perturbation. 


IT. — Équations intégrales 
de l’hydrodynamique. 


3. Nouvelle forme (tensorielle) des équations 
aux dérivées partielles de l’hydrodynamique. — 
Pour déduire des équations (1) les équations inté- 
grales correspondantes on a dû recourir jusqu'à 
présent à des conditions assez restrictives : mouve- 
ments lents ou à tourbillons négligeables, fluides 
incompressibles, viscosité constante. Les équations 
intégrales obtenues dans ces conditions (équations 
d’'Oseen) sont très importantes, parce qu'il est clair 
que tout essai pour se libérer des conditions restric- 
tives ne peut être qu’une généralisation des résultats 
classiques d’Oseen [1]. 

Pour arriver à des équations intégrales qui soient 
applicables aussi à des fluides compressibles à 
viscosité quelconque et en mouvement tourbil- 
lonnaire quelconque, prenons des coordonnées carté- 
siennes rectangulaires æ*(k — 1, 2, 3) dans le réfé- 
rentiel par rapport auquel on étudie le mouvement 
du fluide et définissons par 


op 


pre == ENS (k, NM T,"2, 3) (8) 
un tenseur V”. En posant 
> _ DU ME D 'ErAE 
AE ANR EE HQRV EE V NY, (9) 


les équations (1 à) s’écrivent 


à > 
_ LAVE À (10) 
d’où l’on déduit x 
14 JAM 
D AP = Tr (11) 


Le système (1) des équations fondamentales de 

lhydrodynamique prend donc la forme tensorielle 
OV?} ; GARE 

a D “AVE rs | 

; (42) 


et nous allons voir qu'il est possible d'écrire les 
équations intégrales qui correspondent rigoureu- 
sement au système (12 a), c’est-à-dire les équations 
qui expriment les grandeurs V’;, en un point quel- 
conque M du domaine D(t) sine par le fluide et 
pour toute valeur de £{ > {,, en fonction des valeurs 
initiales 

(WE )o(æt, to) (13) 
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dans tout le domaine D(£,) et des valeurs aux limites 


ne OP 
J ee , . Pr Ne 


14) 
AOF. ee 


ou, ce qui revient au même, des valeurs aux,.limites 


PAS Los 
RARE Dre ; 15 
De On TON F 52) 


pour {> {, sur la ou les surfaces frontières de D(P). 
Désignons par À, le laplacien bidimensionnel relatif 
à ces frontières 5. On sait que l’on a sur S 


> A ) 
NAN AC >, 


8 
On (9) 
x étant la courbure moyenne de $S. D'autre part, 
en utilisant les formules de Navier-Stokes pour la 
force de frottement, les expressions (9) deviennent 
> > > > > 
A ES aNP he 20 VU PNP 

_ + 

+ (au — L)AV + a(u +2) v(v. ÿ) 
ET = > . 
, VEN) PAYATTES (y. F) METAL) 
u et À étant les deux coefficients de viscosité. On 
voit donc, grâce à (16), que les valeurs aux limites 


- 
de À sont connues quand on connaît les valeurs aux 
limites suivantes : 


ROMANS 0P 

F, Frs) 

on on 7 on sur S(2) (18) 
du. d?. 

PS OAR NOR 


Indépendamment d’ailleurs de la définition phy- 


sique de À par les équations de l’hydrodynamique 
et les formules de Navier-Stokes, les relations (10) 
et (16) montrent que l’on a 

a? Mo 
on? F4 0n 


É 


A = (a+ 


Le 1V 

oV 

Eure S(t 219) 
L ds MN S(£), (19) 


de sorte qu'on peut remplacer les valeurs aux 
limites (18) par les valeurs aux limites suivantes : 


SES ot 
7 se “ sur S(4) (20) 


[remarquons que si les frontières sont invariables, 
= 


> 
les valeurs de = sur S({) se déduisent des valeurs 


= . 
aux limites de V et l’on a alors simplement 


= D; 
Se 


n° de (0) 


sur $, 


comme valeurs aux limites du problème]. 
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4, PUS adjointes du Siteine (12). - _. En 
désignant par T la variable temps, nous écrirons 
ces ‘équations sous la forme 


DE FA / dx? 
L. FER: à 


d. 


W# étant un tenseur symétrique qui sera déterminé 
ci-dessous. Pour trouver la solution élémentaire du 
système adjoint posons 


(23) 


x! étant les coordonnées cartésiennes rectangulaires 
d’un point M' de D(r) et x, les coordonnées d’un 
autre point M de 1)(r), pôle de la solution élémen- 
taire. Considérons la fonction 


r= MR mi TE (a mi +(R = ri, 


D(r, 4 1)= — e US) (24) 
NES 
qui satisfait comme on sait à l'équation 
Fee. =D. ::(28) 
Formons la fonction » 
Ÿ—= fo dE, (®) 


7: 


(2 


r' étant une constante indépendante der Elle 
satisfait à l’équation 


Considérons alors le tenseur symétrique 


GE AE CE . ET SA PEUR 
W (AW) 04 + UT 
qu'on peut appeler tenseur d'Oseen et qui satisfait 
à l'équation (22 b). Les équations (22 a) seront 
aussi satisfaites par W}, grâce à (27), en définissant 
le vecteur B# par la relation 


(29) 


Par suite de (27) on voit que ce vecteur satisfait à la 
condition 


ABk— 0, (30) 
5. Équations intégrales pour le tenseur V”. 

— Considérons les équations (12 a) où nous désignons 

maintenant par 7 la variable temporelle : 


oV 7 7 dAM 

RENE cn à 
Multiplions ces équations par W}, les équations (22 «) 
par V',, additionnons les résultats obtenus et 


DE PHYSIQUE. 


< “k WE Le 
Pt 


une fonction continue, on peut écrire 


sommons par rapport à ira 


es | 
de la symétrie des tenseurs W! et = on obtient 


l'expression 
De OPEPE) x OPE AT = FR AIPE) 
ne = PES + PRE Co 


dont nous allons prendre l'intégrale de Sie 
étendue à tout le domaine D(r) compris entre la ou L 
les surfaces frontières S(r) du fluide et la surface … 


d’une sphère £ de rayon r’ centrée au point M, 
de D(r) 


É AUTANT 
D(H 2 D(7 


LS £ Cri JAM Fe p} DE) do. 
D(T) ER 


Oxh 1210) 


(FHAVR = VRAI) du © 


Les V” étant supposés satisfaire à des conditions de 
régularité convenables, la formule de Green donne 


Eee (WAVE ya AWE) do 


Dr) TRRÉERER 
> Sie mm IT # À sl DE 
L [CRE rat ja 
Ts on on : 
ne FREE 
a 1% Sr Ve ue " 64 | 


où les normales n à S(r) sont dirigées vers lntéieus 
de D(r). D'autre part, les conditions (22 b) et (30) 
permettent d'écrire, en appliquant la Pre de 
Gauss, les transformations & à See 


LE ! = 
de . nx WE A = fn ae de. 
S(t) n 


= 


.. Are = e on 0B; Bis, 
59 D) Pr dœk 


=) pu À as sf pal D des 
S(T) 


les n; étant les composantes covariantes des normales 

unitaires (remarquons qu'il n'y à ici aucune diffé- 

rence entre covariance et contravariance puisqu'il 
s'agit de coordonnées cartésiennes rectangulaires 


dans l’espace euclidien). Finalement, W{i Le étant 


de 2 Grért)de 

Bt: 9 à 
> 2 1 PHP do — WA TAT, ds, 
2 Via S (7) : 


Va désignant la composante normale [positive vers … 
l'intérieur de S(r)] de la vitesse de déplacement des. +4 


ESS er RES RÉ PRÉ LR SES : 
2 N°5. SUR LES ÉQUATIONS INTÉGRALES DE L'HYDRODYNAMIQUE 293 


éléments de %$S. Ces différentes transformations 
donnent donc à (31') la forme 


HAE 
fl (rt pee }+ml nt4n [ac 
or } 
= — = Wiv'} du 
Ce D (T) 


Sfr CPC x die 
S(:) on VE The 


= n xWEAm_E WE V# | ds 


ie pr El as = (rm 
ST 
Intégrons cette expression par rapport à r der ={1, 


à Tr —{, en tenant compte du fait que W s’annule 
pour 7 —/{, d’après sa définition (28) et (24). Il 
vient 


) fe 
fe as ee Les aa enr | ds 


= WT de 


D (4) 
É Are 7h 
af [rs 8) 
lo st Ge ET On 


+ miPtAn PET | ” 


l ÿ : 
— f ds f. es de (rrn  de (32) 
pa st & 


Il faut chercher maintenant la limite de cette 
expression lorsqu'on fait tendre vers zéro le rayon r' 
de la sphère X. Désignons par W° le tenseur défini 
par l'expression (28) de W* en. y remplaçant la 
fonction W par la fonction W — lim, W. Un calcul 
facile basé sur l’expression (28) et tout à fait analogue 
à un calcul bien connu d’Oseen [2], permet de 
démontrer : 


1° que la limite pour r'-> 0 du premier membre 
de (32) est nulle; 

20 que les limites pour ro des intégrales du 
second membre où intervient WÂ s’obtiennent 
simplement en y remplaçant W# par Wf. 
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De plus, un calcul élémentaire basé sur (29) 
donne 
“Cr (æi) PU D, ES) 


0 bd t—<$ 


B;= 
d’où l’on déduit sur È 


2By (ai (ar) D(r', 5) 


or F5 L—7 


En développant V’* en série de Taylor dans le 


voisinage du pôle M(xi, x?, x), on aura donc 


nf dr [rude 
150 


re À v Dr", +) £ 
= — i lun 7 fe (CV) : — dr. * 


Or, on a démontré [2] que l’on a 


frot ( SUD ge = » ar 


lo 


lim 7’ z2h(t) (34) 


> 0 


pour toute fonction continue A(:), de sorte que 


s & [rn re do 
0, 


rx [F'R]C2T, æè, 


La 0): 


: 


D'autre part, la même expression (33) donne 


dB, Br',r) 9 [(æi—(æih) 
on Th TL EET) on Le 


et en tenant compte de (34) on obtient immédia- 


tement 
l p. 
lim f° ds f pi 2Bi js 
> 0y S(r) on 


’ > d Fr—(x; 
= — Vr4 Vn = Bi (és) ds. 
sû on 1 : 


(36) 
Ces différents résultats, introduits dans (32), donnent 
finalement [3] les équations intégrales pour les V” 


BR, 
7 VrX VPN, t) 


ie (WE) (6, F=.t5) (PE) (WP, to) dw 


ov'? owW# 
k PA Pi ü 
. af LSE DE — VA Fe) 


+ WE (nr AM EVA, | ds 
= Ver | Fri d [Ex | ds. ( 31) 
st on à 


Elles expriment les V*, en un point quelconque 
M{(x') du domaine D(t) pour tout instant £ > {, en 
fonction des conditions initiales (13) et des Cond. 
tions aux limites (14) ou, ce qui revient au même, 
des conditions aux limites (20). On peut montrer 
inversement qu'étant données des conditions initiales 
et aux limites continues, les fonctions (37) satisfont 
aux équations (12 a). 


. Équations intégrales pour le vecteur Ÿ. 
+ ad déterminé, par les équations intégrales (7, 
le tenseur V'{ pour un fluide quelconque, la déter- 


{= 
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mination du champ de vitesse devient un Sroblème 
> . 
classique. En effet, le rotationnel rot V et la diver- 


ce « 14 
gence © — div V de la vitesse sont donnés par 


>\i 
= TA ri 
(rar) = VA = Vi, 
(ë, 7, k = permutation circulaire dé 1; 2; 8) 


et 
OÙ VE, 


et sont donc déterminés par les équations géné- 
rales (37). Du rotationnel et de la divergence on 


déduit la vitesse Ÿ elle-même quand on connaît 
sa composante normale sur les frontières S. Cette 
donnée aux limites fait précisément partie du 
système de conditions aux limites qui intervient 
dans les équations intégrales (37) pour les V#. Nous 
croyons inutile de reproduire ici le sue classique 


en HQE du potentiel, qui donne ÿ en fonction 


de rot ÿ, 6 et V,. 


7. Le problème général aux limites pour les 
fonctions continues spatio-temporelles. — La 
détermination des équations intégrales des V7, 
dont nous nous sommes occupé dans les para- 
graphes précédents, est en réalité purement phéno- 


5 
ménologique, c’est-à-dire indépendante du fait que V 
est un champ de vitesse des particules d’un milieu 
fluide. En effet, pour une fonction continue spatio- 
temporelle quelconque V(x',{) ayant des dérivées 
premières et secondes on peut écrire l’équation 


oV 


AP = AlEL (38) 


4 étant une constante arbitraire, et ceci d’une 
manière complètement indépendante de la définition 
physique éventuelle du second membre A(æ, ?) 
de (38), par exemple la définition (17) par les équa- 
tions de l’hydrodynamique. Si V possède des dérivées 
du troisième ordre, on déduit de (38) 

OV ; A 


De AE or 


Si la fonction continue V a des dérivées premières 
et secondes continues et des dérivées du troisième 
ordre bornées et intégrables, alors les équations 
intégrales (37) s'appliquent sans aucune modifi- 
oV 
cation aux V; = -——. 
dx 


nous pouvons en déduire 


ENT OS TE 
sr nie 


et une application d’une formule classique donnera 


Ayant déterminé les Y.,, 


LEE 
de: Ÿ t) 


EX | ; CPS ON PRE To EU dard ee AT ES AE de. 
BALL PMP DE RENE Se Re 


en chaque oi M de Do 


ir VOL = à GO, M; t)AV(M', dom, 


*D(1) 
. VO, Ur M; GW, M5 t) 3, 


G(M',M;t étant la fonction de Green du 
domaine D(f), de pôle M, et s’annulant sur les 
frontières S(). 


Les équations intégrales (37) bte done 


de résoudre phénoménologiquement le problème 
général aux limites pour toute fonction continue 
spatio-temporelle V ayant des dérivées premières 
et secondes continues et des dérivées du troisième 


ordre bornées et intégrables dans le domaine d’exis- 


tence D(f de V, c’est-à-dire-le problème de la déter- 
mination de V, en chaque point de D(f) à tout 


instant #>1{,, en fonction des valeurs ïini- 
tiales V(x, {,) dans tout le D({,) et des valeurs aux 
limites 


> nd ar 
SH an dre 


; ; ÿ 0 
Si l’on exigeait que 7 
ne fasse pas partie des valeurs aux limites, alors il 
faudrait recourir, si S dépend de f, à une définition 
de À indépendante de (38), comme la définition (17); 


les valeurs aux limites seraient alors 


sur les S(f) pour £ = {,. 


PU eue 
on ti One 
mais le problème perdrait ainsi un peu de sa pureté 
phénoménologique. 


8. Détermination de la densité et de la pres- 
sion. — D'accord avec les résultats précédents, 


cette détermination, en fonction des valeurs initiales 
et aux limites, 
quement, sans utiliser les équations de l’hydro- 
dynamique. Mais la détermination de la densité 
résulte aussi très simplement de l'équation de conti- 
nuité (1 b), quand on a déterminé préalablement le 
ne de vitesse par les équations (37). Connais- 


sant Ÿ on peut déduire des équations 


dæi 


= Pi(ai; at, a, à), 


les coordonnées a’ des particules en fonction du 
temps et de leurs coordonnées initiales al (coordon- 
nées de Lagrange) 


Lu tata dit). 


Alors, l’équation de continuité donne 


t 
(O] (at, a?, at, 0) dt 
a( ai, a, a, t) = go(at, d?, as, to) e*h ; 


peut être faite phénoménologi- - 


À, » 
one 


pi 


PASS LS CE AT 7e . OUR enr e QUAI RÉ, 
CNE 7 : < 
ADS 2 


CLEAN ' re F- ERP 3 RR, 3 


N° 5. 


l’intégration étant effectuée sur les trajectoires des 


particules. 


Après avoir déterminé Ÿ et x, le calcul hydro- 
dynamique de la pression peut être fait par l’équation 
du mouvement (1 a) quand on connaît les coefficients 
de viscosité qui interviennent dans l’expression des 
forces de frottement en fonction des dérivées 


spatiales de 2 (par exemple dans les formules de 


. Navier-Stokes). 


III. — Application 
à la théorie de perturbations. 


9. Examen d'un cas particulier. — Considérons 
un milieu fluide adhérant à ses surfaces frontières S (?). 


2 
La composante tangentielle de V sur S() s’annule 
et la composante normale est donnée, comme 


>- 2 
toujours, par V.n— V,. Envisageons un phéno- 
mène hydrodynamique auxiliaire décrit par des 


champs de vitesse (vx ), de pression (P*), de volume 


spécifique (x*), etc. La relation (9) permet de définir 


pour ce mouvement un vecteur A* et en partant 
de conditions initiales et aux limites arbitraires 
les équations (37) donnent des grandeurs (V*)" 
satisfaisant aux équations (12 a). Supposons que lés 
valeurs aux limites satisfont aux conditions 


ce 
or ss av 
on æ on’ £ S 36 à 
‘4 : sur. S(t). (39) 
HR re Li 
don? RS, #9ns ve 


En tenant compte de la relation (16), de l'équation 


S 
du mouvement (1 a), de la définition (9) de A* et 
de la condition d’adhérence sur les S(), on. voit 
facilement que les conditions (39) entraînent 


ce _ 
T 


F*— GÉNP=FS — sVP sur S(t). (40) 


Si le mouvement réellement observé (Y, P, à) est 


en régime de Poiseuille (F —a«yP) dans le voisi- 
nage des parois S(), les conditions (40) exigent 
donc que le mouvement auxiliaire soit aussi en 
régime de Poiseuille dans le voisinage des parois. 


> 
Définissons une perturbation (o, P; "à par 


SRE RS 


v=V—F"; p=P— Pr; = 4 — a" 


el posons 


Les équations (37) étant linéaires, en les appliquant 
successivement au mouvement réellement observé et 
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au mouvement auxiliaire on trouve simplement pour 
la perturbation, grâce aux conditions (39) et à la 
condition d’adhérence 

87 /— & 
s V TX ve —= 


D (4) 


) WÉ(E, 5 =) [v'E] (M, b)do, =. (A) 
On voit donc que la perturbation, définie par (41) 
par rapport au phénomène auxiliaire envisagé ici, 
a une évolution déterminée complètement par les seules 
conditions initiales. 


10. Les équations intégrales de l'hydrody- 
namique et le principe fondamental de la 
théorie des perturbations spontanées. -- Dans 
cette théorie [4], on décompose le mouvement 
réellement observé dans un fluide en deux parties : 


19 le mouvement entrelenu; 


20 la perturbation spontanée. 


En désignant par des astérisques les variables 
relatives au mouvement entretenu, la perturbation 


spontanée sera définie par les relations (41). Le 


mouvement entretenu correspond à l'effet des 
actions extérieures qui agissent sur le fluide et repré- 
sente ces actions par des champs bien définis en 
chaque point du fluide au lieu de les représenter 
par des valeurs aux limites (cf. $2). Le principe 
fondamental de la théorie des perturbations spon- 
tanées nous apprend que la connaissance du mouve- 
ment entretenu est nécessaire et suffisante pour déler- 
miner hydrodynamiquement l’évolution de la pertur- 
balion. En d’autres termes, on peut toujours définir, 
dans un fluide soumis à des actions extérieures, un 
phénomène hydrodynamique auxiliaire satisfaisant 
aux équations de l’hydrodynamique, et dont la 
connaissance à tout instant permet de prévoir 
évolution du fluide. Ceci signifie que la pertur- 
bation spontanée, définie par (41) par rapport au 
mouvement entretenu de la théorie des perlurbalions, 
ne dépend pas, dans son évolution, des valeurs aux 
limites, et doit être complètement déterminée par les 


‘.seules conditions initiales. Les équations (42), que 


nous avons déduites des équations intégrales (37) 
de l’hydrodynamique, en admettant un phénomène 
hydrodynamique auxiliaire adhérant aux surfaces 
frontières et satisfaisant aux conditions (39), peuvent 
donc être écrites en général. D’après ces équa- 
tions (42), on voit que si la perturbation spontanée 
est nulle à un instant donné dans tout le fluide, 
elle sera nulle à tout instant ultérieur, à la condition 
évidemment que les conditions de régularité nécessaires 
pour déduire les équations intégrales (37) soient 


‘valables à tout instant t > {5 


Si les actions extérieures sont constantes, le 
mouvement entretenu correspondant sera permanent, 
et alors toute l’évolution du fluide est une évolution 
de la perturbation spontanée. Elle résulte donc, 
dans ce cas important, uniquement des variations 
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temporelles des intégrales 


[ WA(E, z=t)[v#] OP, )du;,, 


/D(1) 


c’est-à-dire en dernière analysé des variations tempo- 
relles des fonctions ® et W qui interviennent dans la 
définition (28) du tenseur d’Oseen. De cette défi- 
nition et des définitions (24) et (26) des fonctions D 
et W on déduit sans difficulté pour WÂ(4 5 —1,) 


AE, 5 = 40) 
D pee en 


= o | = DER =) 
+ > [si = (œil 2i) (aie “.] Dr, 5 = to). 


2% r? lt; 


Par suite de propriétés évidentes des fonctions ® 
et W on obtient immédiatement le résultat 
lim  WE(t1, = = h)=0, 
({— to) + & 

qui entraîne, d’après (42), la conséquence impor- 
tante que la perturbation spontanée tend toujours 
vers zéro à la limite (pour 1{—1,-+ æ). Ce résultat 
n’est évidemment valable que si les V? satisfont 
aux conditions de régularité qui permettent d'écrire 
les équations intégrales (37), c’est-à-dire s'ils sont 
continus, ont des dérivées premières continues et des 
dérivées secondes bornées et intégrables. Supposons 
que les valeurs de { pour lesquelles ces conditions 
de régularité ne sont pas satisfaites partout dans 


n° + 


“ 


le domaine D(1) tone un ensemble désombrable £ 
d’instants que nous désignerons par £,(n — 1,9, …, ©). 
Les équations intégrales (37) et les équations corres- 
pondantes (42) pour la perturbation ne pourront 
alors s'appliquer que pour {,<t<14 et dans les 
intervalles 4, < {<1,.,;et donneront pour chacun de 
ces intervalles les valeurs de la perturbation en 
fonction des conditions « initiales » v'i(M', d) 


EPS 
Vrxv? (M1) 


= WE (4; 
D (tn) 


LE) [ot] CWr, In) do, 


de sorte qu’il faut opérer chaque fois avec de 
nouvelles conditions initiales; c’est précisément ce 
fait qui correspond au non-amortissement de Ja 
perturbation pour {->, les valeurs absolues des 
composantes du tenseur d’Oseen W/({; 7 —=1,) 
augmentant avec { dans le voisinage de { — {, avant 
de tendre vers zéro pour {— . Remarquons que 
les singularités de la perturbation qui ont lieu aux 
instants {, sont évidemment le résultat de singu- 
larités équivalentes dans le mouvement entretenu, 
car d’après le principe central de la théorie des pertur- 
bations l’évolution de la perturbation est complè- 
tement déterminée par la connaissance du mouve- 
ment entretenu. Remarquons enfin. qu'on peut 
avoir un mouvement total sans singularités par là 
compensation des singularités du mouvement entre- … 
tenu et de la perturbation. 


Manuscrit reçu en juillet 1940. 
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SUR LES RAYONNEMENTS 6 ET y DE L'ACTINIUM ET DE L'ACTINIUM X 


Par M. LECOIN, M. PEREY, M. RIOU et J. TEILLAC. 
Institut du Radium, Paris. 


Sommaire. — Le rayonnement y de l’actininm et de l’actinium K a été étudié au moyen d’un 
compteur à xénon. La désintégration de l’actinium s'accompagne de l’émission d’un rayonnement y 
de 37 keV (2,2 photons par 1 000 désintégrations), d’un rayonnement de 15 keV (5 photons par 1 00-désin- 
tégrations) caractérisé comme résultant de la conversion du rayonnement de 37 keV dans le niveau L du 
radioactinium. L'existence d’un rayonnement y de 300 keV (2 photons par r 000 désintégrations) 
reste douteuse. 

Lors de sa désintégration l’actinium X émet : un raYonnement y de 350 keV (6 photons par 100 désin- 
tégrations), un rayonnement y de 48,6 keV (27 photons par 100 désintégrations) et un rayonnement L 
(25 photons par 100 désintégrations) . 

L'étude du rayonnement 8 émis par l’actinium a été faite au moyen d’une chambre Wilson à basse 
pression. Elle a montré la’ présence d'un rayonnement électronique de faible intensité (17 électrons par 
100 désintégrations, au maximum). La répartition spectrale et l'intensité du rayonnement correspondent 
approximativement aux électrons émis lors de la conversion du rayonnement y de 37keV dans les 


L 
niveaux L et M. L'énergie mise en jeu lors de la transition Ac —> RAc semble donc être au moins 
de 37 keV. Dans 85 pour r00 des désintégrations de l’actinium, il devrait done être émis un spectre 


continu & dont la limite supérieure serait de 37 kéV, ce qui semble exclus d’après ces expériences. 


: ion. — Une première étude des 

1. Introduction U tude d 

rayonnements & et y émis lors de la désintégration 
. 


de l’actinium Ac > RAc a mis en évidence l’exis- 


tence d’un rayonnement @ peu important : 8 élec- 


trons de 10 à 20 keV sont émis pour 100 désinté- 
grations d’actinium [1]; un rayonnement y de 37 keV 
- de très faible intensité (0,27 photons par 100 désin- 
tégrations), un rayonnement X ou y de 15 keV 


très important (37 photons pour 100 désintégrations) | 


sont également émis [2]. Il s’est avéré impossible 
d'établir avec ces résultats un schéma de désinté- 
gration satisfaisant pour l’actinium et de discuter 
utilement l’existence de phénomènes d'interaction 
entre le noyau et l'atmosphère électronique de 
Vatome pouvant expliquer Flabsence d’électrons 
nucléaires (création e). : : 

Nous avons repris l'étude de ces rayonnements 
en bénéficiant d’une part de nouvelles -techniques 
de préparation de sources d’actinium beaucoup 
plus concentrées en actinium que celles utilisées 
jusqu'alors et d’autre part, des avantages qu’offrent 
les compteurs de Geiger-Muller par rapport aux 


chambres d’ionisation pour étude des rayon- 
nements X et peu intenses. “| 
2. Préparation des sources. — Dans toutes 


ces expériences il importe d'utiliser de lactinium 
privé de tous ses dérivés. Nous avons déjà exposé 
par ailleurs la méthode utilisée pour arriver à un 
degré de pureté suffisant : cette méthode comporte 
de très nombreuses opérations chimiques. Pour 
effectuer correctement ces opérations il est néces- 
saire d’opérer en présence de quantités relativement 


importantes de lanthane (quelques centigrammes 
de lanthane pour : mC d’actinium). Même dans ces 
conditions il y a déjà perte d’actinium lors de 
Vélimination de l’actinium X : cette élimination 
se ‘fait en ajoutant dans une solution d’actinium 
et de La,O, un peu de baryum. On précipite La, 0; 
par l’ammoniaque  décarbonatée. L'expérience 
montre que dans ces conditions une partie de l’acti- 
nium reste en solution, alors qu’une faible pro- 


‘portion de lactinium X présent est précipité sur 
Je lanthane. Pour se débarrasser complètement 


de l’actinium X, il est nécessaire de répéter plusieurs 
fois l’opération et la teneur en actinium de la pré- 
paration décroît à chaque tour de purification, cette 
perte d’actinium étant d'autant plus grande que la 
concentration en actinium est plus élevée. Il est 
donc impossible d'envisager la purification chimique 
de préparations d’actinium plus concentrées. Mais la 
présence de telles quantités de lanthane apporté 
une gêne importante dans létude des rayon- 
nements X et y très mous. En particulier la déter- 
mination de l'intensité absolue des rayonnements L 
de 15 keV sur de telles préparations ne peut s’effec- 
tuer qu'à l’aide de corrections importantes. Dans 
ces nouvelles expériences nous avons utilisé une 
nouvelle méthode très ‘rapide de fractionnement 
des terres rares imaginée par l’un de nous [3] : elle 
consiste en des précipitations dans des conditions 
bien déterminées de concentrations, des terres 
rares actinifères par l’'ammoniaque en excès, en 
présence de quantités variables d’acétate d’ammo- 
niaque. 

En opérant à la concentration de 70 g d’acétate 


d'’ammoniaque pour 1g de lanthane actinifère 
on obtient, pour certaines concentrations d’ammo- 
niaque, un précipité contenant 60 pour 100 du lan- 
thane initial rigoureusement inactif; il reste en solu- 
tion {o.pour 100 du lanthane avec tout l’actinium. 
En 2 h, après un certain nombre de tours de purifi- 
cation nous avons pu obtenir à partir de lanthane 
actinifère soigneusement purifié des dérivés de 
l’actinium, un produit contenant 0,5 mC d’acti- 
nium pour moins de 0,1 mg de La,O, (masse de : mC 
d’actinium — o,o12 mg). La quantité de radioacti- 
nium et d’actinium X formée pendant ces deux 
heures est inférieure à 1 pour 100. Cette concen- 
tration en actinium est la plus forte qui ait été 
réalisée jusqu'à présent. 

Les sources utilisées dans les expériences sur le 
rayonnement > se présentaient sous la forme d’oxyde 
de lanthane actinifère calciné placé dans de petites 
coupelles de pyrex, hermétiquement closes par des 
lamelles de mica de 2 à 3 mg : cm? Are 

3. Étude du rayonnement X et de l’actinium 
et de l’actinium K. — 31. Dispositif expérimental. 
— Pour l'étude du rayonnement ; mou (inférieur 
à 1ookeV) nous avons utilisé un compteur de 
Geiger-Muller à efficacité élevée pour de tels rayon- 
nements. Ce compteur, déjà employé par l’un d’entre 
nous pour d’autres études [4], est un cone en 
bout fermé par une feuille d'aluminium de 5 mg : em? 
rempli de xénon sous une pression de 1/4.em de 
mercure et d'alcool sous une pression de 2m. On 
place les sources à une distance variant de r0 à 20 cm 
de la fenêtre, suivant l’intensité des sources réalisées. 
Le rayonnement est canalisé par deux diaphragmes 
de plomb et les rayons sont éliminés par un champ 
magnétique de 1 700 gauss produit par un aimant 
permanent. 

Pour les études des rayonnements 7 be 
supérieure à 100 keV, nous avons utilisé un compteur 
de Geïger-Muller ordinaire en magnésium de o,1 mm 
d'épaisseur entouré d’une feuille de plomb de 
170 mg :em? servant de radiateur. Dans certaines 
expériences nous avons utilisé également un 
compteur analogue tapissé intérieurement d’une 
feuille d’or de 19 mg : cm2. 


32. Rayonnement de l’aclinium et de l'actinium K° 
— L'étude du rayonnement + de l’actinium initia- 
lement privé de ses dérivés est rendue complexe 
par la croissance des dérivés d’actinium au cours 
de l'expérience. La durée des purifications et concen- 
trations diverses de l’actinium variant de 2 à 3 h, 
on peut admettre que l’actinium est en équilibre 
avec lactinium X au cours de toutes ces expériences. 
Comme nous le verrons au paragraphe 6, le nombre 
de désintégrations de radioactinium et d’actinium X 
formées pendant ce temps est inférieur à 1 pour 100 
du nombre de désintégrations d’actinium, et les 
corrections à effectuer pour connaître le rayon- 
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nement propre à Ac + Ac X sont peu importantes, 
tout au moins pour les composantes intenses. La 
figure 1 représente une courbe d’absorptiondans 
le cuivre obtenue avec le compteur à xénon et une 
source d’actinium de :,9mC dans {mg d'oxyde 
de lanthane, 

On a commencé les mesures 3 h après la fin des 
purifications. La durée des mesures a été de 1h 
environ. 

D'autres sources d’ intensité variant de 0,5 à 4 mC 
(poids de matière o,1 à 20 mg de La,O,) ont donné 
des courbes identiques, mises à part les différences 
sur les intensités des rayonnements mous pour 


lesquels on doit tenir compte de l’absorption dans la 


matière constituant la source. 


Le dispositif expérimental permettait d'étudier 


en même temps que l’absorption du rayonnement 
mou dans le cuivre au moyen du compteur à xénon, 
l'absorption dans le plomb au moyen du compteur 
à feuille d’or, le rayonnement y ayant traversé le 
fond de la coupelle de pyrex. 


Ing À log 1, - 
1 4,8 
3.6 4,67 
3,4 a+ 
8,2 nt 
8,0 l | 
2,8 3,8 
Q 
sr 860 $5 100 180 : 
mg/cm 
2,4 ee : 
20 AB LESE ST TU 
OCR RE OU TND M/S ne 


Masse superficielle des écrans de Cu en g/cm? - 


Fig: 1.— Absorption des rayons y de Ac -- AcK 
dans le cuivre. 
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L'analyse graphique de toutes les courbes d’ab- 
sorption obtenues donne des résultats cohérents 
et permet de GÉtnEseS les quatre composantes 
suivantes : ) 


Coefficient d'absorption : 
Dans Cu (en em? :g)... 0,14 
Dans Pb (en em? :g)… 
Énergie (en keV)...... Ë 


Les courbes d'absorption dans le cuivre semblent . 
donner une indication peu certaine de l'existence 
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I. CALE 


d’un groupe d’énergie voisine de 90 keV. D'autre 
part, la variation notable d'intensité dans le temps 
pour les composantes de 240 keV ne permet pas 
a priori d'attribuer cette composante à (Ac + AcK) 
nous reviendrons sur cette question plus loin en 
étudiant les courbes d'évolution des divers rayon- 
nements. De plus, le rayonnement pénétrant semble 
être complexe et les courbes d'absorption dans le 


| plomb obtenues avec le compteur à xénon et le 


compteur à feuille d’or indiquent qu'il existe au 
moins deux groupes pénétrants de coefficients 
d'absorption (en cm?:g) 0,22 0,02 (400 keV) 
0,9 © 0,2 (220 keV) et peut être un groupe inter- 
médiaire isolé avec le compteur à feuille d’or et de 
coeflicient d'absorption 0,3 (340 keV). 

Pour préciser la valeur des énergies des raies les 
moins pénétrantes, nous avons employé la méthode 
d'absorption sélective, le détecteur étant le compteur 
à xénon. Nous n'avons pu isoler l’hypothétique 
rayonnement de ookeV par absorption dans les 
éléments Pt, Au, Fo, TI, Pb et’:Bi (9): Par contre 
l'absorption dans les éléments allant de liode 
(Z=53) au tungstène (Z —74) nous a permis 
de caractériser une raie encadrée par les discon- 
tinuités X de ,,1 et ,çBa qui correspond à la raie 
de 36,8 + 0,7 keV mieux isolée précédemment [1] 
et une raie nouvelle encadrée par les discontinuités K 
de Eu et Cd, ce qui donne une énergie égale 
AMIS. Det 61,0 keV en filtrant le rayonnement 
par 99 mg : em? de cuivre et en ajoutant un écran 
d'oxyde de lutécium, puis un écran d’oxyde de 
gadolinium (12o0mg:cem? en élément), le nombre 
d’impulsions par minute passe de 5ro à 790. Cette 
raie s'identifie à la composante de 60 keV déduite 


des courbes d'absorption dans le cuivre et le plomb 


si l’on tient compte de l’imprécision que comporte 
l’analyse précédemment faite et du fait que l’absorp- 
tion dans les éléments compris entre &4Gd et :;W 
n'indique aucune raie d'intensité notable entre 5o 
et 50 keV. 

L'absorption dans les éléments allant de ,,Fe 
à Cd permet d’analÿser la composante de 14 keV; 
les écrans sont constitués soit de feuilles métalliques, 
soit de carbonates ou oxydes dilués dans la paraffine 
et ont des masses superficielles voisines de ro mg : em?. 
Les raies isolées s’identifient aux raies L caracté- 
ristiques du thorium comme le montre le tableau 
ci-dessous, où l’on a indiqué les longueurs d’onde 
en U.X. des discontinuités X entre lesquelles se 
placent les raies observées et les raies L du thorium 
correspondantes dont la longueur d'onde est expri- 

; s $ 12 ,4.10° 
mée en U. X. Aux.) 7 EtkeV) |: 

. D'autre part nous avons cherché s’il existait des 
rayonnements d'énergie inférieure à rokeV en 
traçant les courbes d'absorption dans l’aluminium 


() Si ce rayonnement existe, il est beaucoup moins intense 
que les autres. 
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avec ou sans filtre de cuivre de 8,2 mg : cm? qui sert à 
renforcer les raies d'énergie inférieure à celle de la 
discontinuité X du cuivre (9 keV) : les deux courbes 
d'absorption sont identiques et donnent un coeffi- 
cient d'absorption égal à 8 (15 keV), qui correspond 
bien au rayonnement L : nous pouvons donc dire 
qu'il n'existe pas de rayonnnement entre 5 et ro keV 
d'intensité comparable à celle du rayonnement J. 
Les rayonnements plus absorbables comme le rayon- 
nement M ne pouvaient être observés avec notre dis- 
posili]. 


O7 (34Se) (Lai. - 954 
DIS (35 Br) b Las 966 
See (37Rb) °ZLf, 591 
TOSPRELE (38Sr) CLS 592 
FOBAR er (38 Sr) Lie =764 
687: 5 (4oZr) L Lx: 753 
652 (41 Nb) fs Le 631 
Giga (42Mo) ULY3! 626 


4. Étude du rayonnement } de l’actinium X. 
— Malgré la courte période de l’actinium Æ (21 mn) 
nous avons pu, suivant une technique déjà exposée, 


log Î log Î 
2,9 2,8F [ Pb 
28/[ PE | 
2,71 1,9L AC 
2,6 À 1,7 | U 
DS 1,5 S 
à: 0 
2,#P + 
QE Rd Cu 
2SPE & Br2e | 
= = 83 Lg ; F 
ES ER 
2,0 £ i 1 L 1 1 a L 
0 50 100 150 200 250 300 850 400 


Masse superficielle des écrans de Cu en mg/cm° 


Fig. 2. — Absorption du rayonnement y de l’actinium K 


dans Cu et Pb. 


étudier le rayonnement y en préparant une dizaine 
de sources et en disposant de part et d’autres de 
chacune d’elles deux compteurs : du côté de la feuille 
de mica le compteur à xénon pour étudier les rayon- 
nements mous, du côté du culot de la coupelle de 
verre, le compteur à feuille d’or ou celui à feuille 
de plomb pour les rayonnements pénétrants. Toutes 
les mesures effectuées ont été ramenées à la valeur 
qu'elles auraient eu à l'instant zéro. 
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La courbe d'absorption dans le cuivre qui est - 


représentée sur la figure 2, permet de distinguer 
nettement deux composantes de coeflicients d’ab- 
sorption 2,5 cm : g (56 keV) et 83 cm : g (14,5 keV). 
Nous avons pu isoler les raies correspondantes : 
le saut d'absorption entre les éléments Eu et Gd 
permet de retrouver la raie de 48,6 + 1,6 keV déjà 
observée avec Ac+AcX; labsorption dans les 
éléments yAsSsSe et ,45r entre lesquels se placent 
les raies de la composante de 14,5 keV indique 
que celle-ci se compose très vraisemblablement 
des raies L caractéristiques du radium, les courbes 
ayant un aspect voisin de celles obtenues par ail- 
leurs [4] avec plus de précision pour l'ionium et le 
radioactinium qui émettent tous deux le rayonne- 
ment L du radium. 

La courbe d'absorption dans le plomb (fig 2) 


AcKLE 

a7 NB 
Ac FAC AE NENC AE EAEAC 

S ÿ 

a 


x R Ke ; 


Le nombre de désintégration AN de chacun des 
trois corps principaux Ac, RAc, AcX, au bout 
d’un temps { à partir d’un nombre d’atomes d’ac- 
tinium N complètement privé de ses dérivés sera 
donc, les constantes radioactives étant exprimées 
en 7-1: à 


AC TE==>20:000 14; 

AC KT NT, 

RAc— 0,0366, 

ACA= 0,0617; 
EAN jaox = (A No)ac XX 0,071(1— 67 FLAT), 
( AN )rae — (À No ac >< 1,0038 { ( e—0,000146 __ p—0,03664 je 
(AN aex = (X No)dac[Tr,006r2(e0,0001F 1, 459 et: 08667 


+ 1,426 e:0,0617/) 
+ 0,011(1,0016 ev;0001#4 
+ 0,0613014 e#7,#7t 
— 1,0029 e—0,06174 )]. 


Le deuxième terme dans l'expression de (AN)ax 
correspond à la formation d’actinium XK par. l’inter- 
médiaire d’actinium ÆX, terme non négligeable au 
bout de quelques heures. Les courbes de formation 
sont données sur la figure 3. 

On peut admettre, étant donné la courte période 
du dépôt actif devant celle de l’actinium X, que 
celui-ci est constamment en équilibre avec le dépôt 
actif, si l’on se trouve en vase clos comme c’est 
le cas dans toutes ces expériences. Le rayonne- 
ment y correspondant est donc celui de Pacti- 
nium X en équilibre avec tous ses dérivés. La 
courbe n° 3 représente la croissance de l’actinium X 
dans l’actinium primitivement pur et permet de 
calculer le nombre d’atomes de radioactinium et 
d’actinium ÀX désintégrés par seconde au bout 
d’un temps donné. 

L'étude de Flabsorption dans différents écrans 
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indique elle aussi par la composante de 0,6 g : cm? 
la raie de 48,6 keV et en plus montre qu'il existe 
un rayonnement plus pénétrant de coefficient 
d'absorption d'environ o,3cm:g (33okeV), ce qui 
est en bon accord avec la valeur plus précise de ce 
coefficient obtenue par Mile Perey et M. Lecoin [5] : 
0,2 CM? :g (400 keV). Il n’est du reste pas impos- 
sible que ce rayonnement énergique soit complexe. 


 d'ééltts 


5. Évolution de l'activité des sources d’ac- 
tinium. — Avant d'affirmer l'attribution des 
différents rayonnements X et 7 précédemment 
mis en évidence à la désintégration de l’actinium, 
il importe au préalable d'évaluer la part que peuvent 
y apporter la présence. des dérivés de l’actinium 
dont l’activité croît avec le temps. Rappelons le 
schéma de désintégration de l’actinium : - 


AUS An de DE 


8 NC C NE 


A Ne 
AcC{ À 
ON 7 


NE 8 is 
NAec'/ P 5 


(- 


© f 
= AB Ac). 


du rayonnement produit par des sources d’acti- 
nium X et de radioactinium d'intensité connue, 
placées dans les mêmes conditions expérimentales 
que les sources d’actinium permettra d'effectuer 
toutes les corrections désirées connaissant l'intensité 
de la source d’actinium et l'instant de la mesure. 


2.074 


. de 0 à 180 mn 
Echelle 


10h 20h 80h . 40h temps 


Fig. 3.— Courbe de croissance de Tactinium X dans l’actinium. 


Nous avons suivi l’évolution de l'activité des 
sources d’actinium au trâvers de différents écrans, 
en utilisant le dispositif déjà décrit pour l’actinium K, 
(utilisation simultanée de deux compteurs, l’un à 
xénon, l’autre ordinaire avec radiateur de plomb). | 
Avec le compteur à xénon on a suivi l’évolution de 
l'activité. au travers de deux écrans de cuivre | 


de 0,542 g : cm? etde 2,52 g : em? laissant passer. 
respectivement la raie de 48,6 keV et le groupe 
d'énergie voisine de »4o keV. Avec les compteurs 


| are nous Phvons utilisé les écrans de plomb 


de 1,53 et 4,4, g : cm?, favorisant respectivement le 


| groupe de 240 et de 400 keV. 


©1100 À 
Ë 10 Cu 0,54 q/cm 
900 Cu 2,52 g/cm 


Cu 0,54 g/cm? Ac 


520 Pb 4,4 g/cm° 
S 200 é 
100 Pb kL /cm?Ac 
0 
2 8 4 5 6 7 8 9 10 11 12 13 14 15 16 1 
Temps en heures 
Fig. 4 — Évolution de l’activité y d’une source d’actinium. 


Nous donnons dans la figure 4 un exemple de ces 
courbes d'évolution pour les différents écrans. Les 
courbes en pointillé correspondent aux courbes 
expérimentales directes (en nombre d’impulsions 


minutes), les courbes en trait plein correspondent : 


aux courbes corrigées de l’effet du à la présence 
d’actinium X, de radioactinium et d’actinium X. 


Les courbes de Pb 1,53 g : em? et Cu 2,52 g : en? 
coïncident presque exactement, nous avons fait 
figurer une seule de ces courbes. On voit que l’ac- 
tivité corrigée reste approximativement constante 
au cours du temps. Deux hypothèses sont alors 
possibles, ou bien ces rayonnements appartiennent 
en propre à l’actinium, ou bien, la purification de 


_ l’actinium n’était pas complète et ce rayonnement 


est produit par des traces d’actinium X ou de 
radioactinium restant présent à la fin de la purifi- 
cation. 


6. Détermination. des intensités absolues. 
Pour obtenir l'intensité des raies y ou L de l’acti- 
nium et de l’actinium X nous devons déterminer 
l'efficacité des compteurs utilisés. Pour le compteur 
à xénon cette détermination a déjà été faite par 


l’un de nous [4]. Pour les compteurs à feuille d’or 
ou de plomb nous avons utilisé les courbes d'efficacité 


en fonction de l'énergie, obtenues par G. À. Renard{[7]. 


D'autre part il faut tenir compte de l'absorption 
des rayons L dans la fenêtre d'entrée du compteur, 
la feuille de mica fermant la coupelle contenant la 
source et surtout dans la matière même de la source 
qui contient des éléments lourds comme le lanthane 
pour les sources d’actinium et le cérium pour celles 
d’actinium X. Pour tenir compte de ce dernier effet, 
il suffit de multiplier l'intensité observée par le 
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facteur 


> 
Re ne pr 


7 
où p. est le coefficient d'absorption du rayonnement 
dans la matière constituant la source et m la masse 
superficielle de cette matière. Nous avons pu du 
reste obtenir une source d’actinium contenant 
moins de o,r mg:em?.de lanthane et plusieurs 
sources d’actinium Æ avec également peu de 
matière : les résultats sont en accord avec ceux 
obtenus avec des sources de quelques milligrammes 
par centimètre carré en faisant la correction ci-dessus. 

La valeur des sources utilisées, en millicuries, 
a été déterminée pour celles d’actinium par le dosage 
en actinium ÆX. La valeur des sources d’actinium K 
a été déterminée par une chambre GOSanon 6 
préalablement étalonnée. 

Les moyennes des résultats obtenus avec diverses 
sources conduisent aux valeurs des intensités abso- 


lues des rayonnements L et y de l’actinium et de 


lPactinium ÆK qui sont reportées dans le tableau ]; 
on a indiqué le compteur utilisé et son efficacité 
(les intensités sont exprimées en nombre de pho- 
tons par 100 désintégrations). 


TaBLeAu |. 


Intensité absolue des rayonnements L et 7 
de l’actinium et de l’actinium K. 


Rayonnement, Compteur. Intensité. 

Actinium. 

Rates LME SES Xe(o,25) (5 x). 107? 

MO S  o 7 eV Re Xe{o,14) (2,240,4)107 

Ravonn. pénétrant (*). Pb(o,ot4) (2Æ 957).10— 
fctinium K. 

ties Li, Reese Tate Xe (0,126) 0,25 0,07 

6 = 1, 6-KeVE = de Xe (0,070) 0,27 0,05 

Mo keV. AE Pb(o,014) 0,06 + 0,03 


(*) Pour faire le calcul on a admis que ce rayonnement 
hypothétique avait une énergie moyenne de 300 keV. 


L'intensité de la raie de 48,6 keV de l’actinium K 
déterminée avec Ac + AcX est (0,25 Æ 0,04).r10 ? 
quanta par désintégration d’actinium, ce qui donne 
comme rapport d'embranchement pour la désin- 
tégration de l’actinium vers lactinium ÆX envi- 
ron 0,009, en bon accord avec la valeur la plus 
précise de Me Perey : o,o11. 

L'accord avec les résultats obtenus précédem- 
ment par M. Lecoin, M°e Perey et Tsien San- 
Tsiang [2] est très bon pour la raie de 36,8 keV de 
l'actinium, puisque ces auteurs obtiennent une 
intensité de (2,7 £1).10 # quanta par désintégration. 
Le désaccord pour la raie de 48,6 keV de l’actinium K 
s'explique par le fait que les auteurs précédents 
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attribuaient à cette raie une énergie trop élevée. 
Il est par contre difficile d'interpréter la divergence 
pour l'intensité des rayonnements L, les intensités 
données par les auteurs étant 5 à 9 fois plus élevées 
que les nôtres. 

Ces résultats avaient été obtenus avec une chambre 
d’ionisation à xénon, le courant d’ionisation étant 
mesuré avec un électromètre Compton. La sensi- 
bilité d’un tel dispositif était environ 100 fois infé- 
rieure au dispositif à compteur précédemment 
décrit. Pour obtenir des courants mesurables il 
avait été nécessaire d'utiliser de très grands angles 
solides : l’élimination du rayonnement fB pouvait 
dans ces conditions être incomplète. D'autre part 
la concentration en actinium des sources utilisées, 
était environ 100 fois plus faible que dans ces 
dernières expériences (8o mg de La, O;, par milli- 
curie d’actinium au lieu de 2 mg :mC. Les cor- 
rections d'absorption du rayonnement L dans la 
matière du produit étaient considérables et ont pu 
facilement être surestimées. 


7. Rayonnement $ de l'actinium. — Nous 
avons repris l’étude de l’émission B de l’actinium 
par la méthode de la chambre de Wilson à basse 
pression. L'avantage de cette méthode est de per- 
mettre la détermination du nombre de rayons £ 
par désintégration, le nombre de désintégrations 
étant déterminé par la numération des + de l’em- 


branchement Ac + AcK; la mesure du parcours 
des rayons dans une chambre à basse pression 
permet d'établir avec une précision suffisante la 
répartition spectrale de’ ces rayons. 

Rappelons rapidement les conditions nécessaires 
pour® obtenir une détermination correcte : ces 
conditions ont été exposées dans le travail antérieur 
de deux d’entre nous; puisque l’évaluation du nombre 
de désintégrations repose sur ‘le nombre d’« émis 
par lactinium, il est nécessaire d'utiliser. de l’acti- 
nium soigneusement purifié de ses dérivés (R Ac, 
Ac X, etc., la limite supérieure permise étant ini- 
tialement au plus de 1/10 ooot). 

D'autre part, il faut tenir compte du fait que le 
temps d'efficacité de la chambre à pression réduite 
n'est pas le même pour les rayons «& et pour les 
rayons $ de faible énergie. On a levé cette indéter- 
mination en ne comptant qué les rayons « présentant 
des rayons 9 bien nets. Enfin si la durée d'utilisation 
des sources excède 3h à compter de la fin des 
purifications il faut tenir compte de la formation 
du radioactinium de l’actinium X dans la source. 

Dans cette nouvelle série de déterminations, 
nous avons essentiellement bénécifié de la nouvelle 
méthode de préparation de sources de La,O, très 
riches en actinium : alors que la concentration 
de l’actinium dans La,O, était de 1 mC pour 8o mg 
de La,O; dans les expériences précédentes, ‘nous 
avons utilisé couramment des préparations corres- 
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pondant à 2 mg-de La,O, pour 1: mC d’actinium. 

Nous avons pu ainsi préparer sur des feuilles 
d'aluminium de o,5 mg:cm? des sources ne pré- 
sentant aucune trace de matière visible : les quantités 
d’actinium à introduire dans la chambre étant 
inférieures à 10 %mC, la masse de La,O,; corres- 
pondante est inférieure à 2.10-*mg répartis sur une 
surface de ro mm? en viron; l’absorption à craindre 
est inférieure à 2.10 ? mg : em? (o,r mm d'air 
environ). D'autre part, la dimension réduite du 
support de source (2 mm x 5 mm) réduit la zone 
de perturbation qu'il apporte : cette zone a des 
dimensions a peu près analogue à cellé du support, 
et il est souvent difficile d'y voir les trajectoires. 

La chambre de Wilson utilisée est la même que 
celle utilisée pour l'étude du spectre 8 du MTh;; 
elle comporte un dispositif d'éclairage par décharge 
condensée dans le xénon. La pression initiale dans la 
chambre était suivant les séries de clichés de 25° 
ou de 15 cm. 


n 
N 
20 


RL Te ES - 
0 10 20 30 kev 
Fig. 5.— Spectre 6 de l'actinium. 


Une quinzaine de séries de clichés ont été eflec- 
tués; le nombre des rayons B observés pour 1 @ 
présentant des rayons d nets a varié suivant les 
séries de 8 à 17; cette variation n’est pas attribuable 
aux conditions initiales de pression; le spectre de 
vitesse est sensiblement le même pour les différentes 
séries de clichés. Il n’y a pas une proportion plus 
importante de rayons de faible énergie '(< 10 keV) 
sur les séries de clichés effectués à 15 cm de pres- 
sion). La variation du nombre de B suivant les 
séries semble plutôt devoir être imputée à des 
traces d’actinium Æ restant dans le produit d’ac- 
tinium à la fin de la purification. Nous avons déjà 
eu l’occasion de signaler les difficultés de purifi- 
cation absolue dans dés sources très concentrées 
en actinium. | 

La détermination du parcours et l'appréciation. 
de l’énergie d’après le parcours a. été effectuée de la 
même façon que pour M Th,. 

La figure 5 donne le spectre correspondant à 
l'ensemble des séries. A: côté d’un maximum très 
accentué à 12 keV il semble exister une probabilité 
d'émission plus forte vers 6 et vers 25 keV. 


8. Discussion des résultats. — Le rayonne-. 


ment y de 37keV émis par Ac doit être fortement 
converti dans la couche L en donnant d’une part 
des électrons de conversion de 17 keV (L,) et d'autre 
part un rayonnement JL, dû au réarrangement du 
cortège après la conversion. En admettant que le 
rayonnement de 15 keV observé soit uniquement 
dû à la conversion dans la couche L,, l'intensité 
des raies électroniques de conversion correspondantes 
ee 
0,40 
désintégration (le rendement de fluorescence dans la 
couche L étant de 0,40 pour le thorium). Ceci est 
en accord satisfaisant avec l'intensité obtenue par la 
méthode Wilson; il semble y avoir désaccord quand 
à la position des raies (12 keV au lieu de 17 calculé), 
malgré l'importance des erreurs expérimentales 
possibles dans l’évaluation de l'énergie des électrons 
par le parcours, cet écart paraît supérieur aux 
erreurs expérimentales. 


devrait être de 


; soit 0,12 + 0,02 électrons par 


D'autre part le rayonnement de 37 keV doit se 


convertir dans les niveaux M de l'atome. Le rapport 


RE : Las ; Ê 
des probabilités de conversion étant d’environ1o; 


les électrons de conversion M ont une énergie 
d'environ 30 keV. La partie du spectre B dont l’éner- 
gie est supérieure à 20 keV pourrait donc assez 
approximativement correspondre à ces électrons M 
de conversion. Il n’a pas été possible, nous l’avons 
précédemment exposé, d’étudier le spectre X M 
avec le compteur à xénon. Nous avons essayé de le 
mettre en évidence directement au moyen de la 
chambre de Wilson. On délimite un faisceau de 
rayons X ou y émis par une source de 3 mC d’acti- 
nium, ce faisceau pénètre dans la chambre de Wilson 
à pression ordinaire au travers d’une fenêtre mince 
de nylon de 2 mg : cm. Avant de pénétrer dans 
la chambre le faisceau passe dans un champ magné- 
tique suffisant pour éliminer le rayonnement $ de 
l’actinium X. Le long du parcours du faisceau dans 
la chambre apparaissent des photo-électrons produits 
par l’absorption dans les niveaux de l'azote ou de 
l'oxygène des différents rayonnements X ou L. Sur 
22 clichés le nombre de photo-électrons dont le par- 
cours est supérieur à : em est de 8 (+ 22 keV), 62 ont 
des parcours compris entre 3 à 6 mm (12 et 17 keV) 
et 46 ont un parcours inférieur à 2 mm (6 à 9 keV). 
Les trajectoires de plus de r cm correspondent proba- 
blement à l’absorption des rayonnements L et y 
de 37 et de 48 keV d’actinium et d’actinium XK. 
Les trajectoires de 3 à 6 mm sont dues à l’absorption 
du rayonnement ZL d’actinium et d’actinium X. 


Pour les trajectoires de longueur inférieures à 2 mm 


la discrimination de l’origine est assez complexe. 
Il peut s'agir de photo-électrons dus à l'absorption 
du rayonnement L émis dans une direction perpen- 
diculaire au plan lumineux éclairé ou de photo- 
électrons dus au rayonnement M, les conditions de 
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l'expérience ne permettent pas de distinguer les 
deux catégories. 5 

En effectuant le calcul du nombre de photons L, qui 
sont nécessaires pour produire les photo-électrons 
de > à 6 mm observés, en tenant compte du &emps 
d'efficacité de la chambre, de l'angle solide du 
faisceau, et de l'absorption du rayonnement L 
dans la matière de la source, dans l’air avant la 
chambre et dans la feuille, on trouve un ordre de 
grandeur satisfaisant. Mais si l’on admet que seu- 
lement la moitié des électrons ayant un parcours 
inférieur à 2 mm sont réellement des photons M, 
en effectuant le même calcul on trouve que le 
nombre de photons M émis par la source devrait être 
égal au nombre de photons Z. 

La partie du spectre électronique ayant une énergie 
inférieure à 7 keV peut être expliquée par les élec- 
trons Auger dus aux transitions M, N —-L. 

En résumé, il est donc possible d'expliquer la 
majeure partie du rayonnement X et du rayon- 
nement $ mis en évidence dans la désintégration 
de l’actinium par la conversion du rayonnement 
de. 37 keV: | 

Compte tenu de cette conversion, la probabilité 
d’excitation nucléaire du niveau de 37keV est 
environ 15 pour 100. Ce niveau est donc excité 


Fe 5 à 
lors de la transition Ac -> RAc puisque la tran- 


La 
sition Ac—>AcX ne représente que 1,1 pour 100 
des désintégrations. 

Dans 8o pour 100 des désintégrations il devrait 
aussi y avoir émission d’un spectre continu B dont 
la limite supérieure serait de 4o keV environ. Les 
expériences qui viennent d’être exposées montrent 
que pour 100 désintégrations, il y a au maximum 
17 rayons & de plus de 3 keV émis, et que sur ces 
17 électrons on peut en attribuer au moins 15 aux 
phénomènes : de conversion du rayonnement 7 
de 37 keV dans les différents niveaux, ou à des phé- 
nomènes connexes (effet Auger). Dans ces condi- 
tions 2 pour 100 au maximum des rayons & repré- 
senteraient la partie du spectre continu comprise 
entre 3 et 4o keV et 95 pour 100 des électrons émis 
lors de la transition & directe Ac —R Ac auraient 
une énergie inférieure à 3 keV. 


9. Conclusion. — Les expériences antérieures 
et nos expériences conduisent aux conclusions 
suivantes : : 


Aclinium K. — Le spectre 8 examiné añtérieure- 
ment à la chambre de Wilson [9] a une énergie maxi- 
mum .de 1,2 + o,1 MeV; mais ces expériences 
n’ont pas permis de préciser le degré d’interdiction, la 
complexité du spectre étant néanmoins certaine [10|. 

L'observation d’un rayonnement + de 330 keV 
de 6 photons par 100 désintégrations confirmerait 
que le spectre est complexe. 

Quant au rayonnement y de 49 keV, il est beaucoup 


‘ 


16. 
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At 

plus intense (0,27 photons par 100 désintégrations) 
et doit être très converti; comme le montre d’ailleurs 
l’observation d’un spectre L de 0,25 photons par 
désintégration. Si l’on admet que ce rayonnement L 
est dû uniquement à la conversion du rayonnement 
de 49 keV, le calcul à partir des intensités observées 
montre que la probabilité d’excitation du niveau 
nucléaire de {9 keV est voisine de 1. Il semble donc 
que le rayonnement de 49 keV soit émis en cascade 
avec le rayonnement £. 


Actinium. — Dans la désintégration Ac + R Ac, 
on observe seulement un rayonnement L, un rayon- 
nement y de 37keV et un spectre électronique 
qui semble attribuable à la conversion de ce rayon- 
nement. 

À partir des résultats obtenus pour les intensités, 
on calcule que la-probabilité d’excitation du niveau 
nucléaire de 37 keV est environ 15 pour r00 (ce 
niveau est donc bien excité lors de la transi- 


[CI = 
ë P . pie x » 
tion Ac_> RAc, puisque la transition Ac + AcK 
ne représente que 1,1 pour 100 des désintégrations). 


L'énergie de la transition Ac de RdAc est donc 
au moins de 37 keV ‘(?). Si l’on admet le chéma 
de désintégration représenté par la figure 6, dans 
85 pour 100 des désintégrations, on devrait alors 
observer un spectre continu $ dont la limite supé- 
rieure serait, au moins de 37 keV. 

En calculant à partir des formules classiques 
la forme d’un spectre ayant une énergie maximum 


ACKQ 

a A PNB 
(2) Ac AC Y. En faisant le bilan énergétique 
PNR AA 


Se p 
dans le cycle on trouve pour la transition Ac —— R'Ac une 
énergie de 100 + 100 keV. 
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de 37 keV, [11] on ou que 20° à 50 pour 100 AE 
électrons nucléaires ont une énergie comprise entre 
10 et 37 keV, or l'expérience montre que seulement 
quelques électrons par ro0 désintégrations PES C4 
être attribués au spectre. 


Ac 


(à 85%) TS 


RAc 


Fig. 6. 


Nous n'avons pas pu conclure à l'existence d’un. 
rayonnement y de 300 keV environ émis par l’ac- 
tinium, en tous cas son intensité est suffisamment * 
faible (0,2 par 100 désintégrations) Fo qu’on 


puisse l’attribuer à la transition Ac > AcK. 


Ce travail confirme l’analogie der désintégrations 


de Ac et MTh, [8], [12] : existence d’une raie y for= M 


tement convertie dans 15% des désintégrations. 
L'énergie de la transition $ dans les 85 % qui restent, 
est donc au moins égale à l’énergie du quantung y, | 
néanmoins on n’a pu mettre en évidence expérimen- | 
talement de tels spectres G. 


Ce travail a été effectué à l’Institut du Radium, 


Laboratoire Curie. Nous remercions 1. Joliot-Curie, 


Directeur de ce laboratoire pour le bienveillant 
intérêt qu’elle nous a témoigné. Nous remercions 
‘également le Centre National de la Recherche Scien- 
tifique qui nous a permis à ou quatre de poursuivre 
ces recherches. 
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THÉORIE DES GAZ ET ÉQUATION D'ÉTAT. IX. 
Par J. DUCLAUX. 


s 
Sommaire. — La théorie des gaz, fondée sur la notion d'association progressive, appliquée d’abord à 
l'oxygène, rend compte très exactement de sa compressibilité dans des limites très étendues. Elle a 
permis de découvrir une propriété nouvelle qui simplifie beaucoup l'étude de cette compressibilité. En 
même temps elle donne une interprétation de la liquéfaction et la relie à la compressibilité, le point de 
liquéfaction étant un point singulier de la courbe de compressibilité, De cette manière il suffit finalement 
d'une seule courbe pour représenter la compressibilité el la liquéfaction à loutes les températures. 


1. Ce Mémoire est le premier d’une série destinée 


à développer les résultats contenus dans des Notes - 


préliminaires [1]. Les calculs ont été refaits sur de 
meilleures basés, et l'accord entre la théorie et 
l'expérience a pu être poussé dans de nombreux 
cas jusqu’à la perfection, comme on va le voir. 

Le but général du travail est de trouver une 
expression mathématique permettant de calculer à 
toute température la ‘pression d’un gaz en fonction 
de son volume. C’est le problème que l’on cherche 
d'ordinaire à résoudre par l'équation d'état : la 
solution que j’apporte et qui est de beaucoup préfé- 
rable, est de nature différente. Elle repose sur la 
notion d’association ou condensation progressive 
que je soutiens depuis 1924 [2]. 

Cette première partie est consacrée à l'oxygène. 
Les chiffres expérimentaux sont ceux de Kamerlingh 
Onnes [3]. 

Je donnerai immédiatement les résultats numé- 
riques : in exlenso, pour la température de — 1160,01 
choisie à titre d'exemple; en abrégé pour les autres. 


TaBLeau Î, 


Comparaison du calcul et de l'expérience à — 116°. 


44 
BA 
Pression gs. LA RS 
{atm. }. observé, calculé, Différence. 
22,30-.:-0,483D 0,4839 +0 ,0004 
27585 ; 10, 4500 0,4560 0 ,0000 
92:02 %0:4207 0,4296 - —0,0001 
DUT 0 4OTD2 0,4007 —0,0005 
49, 90e400,93541 0,3537 —0,0004 
50,47 0,2826 0,2851 +0 ,0000 
52:44 0,2477 0,2476 —0,0001 


La comparaison est aussi satisfaisante que possible. 
Les écarts ne montrent pas de marche systématique 


et leur valeur moyenne est 0,0003. Ce chiffre est 
d’ailleurs un maximum, car la méthode de calcul 
entraîne des erreurs pouvant atteindre 0,0002. 

Pour les réduire il faudrait se servir de Tables 
de logarithmes à 6 décimales au lieu de 5; les calculs 
étant déjà très longs, j'ai pensé que cette compli- 
cation était inutile, puisque les écarts ne sont pas 
supérieurs aux erreurs expérimentales. 


Pression (atm.) 


20 30 40 50 
Fig. 1, — Oxygène à — 116,99. 

2 pv 

Valeurs de = 


RAT) écarts entre la théorie et l'expérience. 


A toutes les températures supérieures la concor- 
dance est du même ordre : 


Mesures, Écart moyen, 

A RD OS 6 0,0002 
Cet c 9 0,000 

= HO OD RE 7 0,000 
— 80,03.+.:. 9 0,0003 
—102,46..:,, 5 0,000 
— 109,97: :::: 10 0,0007 
— 113,97... 9 0,0003 


Trois écarts seulement sont supérieurs à 0,001, 
sur un total de 58. Leur variation n’est jamais systé- 
matique. 
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Mais à la température la plus basse (—1160,99, 
très près du point critique —1180,8) des écarts 
plus forts apparaissent aux fortes pressions. Ils sont 
donnés par la figure 1. Jusqu'à 5r atm ils restent 
comparables aux erreurs d'expérience, mais au- 
dessus ils augmentent. L'’irrégularité des chiffres 
expérimentaux ne permet pas de donner de chiffres 
précis, mais il semble que l’écart maximum sous la 
pression la plus élevée soit 0,00 environ. Celui 
que donne la formule de van der Waals est 0,232, 
soit 46 fois plus grand. Nous allons voir que la 
méthode par laquelle nous obtenons ce résultat 
beaucoup plus favorable est très simple. 


2. Bases du calcul. —— Comme il s’agit ici d’un 
mode de calcul nouveau, je serai obligé de l’exposer 
en détail. 

Les principes ont été exposés antérieurement [4]. 
Nous négligeons les actions à distance de la théorie 
de van der Waals. Nous admettons que dans les 
gaz réels une certaine fraction des molécules sont 
associées, fraction d'autant plus grande que la 
pression est plus forte ou la température plus basse. 
Nous admettons de plus que l’association est progres- 
sive, c’est-à-dire que les groupes de 2 molécules 
formés par la première association peuvent former 
des groupes plus nombreux de 3, 4, 5, ... molécules, 
sans limite supérieure. Maïs si nous voulions tenir 
compte de tous ces groupes nous serions conduits 
à des calculs inextricables. Pour simplifier, nous 
admettrons que l'association se fait en progression 
géométrique, c’est-à-dire que nous n'avons à Consi- 
dérer que les groupes de 1, 2, 4, 8, 16, ..., et en 
général 2* molécules. Pour simplifier l’écriture nous 


2 » Là * 
désignerons par M, les groupes de 2* molécules 


simples, et par C, leur concentration dans les unités 
qui seront définies tout à l’heure. 

Nous admettrons que, entre GC, et C,:, existe 
une relation de la même forme que celle qui exprime 
la loi d'action de masses 


CHERE CES (1) 

Cette relation permet de calculer de proche en 
proche toutes les valeurs des C connaissant celles 
des K,. 

Nous ne savons pas a priori comment ces valeurs 
dépendent de la valeur de n. J'ai déjà indiqué [4] 
les raisons de croire à une variation régulière de K, 
en fonction de n. L'hypothèse la plus simple est 
qu'il y a entre deux coefficients successifs la relation 


K}: == pr (2) 


le nombre p étant pour chaque gaz une constante 
indépendante de n et indépendante aussi de la 
température. Toutes ces hypothèses ont pour but 
d'introduire le plus petit nombre possible de para- 
mètres. 
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3. Variation du coefficient K, avec le volume. 


— Le coefficient X, est d’après nos définitions celui 
qui correspond à la première associalion (deux molé- 
cules simples se réunissant en une double). Toutes 
les concentrations €; se déduisant comme nous 
l'avons vu des valeurs des K,, et celles-ci se déduisant 
de K, par l'intermédiaire de la constante p que 
nous venons de définir, tout se ramène finalement 
à la valeur de K,. 

La formule (1) a la même forme que celle de la 
loi d’action de masses. Dans les applications que 
l’on fait de cette loi on admet toujours que la valeur 
de la constante K est indépendante du volume, 
c’est-à-dire qu’elle conserve la même valeur si le 
volume du système est changé. C’est alors une 
véritable constante, pouvant entrer comme telle 
dans les calculs thermodynamiques. C’est l’hypo- 
thèse que j'avais faite au début, ayant constaté 
qu'elle était suffisante pour de faibles pressions. 
Mais il est-impossible de la conserver pour des 
pressions élevées et des températures basses, car 
elle ne peut rendre compte de la liquéfaction. 

Il est nécessaire alors [si du moins l’on veut 


conserver la formule (1)] d'admettre que le coeffi- 


cient K, est d'autant plus ‘petit que le volume du 
gaz est plus petit; c’est-à-dire que si nous prenons 


une quantité donnée de gaz, une molécule-gramme 


par exemple, le coeflicient K, dépendra du volume 
occupé par cette masse de gaz ou de sa densité. 
Ne sachant pas a priori de quelle manière il en 


dépend, nous ferons une fois de plus l’hypothèse 


la plus simple, et qui rend les calculs les plus faciles. 
Nous admettrons que pour une masse donnée 
de gaz, le coefficient K, est proportionnel à une 


puissance m du volume occupé par cette masse. 


Si nous désignons par K,, le coeflicient qui fixe 
l'équilibre, dans le volume », des assemblages de 2” 
et de ceux de 2**! molécules simples, nous aurons 
ainsi | 
Kyv = Ko VIRE 
soit 
logK,,= log Ko1+ logo + nlogp. 


(3) 


Je répète pour éviter toute confusion que K5: 
et » dépendent de la température tandis que p est 
indépendant de la température et du volume, 

Les hypothèses faites peuvent sembler très 
étranges et difliciles à justifier. Je ne me suis pas 
arrêté à cette objection. Il est certain que tout 
ce qui précède peut paraître compliqué et difficile 
à assimiler; tout point de vue nouveau apporte 
un inconvénient de ce genre. D'ailleurs tout l'essentiel 
de la méthode tient dans la formule (3), qui est 
très simple et d'application facile, bien qu elle 
conduise à des calculs fort longs. 


4. Covolume et gaz fictif. — Le calcul numé- 
rique montre qu’il est nécessaire de conserver la 
notion de covolume; cette question sera discutée 


quand nous examinerons les propriétés de l'azote. 
Nous emploierons cette notion sous une forme 
différente de la forme ordinaire en faisant intervenir 
ce que nous appellerons le gaz fictif. C’est un gaz 
auquel nous supposerons toutes les propriétés du 
gaz réel, mais dont le covolume est nul. Les coefi- 
cients À, m, p du gaz fictif sont identiques à ceux 
du gaz réel, sous des volumes correspondants. Par 
exemple, pour un volume du gaz réel égal à 1, si 
le covolume est 0,002, les coefficients K, m, p sont 
identiques à ceux du gaz fictif sous le volume 0,998. 
La pression est aussi la même; tout se passe comme 
si le volume des molécules s’évanouissait sans que 
les propriétés élastiques changent. Tous les calculs 
se font sur le gaz fictif : à la fin on ajoute le covolume 
pour passer au gaz réel.- 


Nous continuerons à représenter le covolume par 
la lettre b et le considérerons comme indépendant 
de la pression et de la température. 


5. Exemple de calcul. — En l’absence de tout- 


guide théorique les valeurs de X,, m, p et b ne peuvent 
s’obtenir que par approximations successives. Cette 
nécessité conduit à de longs calculs et a retardé de 
près d’une dizaine d’années la parution du présent 
Mémoire. Il est nécessaire de dresser des Tables 
numériques dont on trouvera un abrégé plus loin. 
Mais une fois que ces Tables sont calculées, le travail 
devient très simple. 


Soit l’oxygène à la température de — 1160,o1 
et sous la pression de 43,947 atm. Nous prenons 
les unités Amagat; le volume 1 est celui du gaz 
sous la pression de 1 atm à o°. Alors l’expérience 


donne pour pv le nombre 0,35408 d’où l’on déduit 
pour v la valeur 0,008057. Pour le calcul théorique 


nous partirons de cette valeur de » pour en déduire 
celle de p. 


Après quelques tentatives, il a été reconnu que 
la valeur la plus convenable du covolume était 0,0021. 
Le volume du gaz fictif défini plus haut est donc 


0,008097 — 0,0021 = 0,005997. 


Il a été reconnu aussi que les valeurs les plus 
favorables des constantes étaient 


log Ko 1 = 2,63524, m — 0,39226, : p—0,8072. 
Pour avoir une unité de concentration nous admet- 

trons que la concentration des molécules simples 

dans le volume 1 serait égale à l’unité, si tous les 


assemblages étaient détruits. 


La valeur de la constante de première associa- 
tion K,, sous le volume du gaz fictif est donnée 
d’après la iormule (3) par la relation 


log K5,= 2,63524 + 0,39226 log 0,005957 = 1,76247. 
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Nous calculons alors de proche en proche 
l0g AK = 196247, Co= 36,85; 
log A1 = 1,669457, (ES EEE 
log Xs = 71,55643, CS: 70; vie 
logkK:= 1,48347, C3; 09; 
log Ki = 1,39047, Cr 049: 
Cy—= 0,008. 
La valeur de C, doit être telle que la concen- 


: 2 : : ; SRE 
tration des. molécules simples soit égale à. en 


Supposant tous les assemblages détruits, c’'est- 
à-dire que lon ait 
Ai . n Q fr A I a. > 
Co+2Ci+ 4Cr+8C3+16C,+..—= = — 167,87. 
Ù 


Pour le calcul de Er il faut tenir compte de ce 


que chaque assemblage, quelle que soit sa complexité, 
forme une unité cinétique unique. Le nombre total 
des unités cinétiques est donc 


C = Ci += CS + C> ET are ES 76,25 


76,25 


PV est 


RT — 0,4542, d’où 


et la valeur de —— 
167,97 


978,19 —=110,01 


ET = {3,88 atm, 
273,15 


P 


au lieu du chiffre expérimental 43,05. 


6. Courbe de condensation. — Ce mode de 
calcul est complètement différent des modes usuels 
qui reposent sur l'emploi d’une équation d'état. 


: Pt 
Celle-ci pe 
données par le quotient de deux séries dont chaque 
terme correspond à un certain degré de conden- 
sation. La notion de pression interne ne correspond 
plus à rien. Nous verrons plus tard qu'il n’y a plus 
en général d'états strictement correspondants, et que 
l'importance attachée à la température critique est 
exagérée. Tout ceci peut paraître assez choquant, 
mais là encore je ne me suis pas arrêté à ces 
objections. 

Dans le système van der Waals, on est naturel- 
lement conduit, si l’on veut tracer une courbe de 
compressibilité, à prendre comme abscisses les 
volumes et comme ordonnées les pressions. Nous 
sommes conduits ici à employer un tout autre 
système qui a permis de découvrir une propriété 
intéressante des gaz, restée inaperçue jusqu'ici. 
Nous prendrons comme abscisse le logarithme du 
volume du gaz fictif (volume que nous désignerons 
désormais par v/) et comme ordonnée la valeur 


est supprimée. Les - valeurs de sont 


de Ee Ce choix est dicté par la nature des calculs. 


Nous obtiendrons ainsi une courbe en forme de S 
avec deux asymptotes, l’une à l’ordonnée 1 pour 
les très grands volumes, l’autre à l’ordonnée zéro 
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pour les très petits. Quelques-unes de ces DES 
sont reproduites sur la figure 2. Elles ont toutes 


US 10 15 =-0 


— Courbes de condensation expérimentales 
à diverses températures. 


Fig,.2. 


la même forme générale et leur pente est d'autant 
plus forte que la température est plus basse. Pour 
un gaz parfait elles se réduisent à une droite hori- 
zontale d'ordonnée égale à l’unité. 


7. Relation nouvelle. — Ces courbes se déduisent 
les unes des autres par une relation remarquable 
qui apporte une grande simplification dans l'étude 
de la compressibilité des gaz. Elle se transforment 
l’une dans l’autre par une translation horizontale 
et une dilatation horizontale. À deux températures 7: 
et 7, la même ordonnée y correspond à deux valeurs x, 
et x de l’abscisse entre lesquelles on a la relation 


La —= A%: RE B, 


À et B étant deux constantes indépendantes de la 
pression. Si donc une courbe est complètement 
tracée pour une certaine température 7, il sufit 
de déterminer deux points à toute autre tempé- 
rature 7, pour calculer les constantes À et B et 
en déduire le tracé complet de la courbe T,. En ce 
sens on peut dire qu'il y a pour un gaz une courbe 
de condensation unique, indépendante de la tempé- 
rature. Mais cette courbe diffère d’un gaz à l’autre. 

L’exactitude de cette loi est parfaite. Dans tout 
le domaine étudié par Kamerlingh Onnes (+ 20° 
à —1170et 20 à 60 atm), les écarts avec l’expérience 
sont inférieurs aux erreurs expérimentales. 

Nous sommes ainsi très loin des formules du genre 
van der Waals. Cette loi remarquable ne pouvait 
en être déduite. 

Il en résulte que la compressibilité de l’oxygène 
peut être représentée par une courbe (ou une Table) 
unique. J’ai calculé cette Table pour toutes les 
valeurs de (logp;) comprises entre — 2 et —#;1 
par intervalles de o,o1. Dans cette Table l’unité 
de volume est arbitraire. Il suffit de connaître pour 


chaque température la dilatation horizontale et la 
translation définies plus haut et qui sont repré- 


“JOURNAL DE PHYSIQUE 


sentées par les nombres suivants : 


- Dilatation. lirinstation ? 


DORE. A) 0,979 
DA EE PL ATTI I,0018 
40500 1,1020 1,0643 
200, 09%: e5E;0070 1,0800 $ 
—102,46 1,025 1,0410 FA 
— 109,97 1,0050 1,0140 
—113,97 0,9881 0,9827 
— 116,01... -0,9840 0,9803 
LTD JON NTO0S 002 0.097094 


Je donne un extrait de la T'able de condensation. 


Table de condensation de l’oxygène. 


(Valeurs à de Le en fonction de log.) : 

log ep ei loge r, Ed oc ÿ 4 AU 

RT Fe RT PRET RE 
2,0 0,004 3:3;..10,7734 4,6.::"0,3629 
3.9 0,940 F2 OST 4.5 0,3030 
3,8 0,9235 SE 0,6751 4,4... .0:2486 
3,7 0,9028 3,012: %6,0177 4,8. 1032003 
3,6 0,8781 4,9: 7 0,5561 4,2 0,1592 
3,5... 0,8485 ASE. 0 4917 4,1... 0,1248 
3,4... :0,8138 4,7... -0,4265 4,0... 0,0968 


8. Emploi de la Table. — Soit la température 
de — 800,03 avec le volume réel 0,01858/4, Le Covo- 


lume étant o,oo21, le volume », du gaz fictif est. 


vf = 0,018584 — 0,0021 = 0,016484 


«dont le logarithme est 2,21706 ou —1,78294. De 
ce chiffre il faut retrancher le facteur de translation 
donné plus haut (1,0800), ce qui donne — 2,86294. 
En divisant par le facteur de dilatation, 1,06%0, 
nous obtenons — 2,6832 qui est la valeur de (log v;) 
à prendre dans la Table. Nous voyons en face 0,7806 


pour la valeur de 5 ce qui donne pour ppy, 0,524. 


En divisant par la valeur 0,016484 de »; nous 
obtenons finalement une DER de 33,51 atm. Le 
chiffre expérimental est 33,475. Tous les autres 
chiffres de Kamerlingh Onnes ont été traités de la 


même manière et les résultats ont été donnés en. 


tête de ce Mémoire. Par le moyen de la Table les 
calculs sont extrêmement simples, beaucoup plus 
que par la formule de van der Waals pour laquelle 
il est impossible de se servir d’une Table unique. 
Il en faudrait une pour chaque température, La 
simplification est done considérable. 


1 


9. Variation des facteurs avec la tempé- 
rature. — Seuls dépendent de la température les 
facteurs K, (constante de première association) 


et m (coefficient indiquant l'influence de la tempé- 


rature sur la condensation). 


x BE 


sen), 21 Les 


MON TE 


_ Pour des raisons qui seront données plus tard, 
il est préférable d'établir la variation des coeffi- 
cients X, non pas sur les valeurs qui correspondent 
à la densité sous la pression atmosphérique, mais 
sur celles qui correspondent à une densité 100 fois 
plus grande (ou à un volume o,o1). Soient K” ces 
coefficients; nous trouvons qu’ils se déduisent d’une 
loi très simple 


log À’ = 2,9082 — one 


(4) 


La correspondance des nombres donnés par cette 
formule avec les nombres expérimentaux est remar- 
quablement exacte : 


L Observé. Calculé. Différence, 
SOU 2,341 2,341 0 
Dar 2,209 2,299 0 
== AO +... 2,191 2,195 +0,004 
O0 Priute 2,045 2,047 +0 ,002 
—102 1,930 1,934 0,002 
TEO ser 1,Sgt 1,809 — 0,00? 
—113,94 1,860 1,564 — 0,001 
-LI0 OT 1,001 1,800 —0,001 
110,99: .. 1, 844 1,844 ) 


Si l’on interprète à la manière ordinaire la for- 
mule (4), elle signifie que la chaleur d'association 
est constante et égale à 760 calories. La constance 
serait une vérification de la théorie, car on pouvait 
l’'admettre a priori. Mais il faudrait examiner si 
cette interprétation est encore possible ici, le coefi- 
cient X n'étant pas une véritable constante puisqu'il 
dépend de la densité. Ce point ne peut être discuté 
ici. En tout cas l'exactitude de la formule permet 
-l’extrapolation vers les basses températures, 


Température °C 
-100 -50 0 


Fig. 3. — Oxygène. 
Variation avec la température du coefficient m. 


Pour ce qui concerne le coefficient m les difficultés 
sont beaucoup plus grandes. Ses variations au 
voisinage de la température critique sont très 
rapides (fig. 3). Elles sont assez bien HOTSRES 
par une formule du type 


log m = À — 7? 
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. mais à peu près aussi bien par une formule en T%. 


Pour choisir entre ces deux formules, et beaucoup 
d’autres, il faudrait un guide théorique. Elles 
donnent par extrapolation des chiffres très diffé- 
rents; par exemple, à 1300 absolus, 0,770 (T%) 
et 0,6646 (1%). Il en résulte que l’extrapolation 
vers les basses températures est fort incertaine, 


dès qu’on s’écarte de plus de quelques degrés de 


la température critique. 


10. Liquéfaction. — La théorie conduit à se 
représenter de la manière suivante le phénomène 
de la liquéfaction. 

Quand le nr diminue, l’association progresse 


et la valeur de 2 diminue comme le montrent les 


Re ï k 
courbes de la figure 2. Dans le calcul de ces courbes 
là pression n'intervient en rien; elle s’en déduit 
seulement par l'identité 


_ Pvr ES 


DER DrE a 


dont tous les éléments sont connus. Si l’on fait le 
calcul on s'aperçoit que lorsque le volume diminue 
la pression, ainsi calculée, n’augmente pas indéfi- 
niment : elle passe par un maximum pour décroître 
ensuite, si la température est assez basse. 

C’est la position de ce maximum qui est le point 
de liquéfaction. En effet, au delà le système est 
instable toute nouvelle diminution de volume 
amenant une diminution de pression, la conden- 


sation doit devenir totale sous pression constante. 


La position du maximum de pression sur la 
courbe de condensation est d’ailleurs facile à déter- 
miner en ce point la pression passe par un 
maximum (dp — 0). On a donc, d’après la relation (5), 


RT É pure RT 
dd RAT di de —0, 
ce qui se réduit à 
© te œ 
d log sr = dlogvr. 


L’abscisse æ des Ar étant (log v;), tandis que 


l’ordonnée y est 2 Es cette relation peut s’écrire 


d'logy = dx. (6) 

Il suffit donc de construire la courbe (log y, +) 
et de prendre le point auquel la tangente à cette 
courbe est inclinée à 459. 

De cette manière la courbe de condensation est 
en même temps la courbe de liquéfaction. C’est là 
un avantage très considérable sur les formules du 
{type van “der Waals qui ne permettent aucune 
prévision de ce genre. 

Cet avantage doit bien entendu être contrôlé 
par l'expérience. Les mesures de Kamerlingh Onnes 
ne permettent pas de le faire directement, étant 
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toutes relatives à des températures supérieures au 
point critique. Mais nous pouvons essayer de les 
extrapoler vers les basses températures et nous 
obtenons alors : 


Pression (atm) 


Température mn 


(°K). calculée. observée. 
TOOL 44 44 
TO a 33 29 
LRO 25 17 


La concordance est très bonne à 1509. Au delà 
les différences augmentent très vite, mais elles sont 
attribuables à l'incertitude de l’extrapolation. Les 
calculs ont été faits avec la formule d’extrapolation 
en 7% et la valeur calculée de m est 0,770, plus du 
double du dernier chiffre expérimental. La vérifi- 
cation de la théorie est aussi bonne qu’on pouvait 


l'espérer. Le principe apparaît comme parfai- 
4 tement sain. 
11. Résumé. — Les idées qui nous ont servi 


de guide au long de cet exposé étant très différentes 
de celles qui ont cours actuellement je résume ici 
la marche suivie et les résultats obtenus. 


10 La théorie que je propose vise à remplacer 
complètement toutes celles qui sont fondées sur 
le principe de van der Waals. Elle repose sur l’idée 
de la condensation progressive. 


29 Nous désignons par C; la concentration des 
concentration étant celle des molécules simples 


dans le gaz à l’état normal, supposé complètement 
dissocié sous le même volume. 


je 30 Toute la théorie tient dans les formules 
suivantes : : 
CE ES KncCrr, 
loge. = logÆXo1+ mlogv+nlogp, 
: 166,2 
log Ko,-0,01 = 2,9082 ° 


ne 


qui permettent de calculer de proche en proche les 

concentrations C et de définir ainsi complètement 

l'état de condensation du gaz. Le paramètre p est 
= une constante caractéristique du gaz. Le paramètre m 
dépend de la température suivant une loi qui n’a 
pu être établie. Finalement l’état du gaz à toute 
lempéralure dépend de quatre constantes fixes et 
d’un paramètre variable, 


[1] C. R. Acad. Se., 1948, 226, 1113; 1948, 226, 2034; 1948, 
227, 1124. 
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assemblages de 2° molécules simples, l'unité de. 


[2] J. Physique Rad., 1924, 5, 3313: 1927,8, 3363 1940, 1, 293. 


4° Par le moyen de ces relations la compres- 
sibilité de l’oxygène est représentée avec la plus 
complète exactitude entre “+ 209 et —- 1179, sauf- 
dans un domaine très petit au voisinage immédiat 
du point critique. L'écart ne dépasse pas 2 pour 100 
de celui que donne la formule de van der Waals. 


59 La théorie a permis de reconnaître une nouvelle 
propriété des gaz qui simplifie beaucoup l'étude 
de la compressibilité, puisqu'on peut la déduire 
d’une Table unique valable pour toutes les tempé- 
ratures. ; 


69 La courbe correspondant à cette Table est la 
courbe de condensation. Elle porte un point singulier 
qui est celui auquel le gaz se liquéfie. Aïnsi une 
Table unique reproduit à la fois la compressibilité 
et la liquéfaction. Si la variation du paramètre m 
avec la température était connue, il suffirait d’une 
seule courbe de compressibilité à une seule tempé- 
rature pour connaître la compressibilité et la liqué- 
faction à toute température. 


12. Pour alléger ce Mémoire, j'ai dû laisser de 
côté un certain nombre de points relativement 
accessoires. Ils seront discutés dans une deuxième 
Partie relative à l’azote. Chaque gaz doit être traité 
isolément, car chacun a sa propre loi de compres- 
sibilité. Les courbes de condensation qui conviennent 
à l'oxygène ne conviennent pas à l'azote, et l'état 
de ces deux gaz au point critique n’est pas le même. 
I n’y a pas entre eux d'états strictement corres- 
pondants. 

J’indiquerai immédiatement que la loi de variation 
des constantes d'association X avec le volume ou la 
densité, qui a été adoptée dans ce Mémoire, n’est 
pas définilive. Maïs pour établir la loi véritable, de 
très longs calculs seront encore nécessaires. De plus, 
la nouvelle loi devra sur un intervalle étendu se 
confondre avec l’ancienne, de telle sorte qu'il ny 
aura presque rien de changé aux résultats qui ont 
été donnés ici. 

Enfin, on pourra dire que les molécules associées, 
sur l'existence desquels la théorie repose, n'existent 
pas dans la réalité. Je pense qu’elles existent réel- 
lement, mais la question n’a pas grande impor- 
tance, Même s’il ne s’agit què d’une nouvelle 
méthode de calcul, ses avantages sont tels qu’elle 
devra être substituée aux anciennes, fondées sur 
des idées dérivées de celles de van der Waals. 


{ } 
Manuscrit reçu le 13 février 1950. 
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REMARQUES SUR LE CALCUL DE LA PRESSION THERMODYNAMIQUE 
! EN STATISTIQUE QUANTIQUE 


Par J. HOROWITZ. 
Fort de Châtillon, Fontenay-aux-Roses. 


Sommaire, — Différents calculs quantiques de la pression thermodynamique ont été présentés 
récemment pour démontrer que celle-ci est ou n’est pas égale à la pression cinétique. La question a été 
reprise ici. Nous concluons à l'égalité des deux pressions et analysons l’origine de certaines inexactitudes, 


Considérons N particules identiques dans un 
volume V. Nous supposons que ces particules n’ont 
pas de degrés de liberté internes. Le potentiel d’inter- 
action est o(r;;) et il n’est pas nécessaire pour la 
suite de préciser à quelle statistique ces particules 
obéissent. L'équilibre étant réalisé à la tempé- 
rature 7°, la pression thermodynamique p est donnée 
par la relation 

ADI RATIO Fe 
IST D TALRS co 


F est l’énergie libre, Z la somme d’état définie par 
En 
AS et, (2) 


la sommation étant étendue sur tous les états 
propres de l'énergie. Les dérivées partielles sont 
prises à température constante. 

Si les particules étaient libres de se mouvoir dans 
tout l’espace nous pourrions, sans précisions supplé- 
mentaires, écrire l’hamiltonien ä€ sous la forme 


HE — RE ser > Sr 9 
x=Y omPi Ts Prix). (3) 
n jk 


Nous devons cependant décrire le rôle des parois 
limitant le volume V. Pour cela, nous pourrions 


introduire un mur de potentiel de grande hauteur 


mais qui aurait l’inconvénient de laisser déborder 
les fonctions d’onde en dehors du volume V et de 
donner pour de grandes énergies un spectre continu. 
Il vaut mieux limiter les fonctions d’onde au 
volume V : il devient alors indispensable d'introduire 
des conditions aux limites. Nous imposons aux 
fonctions d’onde d’être nulles sur la surface limitant 
le volume .V, cette condition se déduisant assez 
naturellement de la considération d’un mur de 
potentiel dont la hauteur tend vers l'infini. 
D’autres conditions sont mathématiquement pos- 
sibles. D'ailleurs le raisonnement qui va suivre et, 


par suite, la. formule finale, sont dans une large 
mesure indépendants de ce choix. 

Les conditions aux limites assurent l’hermiticité 
de l’hamiltonien ëäC et précisent son domaine de 
définition, ses fonctions et valeurs propres. À deux 
systèmes de conditions aux limites différents corres- 
pondent, en général, des opérateurs totalement 
différents (1). 

Modifions maintenant légèrement le volume V. 
Pour simplifier, nous allons lui faire subir une homo- 
thétie et nous ramener au volume primitif par le 
changement de variables 

= 
I1+E 


, 


Nous avons ainsi à considérer dans le même 
volume V, avec les mêmes conditions aux limites 
(nullité sur la surface), les hamiltoniens dépendant 
du paramètre € : 


I I 1 KT 
JC(E) — = Fe o[r;Kx(1+e)l]. (4 
1 Grec) 2m > , Del J e)] ) 
È 


jk 


et les sommes d’état 
77, (€) 


LS RENE (5 


Pour obtenir l'expression de la pression p nous 
devons calculer 


2Z(0). 1 d7Z(0) 
OV Le] dE 
E, (0) 
= peE EP EE (2). 200) 
3VKT dE ) 
n 
V0E,(0) 


Te s’obtient par la méthode de perturbation 


() Voir par exemple J. von NEuMAN, Mathemalische 
Grundlagen der Quantenmechanik, p. 79. 
d7(o) our 2Z(°) | < 
p de |:—0 


() Nous écrirons 
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habituelle 
dE,(0) 928€ (0) 
——— — n ni 
dE dE 
; HNLE IN. (rx) de 
see | ri DL or ;x n). (1) 
i jE 


où |n>> représente le nit"° état propre de Phamil- 
tonien #(o). De (1), (6) et (7) on déduit 


pe I a 
Hart 


de(T;k) = 
l';k D R ne F7 
Or jk 


SO GE I de(r;e) 
5 »m LP 


Î 


(8) 
i jh 

où la barre indique la moyenne sur l’ensemble cano- 

nique. La relation (8) exprime l'égalité des pressions 

thermodynamique et cinétique. Le mode de raison- 

nement ci-dessus a déjà été publié (°). 

Born et Green conduisent les calculs autrement (4 
et concluent à l'inégalité des pressions thermody- 
namique et cinétique, la différence devant surtout 
être sensible à basse température. De Boer (°) 
établit la relation (8), maïs son calcul renferme des 
imprécisions analogues à celles qui conduisent au 
résultat inexact de Green. Examinons maintenant 
ce qui est à critiquer dans ces démonstrations. 

Nous avons déjà insisté sur le point suivant 
c'est seulement après introduction des conditions 
aux limites que les hamiltoniens 4€%(:) (relatifs au 
volume V) se trouvent définis à partir des opéra- 
teurs ä€, (es) (applicables à tout l’espace). La multi- 
plicité des fonctions d'onde pour lesquelles les dy (e) 
sont définis est plus restreinte que le domaine de 
définition des ä€, (e). Il en résulte que certaines 
égalités entre opérateurs cessent d’être vraies et que 
certains opérateurs construits au moyen des 9€,(e) 
cessent d’avoir un sens quand ü€,(e) est remplacé 
par 2x (e). De Boer et Green (f) écrivent - 


ae (=) 


Z{e) = tracee:.T" 


ie (0) 
2Z7(0) I Care 7e 
TANT AT RSE 
€ (0) 
I ere 
=RESrAtRACe AE KT 
3} 5 ( 
äc (0)  ( ZT) I J8(o)|. 
Dr ÿ RE 
LS AT l AT de SF en €(o) E 


dE 


MT (ATY 


ae DEN 
< | (0). | A0). CE +... | CRT 
ù dE \ 


(5) J. Yvox, Colloque de Mécanique statistique de Florence. 
mai 1949. L'auteur utilise les conditions cycliques de Born 
comme conditions aux limites. 


() H.S, GREEN, Physica, 1949, p. 88». 
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sans se préoccuper de savoir si i les opérations précé- 
dentes et ce développement ont un sens. 

De Boer affirme que les commutateurs ne donnent 
pas de contributions dans le calcul de la trace. 


Green, par contre, calcule ces commutateurs en 


À 


utilisant en fait les propriétés des opérateurs 4€, (€) 
et obtient ainsi un résultat différent. 

Nous allons illustrer ces points à l'aide d’un 
exemple simple. Considérons une seule variable x 
dans l'intervalle (0,1). Pour rendre RES et 


définir réellement l'opérateur A (:) — il; ee ex dans 


l'intervalle considéré, nous imposons aux fonctions - 


d'onde la condition aux limites (1) — f(o). 

Les valeurs propres sont 

: £ 
En(e)=2nr+ >: AVEC HN 0,. ÉT, ho 2 


les fonctions propres 
AnCE)=="C aie (E)x—© | 


En conséquence 


9E,(0) 1 
CE 2? 


ce qui résulte aussi facilement du calcul de pertur- 


bation puisque 
dA(0) 


f fito)r fo) a =. CCE 


D'autre part il vient 


)E> ; 
ol, PAC y y 
Si pour calculer ie nous voulions suivre une 


méthode analogue à celle de De Boer ou à celle de 
Green, nous devrions considérer les opérateurs et 
écrire 


A0) 4e To LOT 
dE 
=» At)+ [At 0), | | 
ee ne féto 404 (0 )fto)de 
ROIS ao 


ZE, (0 + fé co Lo) 22] 7, Co). 


Dans un cas analogue de Boer affirme que le 
deuxième terme, qui implique le CORNE UE est 


1949, p. 843. 


(5) Nous suivons ici le développement de de Boer. 
de Green lui est pratiquement équivalent, 


(5) DE BoERr, Physica, 


Celui 


+ af 
N° à. 
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nul, ce qui conduit au résultat correct. Si, en suivant 
l’esprit du calcul de Green, nous écrivons 


[co “1 ë Ê p) 1= à 


dx" 
dE? 2 
nous obtenons pour / Ja valeur inexacte E, (0) +i 
db . 0) 


> 
En réalité, l'expression À (0) 


fn(o) écrite ci- 


dessus n’a pas de sens : la ae 


JA( 
ro) = A Mae 
ne satisfait pas à la condition de périodicité et même 


par prolongement il n’est pas possible de définir 
l'opérateur A(o) — le ; Pour cette fonction. Pour s’en 


rendre compte, 
pement 


il suit de considérer le dévelop- 


Le Ênk 


a —— er?rikx, 
27 KE — | 


SUR LE MÉSON + 243 


où 
I si 


{ 
to si 


d 
- appliqué au second membre-donne 


n = k, 


dE — 
G n £ k. 


L'opérateur ir 


Cette série ne converge pas en moyenne ce qui 
montre l'impossibilité du prolongement. 

Cet exemple montre quelles erreurs peuvent s’intro- 
duire dans les calculs s’il n’est pas tenu compte de 
l’existence de conditions aux limites imposées aux 
fonctions d’onde et des modifications que celles-ci 
apportent aux propriétés des opérateurs. 


Je tiens à remercier le Professeur J. Yvon de 
m'avoir conseillé d'entreprendre ce travail. 


Manuscrit reçu le 13 janvier 1950. 
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CONSIDÉRATIONS SUR LE MÉSON - 
Par Louis LEPRINCE-RINGUET et CHARLES PEYROU. 
Laboratoire de l’École Polytechnique, Paris. 
Sofnmaire. — Discussion des phénomènes observés favorables à l'existence d'un méson 7 et des 


expériences systématiques entreprises pour trouver ce méson et qui ont donné des résultats négatifs, 


Depuis plusieurs années, quelques photographies 
en faveur de l’existence d’un méson chargé (méson 7) 
de masse intermédiaire entre celle du méson 7 
et celle du proton ont été présentées par divers 
auteurs [1 à 6]. 

La première [1] de ces photographies présente un 
électron de collision obtenu dans une chambre de 
Wilson avec champ magnétique; la masse de la 


particule incidente, mesurée en supposant la collision 


élastique, est de l’ordre de r 000 m, (m,, masse de 
l’électron). 

La seconde photographie [2] obtenue dans une 
chambre de Wilson avec champ magnétique s’inter- 
prète comme la désintégration en vol d’un méson 
très lourd chargé donnant naissance à une particule 
neutre et à un méson chargé. 

La troisième indication [3] présente deux étoiles 
observées dans une plaque photographique et reliées 
par un méson 7 négatif : la première étoile donne 


naissance au méson 7 et à d’autres branches habituel- 
lement observées dans les étoiles, et paraît due à une 
particule se dirigeant vers son centre et s’y arrêtant : 
l'aspect général de la trajectoire correspondante 
ainsi que le bilan d’énergie conduisent à attribuer 
à cette particule une masse supérieure à 800 m,. 

La quatrième [4] indication présente la trace 
d'une particule se terminant dans l’émulsion d’une 
plaque photographique après avoir parcouru une 
longueur visible de plus de 2500 et donnant 
naissance à une étoile de trois trajectoires dont l’une 
est due à un méson r négatif lent facile à reconnaître, 
dont une autre paraît due aussi, d’après les données 
expérimentales à un méson r et dont la troisième ne 
peut être identifiée directement; l'étude du bilan 
d'énergie et d’impulsion de ces trois trajectoires, qui 
sont d’ailleurs coplanaires, est en faveur de la tripar- 
tition d’un méson + en trois mésons. 

La cinquième [5] est obtenue à la chambre de 
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Wilson avec champ magnétique dans des conditions 
analogues à celles de la seconde. Elle indique l’exis- 
tence, au sein d’une gerbe provenant du plomb 
situé au milieu de la chambre, d’une particule 
chargée négativement et dont la masse est nettement 
supérieure à celle d’un méson 7. Les indications 
d’ionisation et de courbure sont favorables à une 
particule négative dont la masse est de l’ordre 
de r 000 my. 

Ces photographies ont été discutées dans les 
publications correspondantes et nous ne ferons pas 
à nouveau de discussion générale, mais nous revien- 
drons sur certains points qui nous semblent inté- 
ressants. Depuis la publication de ces photographies, 
de nombreux travaux ont été effectués, notamment 
avec des plaques photographiques, et la quantité 
de documents actuellement dépouillés par les 
physiciens des divers laboratoires, s’est accrue d’une 
façon considérable. 


Expériences systématiques. — En particulier, 
plusieurs séries d'expériences systématiques ont été 
effectuées : les unes par le groupe de Bristol, une 
autre par le groupe de l'École Polytechnique et une 
par le groupe de Berkeley. 28 

Considérons, tout d’abord, cette dernière conduite 
par Brode, Fretter, Retallack [6]. Elle est essentiel- 
lement destinée à mesurer les masses des particules 
du rayonnement cosmique; elle s'adresse à des 
particules électrisées traversant d’abord une chambre 
de Wilson sans écran, courbées par un champ magné- 
tique et se terminant enfin dans un des nombreux 
écrans d’une chambre de Wilson inférieure. La 
mesure, pour chaque trajectoire s’arrêtant dans un 
écrap, de la déviation par le champ magnétique et 
du parcours après passage dans ce champ, donne la 
masse de la particule. 

Alors que, dans la première expérience de Fretter 
avec 30 cm de plomb au-dessus de l'appareil, l’en- 
semble des résultats était cohérent avec une masse 
unique, à savoir celle du méson y, Brode, dans 
une expérience ultérieure sans filtre de plomb, 
observe 15 particules dont la masse est intermédiaire 
entre le méson 7 et le proton : les six dont la masse 
est le mieux déterminée sont positives. Un tel 
nombre correspondrait, dans ces conditions d’expé- 
riences, c’est-à-dire dans un domaine de parcours 
compris entre 1 et ro cm de plomb et au niveau 
du sol, à des proportions 

mésonT 


méson 4 


CARS) 


Pour vérifier ce fait, Peyrou et Lagarrigue ont 
utilisé sans plomb au-dessus de la chambre l'appareil 
décrit en 1948 [7] et fonctionnant à 1 000 m d’alti- 
tude à lArgentière-la-Bessée, dans les Alpes fran- 
çaises, appareil comportant une grande chambre de 
Wilson de 8o cm de hauteur avec champ magnétique 
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de 2 500 gauss, sous laquelle sont installés divers 
trains de compteurs et diverses plaques d’absorbants. 
Cette expérience est, dans son principe, très proche 
de celle de Brode. Les résultats détaillés sont décrits 
dans les Comptes rendus de l’Académie des Sciences [8]. 

La conclusion intéressant le méson + est que l’on. 
n’en a retrouvé aucun,sur 16 protons et sur {5 mésons 
mesurés. Il n’y a donc pas confirmation des résultats 
de Brode (1. Avec un tel dispositif expérimental, la 
présence de masses aberrantes ne peut être envisagée 
avec une certaine garantie que si leur groupement 
est convenable et leur proportion notable; mais le 
fait de trouver un résultat négatif montre que la 
proportion de mésons 7 est probablement inférieure 
à > pour 100 par rapport au méson ordinaire. 

Une autre expérience systématique a été faite par 
Franzinetti [9] à Bristol au moyen de plaques photo- 
graphiques placées l’une en face de l’autre à faible 
distance, un champ magnétique régnant dans l’espace 
séparant les deux plaques selon la technique de 
Powell et Rosenblum. Là encore, le résultat n’est 
pas positif : sur 30 mésons r et identifiés, aucune 
particule n’a été enregistrée comme étant avec certi- 
tude un méson 7, on retombe sur une proportion de 
mésons 7 inférieure à 3 pour 100 environ. 

Une troisième série d'expériences systématiques à 
été faite par le groupe de Bristol [10] (Fowler) en 
mesurant la masse des particules sortant des étoiles 
par l'observation de leur « scattering » dans l’émul- 
sion. Là encore, il n’y a pas de certitude en faveur 
du méson + : l’on trouve bien quelques points qui 
se placent sur la courbe de scattering des mésons 7 
ou -dans le voisinage de cette courbe, mais leur 
présentation n’est pas symétrique par rapport à la 
courbe; il se peut donc que ce soient des protons 
un peu aberrants. s 

Or, on peut développer certaines considérations sur 
ces mesures. Prenons les dernières nommées : la 
statistique est effectuée en prenant des particules 
provenant d'étoiles, assez longues pour permettre 
une bonne mesure du scattering. La plupart des 
impulsions mesurées par cette méthode sont infé- 
rieures à 700 MeV :c, 20 pour 100 seulement des 
trajectoires correspondent à des impulsions supé- 
rieures à cette valeur. 

Figurons les trois courbes donnant la densité de 
grains en fonction de l'inverse de l'angle moyen 
de scattering pour les mésons 7, les mésons 7 et les 
protons (fig. 1). 

Les particules dont les points représentatifs 
tombent entre À et D sont des mésons 7; de plus 
celles dont les points tombent entre D et B vont être 


() Dans une Communication récente (Bull. Am. Phys. 
Soc., 1949, 24, 8-15), Brode a annoncé qu'un examen plus 
critique des erreurs de la méthode rend suspecte l'interpré- 
tation des clichés de trajectoires considérées comme dues 
à des mésons +. Il y a donc accord entre les diverses expé- 
riences systématiques de mesurés de masse faites à la chambre 
de Wilson. 


N°5, 
aussi considérées comme des 7 si l’on exclut l’exis- 
tence de particules de plus faible masse. Si l'on 
envisageait l'existence de telles particules légères 


alors, la partie DB de la courbe serait suspecte et 
seule la partie AD compterait pour les mésons fr. 


ti) 


Fig. 
2, angle moyen de scattering; 
æ est proportionnel à p£: 


p, impulsion de la particule; 5 — 


Oz 


D'autre part, la courbe des mésons 7 n’est utili- 
sable, pour une discrimination, que jusqu’à F : à 
droite de F, on se trouve en compétition avec.les 
mésons ordinaires et aucune discrimination n’est 
possible. Dans cette région de moments, qui s’étend 
jusque vers 500 MeV : c, Fowler trouve 35 mésons 7 
pour o mésons 7. On peut en conclure que la propor- 
tion des 7 est inférieure à 3 pour 100 environ de 
celle des'z.. 

Enfin, la courbe pour les protons est utilisable de G 
en H. Elle sera nécessairement bien balisée car les 
protons observés dans les étoiles, qui ne sont pas 
créés dans la désintégration, ont un spectre d'énergie 
très étalé et l’on pourra accumuler, sur la portion GH, 
un grand nombre de résultats. 

On peut se demander si la comparaison des 7 et 
des 7 dans un même domaine de moments est signi- 
ficative et cela à cause de l'incertitude sur la loi de 
production des divers types de mésons. On pourrait 
imaginer également d'effectuer la comparaison, non 
plus dans un même domaine de moments, mais dans 
un même domaine de moments réduits : dans ce 
cas, les conclusions de Fowler sont à modifier dans 
un séns moins défavorable au méson 7. 

Mais on doit également considérer autre chose : 
comme les mésons sont produits dans les étoiles, et 
que les étoiles de moindre énergie sont les plus 
nombreuses, le nombre absolu des mésons + produits 
dans les étoiles doit être très inférieur à celui des 
mésons 7. Le nombre de phénomènes au cours 
desquels des mésons + peuvent être produits sera 
donc très inférieur au nombre de ceux susceptibles 
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de fournir des mésons 7 avec une probabilité du 
même ordre. Dans le même ordre d’idées remarquons 
que les mésons 7 détectés par Fowler ont le plus 
souvent une énergie inférieure à 150 MeV. On doit 
donc s'attendre à n’observer qu'une proportion 
faible de + par rapport aux . * 

L'ensemble des expériences systématiques et des 
considérations les accompagnant montre que l’on 
n'observe pas fréquemment des mésons 7, mais 
cela n’infirme pas l’interprétalion des cas particuliers 
cilés au début de cet article. 

Dans l'exemple de la photographie obtenue par 
Powell [4], il paraît très difficile d'éviter l’interpré- 
tation de la tripartition d’un méson 7 : en effet, le 
fait d'observer une particule dont la masse r000 
est bien déterminée par la mesure du scattering est 
déjà un phénomène étonnant car la trajectoire a 
plus de 2 500 2: et se termine au sein de l’émulsion, 
ce qui permet une bonne mesure. 

En second lieu, le fait d'observer une coïncidence 
parfaite entre l’extrémité de cette trajectoire et le 
centre d’une étoile est aussi étonnant, et le fait 
que cette étoile soit formée de trois. trajectoires 
coplanaires dont l’une est sans aucun doute un 
méson 7 négatif, et dont l’une des autres, mesurée à 
nouveau tout récemment par la méthode de scat- 
tering, est aussi un méson f, est également très 
étonnant. 

Des considérations analogues à celles qui viennent 
d’être développées pour le cliché de Powell peuvent 
aussi être données pour le cliché obtenu par notre 
groupe (?). 

Sur ce cliché, la longueur du méson + est nota- 
blement inférieure à celle du méson + de Powell 
mais, d’une part, la coïncidence spatiale est excellente 
et, d'autre part, l’on observe simultanément, dans 
le même phénomène, une trajectoire de proton et 
une trajectoire de méson 7, ce qui donne du poids 
au comptage de grains particulier du méson 7; 
l’argument fondamental de l'interprétation de ce 
cliché est d’ailleurs celui du bilan d’énergie. 

Les deux aspects de notre cliché et de celui de 
Powell apparaissent d’ailleurs comme ne présentant 
pas de contradiction. Ils peuvent même être consi- 
dérés comme complémentaires et correspondant, 
l’un à un méson 7 négatif donnant une étoile, l’autre 
à un méson 7 positif se décomposant en mésons 
plus légers, 

Il est toutefois très curieux que l'accumulation 
considérable des résultats n'ait pas permis de 
retrouver à plusieurs exemplaires les mêmes phéno- 
mènes, ce qui aurait dû se produire dans le cas où 
il s'agirait effectivement de diverses manifestations 
du méson +. Harding a présenté au Congrès d'Oxford 


() Il n’est pas question, dans cet article, de Ia discussion 
des résultats obtenus par Alichanov et collaborateurs, discus- 
sion qui sera faite par J. Daudin dans un.prochain numéro 
du Journal de Physique. 
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en mars 1950 deux cas de mésons 7 analogues à celui 
de Powell. 


Conclusion. —— La position qu'il paraît raison- 
nable d'adopter en regard du méson 7 nous semble 
être la suivante : 


19 Sélectivité des méthodes. — Les méthodes de 
séparation des particules doivent être examinées 
avec beaucoup de soin car chacune d'elles est, en 
général, sélective : ainsi la discussion indiquée plus 
haut montre que, dans l'expérience de mesure des 
masses par « scattering », la méthode risque de 
défavoriser les mésons 7 par rapport aux mésons t. 


20 Indépendamment de la sélectivité des méthodes.— 
a. Dans les mesures systématiques de masses faites 
à la chambre de Wilson (Brode, Fretter, Peyrou et 
Lagarrigue), la très grande majorité des particules 
qui traversent la chambre, et dont la masse est 
mesurée, provient de l’atmosphère et n’a pas été 
produite localement. En ce cas, seules sont enre- 
gistrées des particules dont la vie moyenné est assez 
longue. En effet, on ne parvient pas à enregistrer 
une proportion importante de mésons 7. 


b. Les mesures systématiques de masses faites 
avec deux plaques photographiques et un champ 
magnétique (Franzinetti) échappent à cette critique 
puisqu’un grand nombre de mésons 7, probablement 


produits pour la plupart dans la matière environ- 


nante, sont enregistrés. Cependant, ces expériences 
sont faites à 3 500 m, altitude à laquelle les phéno- 
mènes nucléaires sont surtout des phénomènes de 
faible énergie : d’ailleurs, la plupart des mésons 7 
produits à cette altitude -sont lents; il est possible 
que la production des mésons 7, même lents, qui 
exige beaucoup plus d'énergie, soit rare. - 


c. Dans l'expérience systématique de mesure des 
masses par scattering au départ du centre du 
phénomène nucléaire (Fowler), on peut supposer 
que le nombre des mésons r produits par gerbe 
est très inférieur au nombre des r produits par 
gerbe, quoique l’on ne sache rien sur les propriétés 
respectives de production des 7 et des x. 


On voit donc que l'absence d’une proportion 
notable des mésons 7 discriminés par les méthodes 
précédentes ne peut constituer un argument défi- 
nitif en faveur de la non-existence de cette particule. 

D'autre part, les phénomènes individuels rappelés 
plus haut ne peuvent actuellement, pour la plupart, 
être interprétés, ainsi qu'on l’a vu au cours de la 
discussion, qu’en admettant l'existence d’un méson x. 

Il est donc probable qu'une telle particule existe, 
mais sa discrimination échappe, en général, aux 
expériences systématiques, soit en raison du carac- 
tère sélectif des méthodes, soit en raison d’une faible 
abondance à la production, soit en raison d’une vie 
moyenne brève. 


Manuscril reçu le 25 janvier 1950. 
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SUR LA FOCALISATION DANS LES MONOCHROMATEURS A RÉSEAUX CONCAVES = 
DANS LE VIDE EN INCIDENCE NORMALE POUR L'ULTRAVIOLET 


Par Mne SimoNE ROBIN. 
Laboratoire d'Enseignement de Physique, Sorbonne. 


Sommaire. — La focalisation actuelle d’un monochromateur à vide à réseau concave en incidence 
normale est étudiée, indiquant ses avantages et ses inconvénients. Plusieurs autres systèmes de foca- 
lisation sont proposés pour ce même type d'appareil. : 


Un projet de monochromateur à réseau concave 
dans le vide pour l’ultraviolet lointain a été élaboré 
par M. Vodar vers 1943 et l'appareil construit par les 
Établissements Beaudouin est en service depuis quel- 
ques mois. Nous exposons ici les perfectionnements 
- suggérés par l’usage de cet appareil et complétant les 

procédés précédemment décrits [2, 5]. 

Au moment où notre monochromateur a été 
terminé, il était, à notre connaissance, le seul appareil 
de ce type, c’est-à-dire à réseau en incidence normale 
pour ultraviolet à vide (1). ‘Tout récemment 
Parkinson et Williams [4] ont décrit un appareil du 
même genre, et nous savons, d’autre part, qu’un 
autre appareil analogue est en construction aux 

_ États-Unis. Mais aucune discussion générale sur les 
procédés de focalisation n’a été publiée. 

. Cependant, pour les courtes longueurs d’onde 
(de 2 000 à 800 À par exemple) le réseau ne s’écarte 
jamais beaucoup de la normale au faisceau lumineux 
et cela entraîne certaines possibilités ou incompa- 
tibilités mécaniques. Nous avons donc cru utile 
étant donné les nombreuses recherches à l’ordre du 
jour auxquelles un monochromateur pour lultra- 


violet lointain est susceptible de servir (photo- 


conductivité, absorption, photoémission, etc.) d’ex- 
poser brièvement ici le cas.-particulier de focalisation 
d’un réseau, cas qui n'avait pas été considéré dans 
les publications précédentes [2, 5]. 


10 Il s’agit d’un système de focalisation approchée. 
Il a été utilisé sur notre monochromateur qui a été 
prévu à la fois pour être équipé d’un réseau concave 
ou d’un prisme de Féry. Il a l’avantage de s'adapter 
facilement à ces deux cas. 

Une seule manivelle, manœuvrée de l’extérieur, 
assure un mouvement de translation et un mouve- 
ment de rotation au système dispersif, les fentes 
d'entrée et de sortie étant fixes. 

Rappelons que, pour les prismes de Féry, la rota- 
tion se fait grâce à un bras de levier attaché à une 
extrémité à la plate-forme du prisme, à l’autre 


() Pour cette région spectrale, les autres monochromateurs 
décrits sont ceux de Hilsch et Pohl [3] à prismes et celui de 
Baker [1] à réseau tangent. 


extrémité, il s'appuie sur une glissière [5] (?). 

Pour le réseau, il est mécaniquement impossible, 
dans le cas des petites longueurs d'onde, d'utiliser 
ce procédé de focalisation géométriquement parfait, 
car le bras de levier se trouve presque perpendi- 
culaire à la glissière. Le constructeur a tourné la 
difficulté en attachant le même levier à une extrémité 
à la plate-forme du réseau, l’autre extrémité étant 
fixée en un point tel que le bras de levier soit perpen- 
diculaire à la glissière lorsque l’image centrale 
tombe sur la fente de sortie (bras qui est muni d’une 
coulisse de longueur variable). 

En faisant la mise au point pour 1550 À, on 
obtient avec le spectre du premier ordre pour un 
réseau de 57o -traits:mm et rayon R =:1 m, 
f — distance de la fente d’entrée au réseau, i — angle 
d'incidence, les résultats suivants (tableau I) : 


TABLEAU I. 
Focalisation 
) parfaite J réelle Arf 405. 

(49): (em). (em). moyen 
1200: 99,94 100, 16 DAT 
DOTE 99; 92 100,09 367 
ESS 99,90 99,90 () 
1066: 99,98 99,81 0,7 
FOOD 1 99; 86 7 


99,69 I; 


Le calcul montre que la condition de focalisation 
du réseau dans ce cas est Af — k(Ai)? et ici nous 
n’assurons qu'une correction de premier ordre par 
la relation f — kAï. La focalisation n’est acceptable 
que pour un domaine spectral étroit. La précision 
du pointé des longueurs d’onde est suffisante car la 
translation de la vis de commande de la rotation 
est relativement grande. 


29 On fixe le réseau’ pour une valeur moyenne 

- de 2 000 À de façon à assurer une focalisation rigou- 
reuse pour cette longueur d’onde et l’on fait tourner 

le réseau autour de cette position. C’est le procédé 

utilisé par Parkinson et Williams [4]. On trouve 


() Erratum. — Dans l'article côté [5], le point B (fig. 13) 
doit se trouver évidemment sur le cercle, 
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alors les résultats tableau II : 


TasLeau I. 

À M = 67, 
(À ). Jmoyeu 
TOPOS RENTE V2 
BOURSE 0 


© 


FO0O rene 


La précision du pointé des longueurs d'onde est 
indépendante de la focalisation et peut donc être 
suffisante. 


30 On. produit une rotation à l’aide d’une vis V, 
qui pousse et fait tourner en avançant un bras 
de levier solidaire de la plate-forme du réseau. En 
même temps, on produit une légère translation de 
la plate-forme du réseau à l’aide d’une vis V, de 
pas différent et solidaire de V.. Le déplacement 
de V, étant proportionnel à Ai, la correction n'est 
pas parfaite. 

7 l'on focalise exactement pour r000 et 4358 À; 

raleur réelle de / est, entre ces longueurs d’onde 


Fa — 9,1.107# À + 100,16 (en cm ). 


De plus, pour une avance donnée de la vis comman- 
dant la rotation du réseau, ce dernier subissant 
simultanément une légère translation, on a une erreur 
d'incidence i due à cette translation. Mais comme 
celle-ci est faible, cette erreur est presque négligeable. 
Tout compte fait, on arrive au tableau IF : 


Tagceau II. 


k f réelle Ar «103, 

(CA (cm). / moyen 
TOO 99:92 (e] 
DO00: 99,7? ,8 
DO0O AR N 99; 54 0,9 
4000...:. 99:33 O,1 
HI DONNE 99, 22 o 


Dans ce procédé, il est diflicile de réaliser un bras 
de levier assez long pour obtenir une bonne précision 
pour le repérage des longueurs d’onde. 


4° Le réseau est sur une plate-forme munie d’un 
bras de levier L. Cette plate-forme est solidaire d’une 
came s'appuyant sur une surface fixe. Le bras de 
levier L s'engage (entre deux rouleaux) dans une 
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rainure d’une Diatetutae p' Fe subit une Trans 
lation +(cm). Cette translation entraîne la ‘rotation 
du réseau et aussi celle de la came solidaire de la. 
plate-forme P. La came peut être calculée pour. 
assurer une focalisation pratiquement parfaite. 
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X cm | 
Lo 
Ü 
Fig. La 


Quant à la précision du repérage de la rotation 
elle peut être extrêmement grande car la trans-, - 
lation æ de P' est considérable et d’ailleurs sensi- 
blement linéaire en fonction de la longueur d'onde. 
Si l'angle du bras L avec la normale da réseau est 
‘de 50, on aboutit au tableau IV : 


TABLEAU IV. 


(À). + 5% æ (cm). 

TOOO NES 6038" 733 b3 

DOOO EE 8 16 78,32 

SO00n ae 054 81,32 | 
f000..::. 11 33 83,29 

4358.24 12 8 83,80 


On constate que la translation entre 1000 
et 4358 À est de 10,27 Cm. PTE 
Manuscrit reçu,le 27 février 1950. 
} ; 
/ 
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